ODE1[f_, x0_List, delta_, n_Integer] :=
	Module[{xi = x0, res = {x0}},
		Do[ xi = xi + delta f[xi]; AppendTo[res, xi],
		   {n} ];
		res
	]
