2.1.318 Problem 325

Solved as second order ode using Kovacic algorithm
Maple
Mathematica
Sympy

Internal problem ID [9490]
Book : Collection of Kovacic problems
Section : section 1
Problem number : 325
Date solved : Wednesday, March 05, 2025 at 07:49:59 AM
CAS classification : [[_2nd_order, _with_linear_symmetries]]

Solve

x2y+(53x+x2)yy3=0

Solved as second order ode using Kovacic algorithm

Time used: 0.418 (sec)

Writing the ode as

(1)x2y+(53x+x2)yy3=0(2)Ay+By+Cy=0

Comparing (1) and (2) shows that

A=x2(3)B=53x+x2C=13

Applying the Liouville transformation on the dependent variable gives

z(x)=yeB2Adx

Then (2) becomes

(4)z(x)=rz(x)

Where r is given by

(5)r=st=2AB2BA+B24AC4A2

Substituting the values of A,B,C from (3) in the above and simplifying gives

(6)r=9x2+30x+736x2

Comparing the above to (5) shows that

s=9x2+30x+7t=36x2

Therefore eq. (4) becomes

(7)z(x)=(9x2+30x+736x2)z(x)

Equation (7) is now solved. After finding z(x) then y is found using the inverse transformation

y=z(x)eB2Adx

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 r and the order of r at . The following table summarizes these cases.

Case

Allowed pole order for r

Allowed value for O()

1

{0,1,2,4,6,8,}

{,6,4,2,0,2,3,4,5,6,}

2

Need to have at least one pole that is either order 2 or odd order greater than 2. Any other pole order is allowed as long as the above condition is satisfied. Hence the following set of pole orders are all allowed. {1,2},{1,3},{2},{3},{3,4},{1,2,5}.

no condition

3

{1,2}

{2,3,4,5,6,7,}

Table 2.318: Necessary conditions for each Kovacic case

The order of r at is the degree of t minus the degree of s. Therefore

O()=deg(t)deg(s)=22=0

The poles of r in eq. (7) and the order of each pole are determined by solving for the roots of t=36x2. There is a pole at x=0 of order 2. Since there is no odd order pole larger than 2 and the order at is 0 then the necessary conditions for case one are met. Since there is a pole of order 2 then necessary conditions for case two are met. Therefore

L=[1,2]

Attempting to find a solution using case n=1.

Looking at poles of order 2. The partial fractions decomposition of r is

r=14+736x2+56x

For the pole at x=0 let b be the coefficient of 1x2 in the partial fractions decomposition of r given above. Therefore b=736. Hence

[r]c=0αc+=12+1+4b=76αc=121+4b=16

Since the order of r at is Or()=0 then

v=Or()2=02=0

[r] is the sum of terms involving xi for 0iv in the Laurent series for r at . Therefore

[r]=i=0vaixi(8)=i=00aixi

Let a be the coefficient of xv=x0 in the above sum. The Laurent series of r at is

(9)r12+56x12x2+56x35936x4+385108x52681324x6+19525972x7+

Comparing Eq. (9) with Eq. (8) shows that

a=12

From Eq. (9) the sum up to v=0 gives

[r]=i=00aixi(10)=12

Now we need to find b, where b be the coefficient of xv1=x1=1x in r minus the coefficient of same term but in ([r])2 where [r] was found above in Eq (10). Hence

([r])2=14

This shows that the coefficient of 1x in the above is 0. Now we need to find the coefficient of 1x in r. How this is done depends on if v=0 or not. Since v=0 then starting from r=st and doing long division in the form

r=Q+Rt

Where Q is the quotient and R is the remainder. Then the coefficient of 1x in r will be the coefficient in R of the term in x of degree of t minus one, divided by the leading coefficient in t. Doing long division gives

r=st=9x2+30x+736x2=Q+R36x2=(14)+(30x+736x2)=14+30x+736x2

Since the degree of t is 2, then we see that the coefficient of the term x in the remainder R is 30. Dividing this by leading coefficient in t which is 36 gives 56. Now b can be found.

b=(56)(0)=56

Hence

[r]=12α+=12(bav)=12(56120)=56α=12(bav)=12(56120)=56

The following table summarizes the findings so far for poles and for the order of r at where r is

r=9x2+30x+736x2

pole c location pole order [r]c αc+ αc
0 2 0 76 16

Order of r at [r] α+ α
0 12 56 56

Now that the all [r]c and its associated αc± have been determined for all the poles in the set Γ and [r] and its associated α± have also been found, the next step is to determine possible non negative integer d from these using

d=αs()cΓαcs(c)

Where s(c) is either + or and s() is the sign of α±. This is done by trial over all set of families s=(s(c))cΓ until such d is found to work in finding candidate ω. Trying α+=56 then

d=α+(αc1)=56(16)=1

Since d an integer and d0 then it can be used to find ω using

ω=cΓ(s(c)[r]c+αcs(c)xc)+s()[r]

Substituting the above values in the above results in

ω=(()[r]c1+αc1xc1)+(+)[r]=16x+(12)=16x+12=16x+12

Now that ω is determined, the next step is find a corresponding minimal polynomial p(x) of degree d=1 to solve the ode. The polynomial p(x) needs to satisfy the equation

(1A)p+2ωp+(ω+ω2r)p=0

Let

(2A)p(x)=x+a0

Substituting the above in eq. (1A) gives

(0)+2(16x+12)(1)+((16x2)+(16x+12)2(9x2+30x+736x2))=013a03x=0

Solving for the coefficients ai in the above using method of undetermined coefficients gives

{a0=13}

Substituting these coefficients in p(x) in eq. (2A) results in

p(x)=x13

Therefore the first solution to the ode z=rz is

z1(x)=peωdx=(x13)e(16x+12)dx=(x13)ex2ln(x)6=(1+3x)ex23x1/6

The first solution to the original ode in y is found from

y1=z1e12BAdx=z1e1253x+x2x2dx=z1ex25ln(x)6=z1(ex2x5/6)

Which simplifies to

y1=1+3x3x

The second solution y2 to the original ode is found using reduction of order

y2=y1eBAdxy12dx

Substituting gives

y2=y1e53x+x2x2dx(y1)2dx=y1ex5ln(x)3(y1)2dx=y1(9ex5ln(x)3x2(1+3x)2dx)

Therefore the solution is

y=c1y1+c2y2=c1(1+3x3x)+c2(1+3x3x(9ex5ln(x)3x2(1+3x)2dx))

Will add steps showing solving for IC soon.

Maple. Time used: 0.049 (sec). Leaf size: 29
ode:=x^2*diff(diff(y(x),x),x)+(5/3*x+x^2)*diff(y(x),x)-1/3*y(x) = 0; 
dsolve(ode,y(x), singsol=all);
 
y=c1x4/3hypergeom([2],[73],x)ex3c2x+c2x

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 
trying a symmetry of the form [xi=0, eta=F(x)] 
checking if the LODE is missing y 
-> Trying a Liouvillian solution using Kovacics algorithm 
   A Liouvillian solution exists 
   Reducible group (found an exponential solution) 
   Group is reducible, not completely reducible 
   Solution has integrals. Trying a special function solution free of integrals... 
   -> Trying a solution in terms of special functions: 
      -> Bessel 
      -> elliptic 
      -> Legendre 
      <- Kummer successful 
   <- special function solution successful 
      Solution using Kummer functions still has integrals. Trying a hypergeometric solution. 
      -> hyper3: Equivalence to 2F1, 1F1 or 0F1 under a power @ Moebius 
      <- hyper3 successful: received ODE is equivalent to the 1F1 ODE 
      -> Trying to convert hypergeometric functions to elementary form... 
      <- elementary form is not straightforward to achieve - returning special function solution free of uncomputed integrals 
   <- Kovacics algorithm successful`
 

Maple step by step

Let’s solvex2(d2dx2y(x))+(53x+x2)(ddxy(x))y(x)3=0Highest derivative means the order of the ODE is2d2dx2y(x)Isolate 2nd derivatived2dx2y(x)=y(x)3x2(5+3x)(ddxy(x))3xGroup terms withy(x)on the lhs of the ODE and the rest on the rhs of the ODE; ODE is lineard2dx2y(x)+(5+3x)(ddxy(x))3xy(x)3x2=0Check to see ifx0=0is a regular singular pointDefine functions[P2(x)=5+3x3x,P3(x)=13x2]xP2(x)is analytic atx=0(xP2(x))|x=0=53x2P3(x)is analytic atx=0(x2P3(x))|x=0=13x=0is a regular singular pointCheck to see ifx0=0is a regular singular pointx0=0Multiply by denominators3x2(d2dx2y(x))+x(5+3x)(ddxy(x))y(x)=0Assume series solution fory(x)y(x)=k=0akxk+rRewrite ODE with series expansionsConvertxm(ddxy(x))to series expansion form=1..2xm(ddxy(x))=k=0ak(k+r)xk+r1+mShift index usingk>k+1mxm(ddxy(x))=k=1+mak+1m(k+1m+r)xk+rConvertx2(d2dx2y(x))to series expansionx2(d2dx2y(x))=k=0ak(k+r)(k+r1)xk+rRewrite ODE with series expansionsa0(1+r)(1+3r)xr+(k=1(ak(k+r+1)(3k+3r1)+3ak1(k+r1))xk+r)=0a0cannot be 0 by assumption, giving the indicial equation(1+r)(1+3r)=0Values of r that satisfy the indicial equationr{1,13}Each term in the series must be 0, giving the recursion relation3(k+r+1)(k13+r)ak+3ak1(k+r1)=0Shift index usingk>k+13(k+2+r)(k+23+r)ak+1+3ak(k+r)=0Recursion relation that defines series solution to ODEak+1=3ak(k+r)(k+2+r)(3k+2+3r)Recursion relation forr=1; series terminates atk=1ak+1=3ak(k1)(k+1)(3k1)Apply recursion relation fork=0a1=3a0Terminating series solution of the ODE forr=1. Use reduction of order to find the second linearly independent solutiony(x)=a0(13x)Recursion relation forr=13ak+1=3ak(k+13)(k+73)(3k+3)Solution forr=13[y(x)=k=0akxk+13,ak+1=3ak(k+13)(k+73)(3k+3)]Combine solutions and rename parameters[y(x)=a0(13x)+(k=0bkxk+13),bk+1=3bk(k+13)(k+73)(3k+3)]
Mathematica. Time used: 0.763 (sec). Leaf size: 52
ode=x^2*D[y[x],{x,2}]+(5/3*x+x^2)*D[y[x],x]-1/3*y[x]==0; 
ic={}; 
DSolve[{ode,ic},y[x],x,IncludeSingularSolutions->True]
 
y(x)(3x1)(c21x9eK[1]K[1]3(13K[1])2dK[1]+c1)3x
Sympy. Time used: 1.515 (sec). Leaf size: 544
from sympy import * 
x = symbols("x") 
y = Function("y") 
ode = Eq(x**2*Derivative(y(x), (x, 2)) + (x**2 + 5*x/3)*Derivative(y(x), x) - y(x)/3,0) 
ics = {} 
dsolve(ode,func=y(x),ics=ics)
 
Solution too large to show