f(x, y) = sin((x * y))
numeriskIntegral1 = SolveODE(f, 0, 3, 5, 0.1)