forward_euler1.m

   1 function y = forward_euler1(x0, y0, f, h)
   2   ## Compute one iteration of the forward Euler method.
   3   ##
   4   ## INPUT:
   5   ##
   6   ##   * ``x0`` - The initial x-coordinate.
   7   ##
   8   ##   * ``y0`` - The initial value y(x0).
   9   ##
  10   ##   * ``f`` - The function y' = f(x,y) of two variables.
  11   ##
  12   ##   * ``h`` - The step size.
  13   ##
  14   ## OUTPUT:
  15   ##
  16   ## The approximate value of y(x) at x = x0+h.
  17   ##
  18   y_prime = f(x0,y0);
  19   y = y0 +  h * y_prime;
  20 end