clc; clear all; close all;
% 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] = eig(A);

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

subplot(211)
plot(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)
plot(t',x')
xlabel('t')
ylabel('x')
title('subespacio INESTABLE')
legend('x1','x2')