// INFLUENCIA DE LAS CONDICIONES INICIALES - SUBESPACIOS "ESTABLE" E "INESTABLE"

// SISTEMA NO ESTABLE

// x1' = 2 x1 + x2
// x2' = 2 x1 - x2

A = [2 1 ; 2 -1]
// R matriz de vectores propios y L valores propios en la diagonal
[R,L] = spec(A)

// subespacio "ESTABLE", asociado al vector propio [0.2703;-0.9628]
t=0:0.1:4;
x=R(:,2)*exp(L(2,2)*t);

clf
subplot(211)
plot2d(t',x')
xlabel('t')
ylabel('x')
title('subespacio ESTABLE')
legend('x1','x2')

// subespacio "INESTABLE", asociado al vector propio [0.8719;0.4896]
t=0:0.05:1;
x=R(:,1)*exp(L(1,1)*t);

subplot(212)
plot2d(t',x')
xlabel('t')
ylabel('x')
title('subespacio INESTABLE')
legend('x1','x2')
