This is an ODE of the form y' + P(x)y = Q(x)
IF = [tex]e^{ \int\limits {P(x)} \, dx }=e^{ \int\limits {2x} \, dx }=e^{x^2}[/tex]
Multiply through by IF
[tex]e^{x^2}y'+2xe^{x^2}y=2x^3e^{x^2} \\ (e^{x^2}y)'=2x^3e^{x^2} \\ \int {(e^{x^2}y)'} \, dx = \int {2x^3e^{x^2}} \, dx \\ e^{x^2}y=x^2e^{x^2}-e^{x^2} + c \\ y=c_1e^{-x^2}+x^2-1[/tex]
But y(0) = 1
[tex]1=c_1e^{-(0)^2}+(0)^2-1 = c_1-1 \\ c_1=1+1=2[/tex]
Therefore, solution is
[tex]y=2e^{-x^2}+x^2-1[/tex]
