% metodo de Euler simples (explicito)
function [t,x] = Euler(t0,tf,x0,dxdt,n)

% atribui o valor de n = 100 como padrao
if nargin < 5, n = 50; end

% armazenado as condicoes iniciais
t(1) = t0;
x(1) = x0;

% calculo do passo
h = (tf-t0)/n;

% loop de integracao
for i = 1:n
     dx = feval(dxdt,t(i),x(i));
     x(i+1) = x(i)+h*dx;
     t(i+1) = t(i)+h;
endfor

end