% 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;

%intervalo de integracao
t0 = 0;
tf = 5;

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

fun = 'f3';

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(3);
%definindo detalhes
title('resultado da integracao')
xlabel('t')
ylabel('x(t)')
%plotando os pontos (circulos) e a funcao (linha vermelha)
text
plot(t,x,'o;text;')

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