j=1;   h=0.02;   

a=[1;2;1];   b=[-sqrt(3);0;sqrt(3)];   c=[sqrt(2);-sqrt(2);sqrt(2)];

for t=0:h:50
    x=-0.088*a*(1/1)*sin(3*t)-2.63*b*(1/1)*sin(3*t)+0.21*c*(1/1)*sin(3*t);
    u1(j)=x(1);   u2(j)=x(2);   u3(j)=x(3);
    j=j+1;
end

t=0:h:50;