ti 10.10.00
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 shgAjon tuloksia