A parte il problema con le parentesi (vedi risposta precedente - nell'equazione che hai scritto, ne manca qualcuna), si può impostare il ciclo "for" come segue:
clear all;
clc;
m=0.65; %massa sistema [kg]
k=8e7; %rigidezza [N/m]
zeta=0.005; %rapporto di smorzamento
x0=0.05; %spostamento iniziale [m]
v0=-2.8; %velocità iniziale [m/s]
omn=sqrt(k/m)
oms=omn*sqrt(1-zeta^2)
A=x0 %costante di integrazione
B=(v0+zeta*omn*x0)/(oms) %costante di integrazione
t=[0 2 3 5]
for i=1:length(t)
x(1,i)=exp(-zeta*omn*t(1,i)*[A*cos(oms*t(1,i)+B*sin(oms*t(1,i)];
end
plot(t,x)
hold on
Il codice riportato sopra dovrebbe darti un messaggi di errore, in quanto come già detto, mancano delle parentesi in
x(1,i)=exp(-zeta*omn*t(1,i)*[A*cos(oms*t(1,i)+B*sin(oms*t(1,i)];
Controlla l'equazione e verifica dove aggiungerle.
Volendo, si potrebbe ottenere lo stesso risultato senza usare il ciclo "for" sfruttando la proprietà di "indxing" di MatLab, ma, per cominciare, usare un ciclo "for" fa lo stesso.