f = @(x,y) y;
x0 = 0;
y0 = 1;
h = 1;
actual_y = forward_euler1(x0, y0, f, h);
expected_y = 2;

unit_test_equals("Forward Euler works for one step", ...
		 expected_y, ...
		 actual_y);
