working example, too high values due to missing h^2 in rhs function