% programa para chamada o metodo de integracao
% condicao incial
x0 = 1;

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

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

%[t,x] = Euler(t0,tf,x0,'f1',n);
%[t,x] = EulerMOD(t0,tf,x0,'f1',n);
%[t,x] = EulerAP(t0,tf,x0,'f1',n);
%[t,x] = RK(t0,tf,x0,'f1',n);
[t,x] = Adams(t0,tf,x0,'f1',n);

% plotando a solucao exata
tt = t0:0.05:tf;
xx = feval('f1_exata',tt);


%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)
plot(t,x,'o;dados;',tt,xx,'r-;exata;')

%fixando o grafico
figure(2);
%definindo detalhes
title('erro da integracao')
xlabel('t')
ylabel('erro(t)')
%plotando 
xe = feval('f1_exata',t);
xe = x-xe;
plot(t,xe,'r-;erro;')




