Add is_positive_definite() and its tests.
[octave.git] / 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