matlab dsolve('D2y=(1/x)*Dy+x*e^x','x') ans = 1/2*x^2*C1+(-1+x*log(e))/log(e)^3*e^x+C2
y=(x-1)*exp(x)+A*x^2+B
