dsolve(x*(x^2-1)*diff(y(x),x)+(x^2-1)*y(x)^2-x^2 = 0,y(x))