solve by MAPLE
first, solve 2nd order differential equation
ode2 := x^2*(diff(f(x), x, x))-2*x*(diff(f(x), x))+2*f(x) = 0
ics2 := (D(f))(1) = 9, f(1) = 4
dsolve([ode2, ics2])
answer is f(x) = 5*x^2-x
convert to a system of 1st ode
sys1ode := diff(y(t), t) = z(t), diff(z(t), t) =...