令z=y^2
dz/dx=2y(dy/dx)=2yy'
所以原方程變?yōu)?br/>z'-xz=xe^x
z(0)=y(0)^2=1
然后利用積分因子
e^(∫-xdx)
=e^(-x^2/2)
兩邊同乘,左邊是一全微分
(ze^(-x^2/2))'=xe^(x-x^2/2)
兩邊積分
ze^(-x^2/2)=∫ te^(t-t^2/2) dt +C
z=e^(x^2/2)[∫ te^(t-t^2/2) dt +C]
代入x=0,z=1
1=C
所以
z=e^(x^2/2)[∫ te^(t-t^2/2) dt +1]
y=e^(x^2/4)[∫ te^(t-t^2/2) dt +1]^(1/2)