% Harj6 LV, pe 26.10.01

% Teht. 2
A=[-14 -160 -40;
    181 5   2;
     96 84  18]

[V,D]=eig(A)
lamc=D(1,1)
w=V(:,1)
     u=real(w),v=imag(w)
     lamr=D(3,3)
     wr=V(:,3)
clf
     plot3([0,u(1)],[0,u(2)],[0,u(3)],'b')
     hold on
plot3([0,v(1)],[0,v(2)],[0,v(3)],'b')
plot3([0,wr(1)],[0,wr(2)],[0,wr(3)],'r')
     grid; shg

% Teht. 3
A=[0 0 1 0;0 0 0 1;-5 2 0 0;2 -2 0 0]

[V,D]=eig(A)
x0=[1; 0; 1; 0]
[t,y]=linsys(A,x0,2);

plot(t,y(:,1:2))
plot(y(:,1),y(:,2))
grid

[t,y]=linsys(A,x0,20,200);

clf
plot(t,y(:,1:2))
plot(y(:,1),y(:,2))
grid
clf
plot(y(1:20,1),y(1:20,2))
hold on
pause
plot(y(20:40,1),y(20:40,2))
grid



plot(y(1:20,1),y(1:20,2)
plot(y(1:20,1),y(1:20,2)
plot(y(1:20,1),y(1:20,2)



