% variabili
clear all;
q=0.3;
b=2;
w=4;
K=b*w/(q^2+w^2);
K1=(q^2+w^2)^0.5/w;
fi=-atan(w/q);
t=0:0.001:10;


x=K*(exp(-q*t)+K1*sin(w*t+fi));
plot(t,x);