desolve(x*(x^2-1)*y'+2*y=0,y)
(c_0*x^2)/(x^2-1)
desolve(x*(x^2-1)*y'+2*y=x^2,y)
((ln(x)+c_0)*x^2)/(x^2-1)