function [t , y] = mieuler( funcion, t0, x0, tf, N ) % funcion mieuler(funcion, t0, x0, xf, N ) % % resuelve utilizando el metodo de Euler la ecuacion diferencial % x' = funcion(t,x) % funcion debe ser el nombre de una funcion .m % t0 debe ser el tiempo inicial % x0 valor inicial x(t0) % tf valor final para t % N numero de pasos h = (tf-t0) / n; t = [ t0 : h : tf ]'; y = zeros(size(t)); y(1) = y0; for i = 1 : n ff = feval(funcion, x(i), y(i)); y(i+1) = y(i) + h * ff; end