clear all
clc
close all

iteraciones = 100;
f = zeros(iteraciones, 1);
deltat = 0.1;


f(1)=1;
n = 1;
while n <= iteraciones - 1
  f(n+1) = f(n) - (f(n) * deltat)
  f(n) * deltat  + f(n)
  n++;
end


figure(1, "visible", "off");
plot(f);
print -deps "leibniz.eps";