% script to compute the data to solve assignment three
diary sol_three_data
% first part
figure
f = inline('[0.08*y(1) + 1.44*y(2); 1.44*y(1) + 0.92*y(2)]','t','y');
xzero = [-4 3];
tspan = [0 10];
[t x] = ode45(f,tspan,xzero);
endx = x(length(x),:)
plot(t,x)
% second part: change initial conditions
figure
xzero1 = [-4.001  3.001];
[t1 x1] = ode45(f,tspan,xzero1);
endx1 = x1(length(x1),:)
hold on
plot(t1,x1);
% third part: interpret
A = [0.08 1.44; 1.44 0.92];
[v,d] = eig(A)
c = v\xzero'
diary off