clc
clear all
close all

%% Entradas para a simulação numérica
y0=[0;0.1]; %vetor de condições iniciais, suposto para um sistema de dimensão 2
A=[0 1;-1 -0.1]; %Exemplo para um sistema de dimensão 2
dt=0.01; %intervalo de tempo nos quais os estados serão avaliados na integração numérica
tf=1000; %tempo final da simulação
t=[0:dt:tf];

%% Integração numérica
f=@(t,y) Integrador(t,y,A);
[tsim,xsim]=ode45(f,t,y0);
%xsim é um matriz com Nt linhas e, no caso, duas colunas. A série temporal
%y1(t) está na primeira coluna da matriz xsim. Já a segunda coluna traz a
%série temporal y2(t)
%tsim é o vetor de tempo salvo. No caso, coincide com t.

%% Exibição de gráficos
plot(tsim,xsim(:,1),'k')
xlabel('t','fontsize',14)
ylabel('y1(t)','fontsize',14)

%% Cálculo dos auto-valores da matriz A
[V,D]=eig(A); % V é uma matriz cuja colunas correspondem aos auto-vetores de A. 
%D é uma matriz diagonal com os correspondentes auto-valores