clear all

A=[-1 10;0 -2];

x1=[1;0];
x2=[1;-1/10];x2=x2/norm(x2);

a1=[1/10;1];a1=a1/norm(a1);
a2=[0;1];

theta=1/10;


dt=0.1;
Prop=expm(A*dt);

Nt=200;
t=zeros(Nt,1);E=t;x=t;y=t;

xt=a2;
E(1)=xt'*xt;
x(1)=xt(1);y(1)=xt(2);

for it=2:Nt;
    t(it)=(it-1)*dt;
    
    xt=Prop*xt;
    E(it)=xt'*xt;
    x(it)=xt(1);y(it)=xt(2);
end;

figure(34);hold on
plot(t,E,'g','Linewidth',2)
xlabel('$t$','Interpreter','latex','Fontsize',24)
ylabel('$E(t)$','Interpreter','latex','Fontsize',24)

figure(35);hold on
plot(x,y,'g','Linewidth',2)
xlabel('$x$','Interpreter','latex','Fontsize',24)
ylabel('$y$','Interpreter','latex','Fontsize',24)








