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

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

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

% calculo do passo
h = (tf-t0)/n;
k = 1000;
% loop de integracao
for i = 1:n
     xs     = xs/(1+h*(1+k));
     t(i+1) = t(i)+h;
     x = vertcat(x,xs);
endfor

end