Onestep Theta Solver for the ODE g1(x)y + g2(x)y' = f(x)
x0
xn
y0
h
theta
f(x)
g1(x)
g2(x)