b = 0.1 (slider from 0 to 1, step 0.1)
c = 0.9 (slider from 0 to 10, step 0.1)
ω: sqrt(4 * c - b^(2)) / 2
g(x) = 1 / ω * exp(-b * x / 2)
h(x) = -g(x)
E2 = SolveODE(b, c, sin(x), 0, 0, 1, 100, 0.1)
E = SolveODE(b, c, 0, 0, 0, 1, 100, 0.1)