Internal
problem
ID
[9996]
Book
:
Collection
of
Kovacic
problems
Section
:
section
2.
Solution
found
using
all
possible
Kovacic
cases
Problem
number
:
3
Date
solved
:
Wednesday, March 05, 2025 at 08:01:57 AM
CAS
classification
:
[[_Emden, _Fowler]]
Solve
Time used: 0.092 (sec)
Writing the ode as
Comparing (1) and (2) shows that
Applying the Liouville transformation on the dependent variable gives
Then (2) becomes
Where
Substituting the values of
Comparing the above to (5) shows that
Therefore eq. (4) becomes
Equation (7) is now solved. After finding
The first step is to determine the case of Kovacic algorithm this ode belongs to. There are 3 cases depending on the order of poles of
Case |
Allowed pole order for |
Allowed value for |
1 |
|
|
2 |
Need to have at least one pole
that is either order |
no condition |
3 |
|
|
The order of
The poles of
Attempting to find a solution using case
Looking at poles of order 2. The partial fractions decomposition of
For the pole at
Since the order of
Since the
The following table summarizes the findings so far for poles and for the order of
pole |
pole order | |
|
|
| | | | |
Order of |
|
|
|
|
| | |
Now that the all
Where
Since
The above gives
Now that
Let
Substituting the above in eq. (1A) gives
The equation is satisfied since both sides are zero. Therefore the first solution to the ode
The first solution to the original ode in
Since
Which simplifies to
The second solution
Since
Therefore the solution is
Will add steps showing solving for IC soon.
ode:=diff(diff(y(x),x),x) = 12*y(x)/x^2; dsolve(ode,y(x), singsol=all);
Maple trace
`Methods for second order ODEs: --- Trying classification methods --- trying a quadrature checking if the LODE has constant coefficients checking if the LODE is of Euler type <- LODE of Euler type successful`
Maple step by step
ode=D[y[x],{x,2}]==((4*(7/2)^2-1)/(4*x^2))*y[x]; ic={}; DSolve[{ode,ic},y[x],x,IncludeSingularSolutions->True]
from sympy import * x = symbols("x") y = Function("y") ode = Eq(Derivative(y(x), (x, 2)) - 12*y(x)/x**2,0) ics = {} dsolve(ode,func=y(x),ics=ics)