clear all;





ksi = [0.5,-0.5,1,0.001,-0.001,22];
omega = [3,-3,5,0.1,3,3];
a= 1;
d=1;
t0=0;
du = 10;
u0=0;
t1=1/du;



x10 = 0;
figure(3);

for i=1:6
    c = omega(i)*omega(i);
    b= 2*ksi(i)*omega(1);
    x0= (b*u0)/(c);
    [t]=sim('lab_5_sim2');
    subplot(6,1,i);
    plot(t,x);
end;
