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