% programa para chamada o metodo de integracao
disp('Metodos disponiveis:')
disp('1 - Euler')
disp('2 - Euler Modificado')
disp('3 - Euler Aprimorado')
disp('4 - Runge Kutta')
disp('5 - Adams')
Metodo = input("Escolha um metodo: ");

% condicao incial
x0 = [1 0 0];
%intervalo de integracao
t0 = 0;
tf = 1;

% numero de pontos
n = 20000;
dt = (tf-t0)/n;
%chamada do integrador

fun = 'Ex3';

switch(Metodo)
case 1
[t,x] = Euler(t0,tf,x0,fun,n); text="Euler";
case 2
[t,x] = EulerMOD(t0,tf,x0,fun,n); text='Euler MOD';
case 3
[t,x] = EulerAP(t0,tf,x0,fun,n); text='Euler AP';
case 4
[t,x] = RK(t0,tf,x0,fun,n); text='RK';
case 5
[t,x] = Adams(t0,tf,x0,fun,n); text='Adams';
endswitch


%fixando o grafico
%figure(1);
%definindo detalhes
title('resultado da integracao')
xlabel('t')
ylabel('x(t)')
%plotando os pontos (circulos) e a funcao (linha vermelha)
text

x1 = x(:,1);
x2 = x(:,2);
x3 = x(:,3);

plot(t,x1,'o-;x1;',t,x2,'o-;x2;',t,x3,'o-;x3;')

%filename = "resultados.txt";
%fid = fopen (filename, "w");
save -ascii resultados.txt t x;
%fclose (fid);
