function [x , y] = mieuler( funcion, x0, y0, b, n )
% funcion  mieuler(funcion, x0, y0, b, n )
%
% resuelve utilizando el metodo de Euler la ecuacion diferencial
%   y' = funcion(x,y)
% funcion debe ser el nombre de una funcion .m
% x0 debe ser el punto inicial
% y0 valor inicial y(x0)
% b  valor final para x
% n  numero de pasos

h = (b-x0) / n;

x = [ x0 : h : b ]';

y = zeros(size(x));

y(1) = y0;

for i = 1 : n
  ff = feval(funcion, x(i), y(i));
  y(i+1) = y(i) + h * ff;
end