--- /dev/null
+function y = forward_euler1(x0, y0, f, h)
+ ## Compute one iteration of the forward Euler method.
+ ##
+ ## INPUT:
+ ##
+ ## * ``x0`` - The initial x-coordinate.
+ ##
+ ## * ``y0`` - The initial value y(x0).
+ ##
+ ## * ``f`` - The function y' = f(x,y) of two variables.
+ ##
+ ## * ``h`` - The step size.
+ ##
+ ## OUTPUT:
+ ##
+ ## The approximate value of y(x) at x = x0+h.
+ ##
+ y_prime = f(x0,y0);
+ y = y0 + h * y_prime;
+end