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.
- ##
+ % 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