Luento 13

ti 10.10.00 

Matlab-skripti trajektorien piirtoon

format compact
A=[1 1;4 1]
[V,D]=eig(A),lam=diag(D)
t=linspace(-10,0.3); 
% t=logspace(-1,0.1); % joissakin tapauksissa kannattaa kokeilla.
clf
axis([-5 5 -5 5]);hold on
varit=['br']
i=0;
for C1=-4:4
   i=i+1;
   for C2=-3:3
     X=V*[C1*exp(lam(1)*t);C2*exp(lam(2)*t)];
     plot(X(1,:),X(2,:),varit(1+mod(i,2)))
   end
end
grid
shg