tests/forward_euler1_tests.m

   1 f = @(x,y) y;
   2 x0 = 0;
   3 y0 = 1;
   4 h = 1;
   5 actual_y = forward_euler1(x0, y0, f, h);
   6 expected_y = 2;
   7
   8 unit_test_equals("Forward Euler works for one step", ...
   9                  expected_y, ...
  10                  actual_y);