X-Git-Url: http://gitweb.michael.orlitzky.com/?p=octave.git;a=blobdiff_plain;f=forward_euler1.m;fp=forward_euler1.m;h=3e1061d48340b1b3db9c4c4c45c1ced56528b82e;hp=0000000000000000000000000000000000000000;hb=24b6a4fe081d46f6a2f7ee2a3161401256d21fb5;hpb=b075adde3bee54eaded3dd161f037e3a8de88618

diff --git a/forward_euler1.m b/forward_euler1.m
new file mode 100644
index 0000000..3e1061d
--- /dev/null
+++ b/forward_euler1.m
@@ -0,0 +1,20 @@
+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