By the way, this is the ODE `eq = x**2*(a*f(x)**2+(f(x).diff(x)))+b*x**alpha+c`