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