Exercises/Experiments

1.
A=[0 0 -2; 1 2 1; 1 0 3]

1a.
[evect, eval]=eig(A)
J=inv(evect)*A*evect

1b.
P=inv(evect)
J=P*A*inv(P)

2.
C=[-1 1 0; 1 -1 0; 0 0 -1]

2a.
[evects, evals] = eig(C)

2b.
% contents of the m-file exer252.m
%   function Zprime=exer252(t,Z);
%   Zprime=[-1*Z(1)+1*Z(2); 1*Z(1)-1*Z(2); -1*Z(3)];
[t,Z]=ode23('exer252',[0 10],[1; 0; 1]);
plot(t,Z)