dxdy=ye−x2
ydy=e−x2dx
∫ydu=∫e−x2dx
But ∫e−x2dx is a non elementary integral. We can express it as a power series
e−x2=1−1!x2+2!x4−3!x6+4!x8−..
We can approximate it as
e−x2=1−1!x2+2!x4−3!x6
So differential equation transforms to
∫ydu=∫ (1−1!x2+2!x4−3!x6)dx
=>ln|y| = x−3x3+10x5−42x7+C
By initial condition y(5)=1
0 = 5−353+1055−4257+C
C = 1584.286
So ln|y| = x−3x3+10x5−42x7+1584.286
=> y = e(x−3x3+10x5−42x7+1584.286)
This is the explicit solution of the given differential equation
ALTERNATIVELY
dxdy=ye−x2
=> ydy=e−x2dx
=> ∫ydu=∫e−x2dx
As ∫e−x2dx is not an elementary integral, we will express it as error function
=> ln|y| = (2π ) erf(x) + C1, C1 is integration constant
Where erf(x) = ∫0xe−t2dt
So y = e(2π)erf(x)+C1
=> y = Ce(2π)erf(x)
By initial condition y (5)=1
1 = Ce(2π)erf(5)
C = e−(2π)erf(x)
So y = e(2π)erf(x)−(2π)erf(5)
=> y = e(2π)[erf(x)−erf(5)]
This is the explicit solution of the given differential equation
Comments