3.72 problem 1073

Internal problem ID [9405]
Internal file name [OUTPUT/8342_Monday_June_06_2022_02_47_10_AM_86850837/index.tex]

Book: Differential Gleichungen, E. Kamke, 3rd ed. Chelsea Pub. NY, 1948
Section: Chapter 2, linear second order
Problem number: 1073.
ODE order: 2.
ODE degree: 1.

The type(s) of ODE detected by this program : "unknown"

Maple gives the following as the ode type

[[_2nd_order, _with_linear_symmetries]]

Unable to solve or complete the solution.

\[ \boxed {y^{\prime \prime }+\frac {\left (11 \operatorname {WeierstrassP}\left (x , a , b\right ) \operatorname {WeierstrassPPrime}\left (x , a , b\right )-6 \operatorname {WeierstrassP}\left (x , a , b\right )^{2}+\frac {a}{2}\right ) y^{\prime }}{\operatorname {WeierstrassPPrime}\left (x , a , b\right )+\operatorname {WeierstrassP}\left (x , a , b\right )^{2}}+\frac {\left (\operatorname {WeierstrassPPrime}\left (x , a , b\right )^{2}-\operatorname {WeierstrassP}\left (x , a , b\right )^{2} \operatorname {WeierstrassPPrime}\left (x , a , b\right )-\operatorname {WeierstrassP}\left (x , a , b\right ) \left (6 \operatorname {WeierstrassP}\left (x , a , b\right )^{2}-\frac {a}{2}\right )\right ) y}{\operatorname {WeierstrassPPrime}\left (x , a , b\right )+\operatorname {WeierstrassP}\left (x , a , b\right )^{2}}=0} \]

Maple trace

`Methods for second order ODEs: 
--- Trying classification methods --- 
trying a symmetry of the form [xi=0, eta=F(x)] 
checking if the LODE is missing y 
-> Heun: Equivalence to the GHE or one of its 4 confluent cases under a power @ Moebius 
-> trying a solution of the form r0(x) * Y + r1(x) * Y where Y = exp(int(r(x), dx)) * 2F1([a1, a2], [b1], f) 
-> Trying changes of variables to rationalize or make the ODE simpler 
   trying a symmetry of the form [xi=0, eta=F(x)] 
   checking if the LODE is missing y 
   -> Heun: Equivalence to the GHE or one of its 4 confluent cases under a power @ Moebius 
   -> trying a solution of the form r0(x) * Y + r1(x) * Y where Y = exp(int(r(x), dx)) * 2F1([a1, a2], [b1], f) 
      trying a symmetry of the form [xi=0, eta=F(x)] 
      trying 2nd order exact linear 
      trying symmetries linear in x and y(x) 
      trying to convert to a linear ODE with constant coefficients 
<- unable to find a useful change of variables 
   trying a symmetry of the form [xi=0, eta=F(x)] 
   trying 2nd order exact linear 
   trying symmetries linear in x and y(x) 
   trying to convert to a linear ODE with constant coefficients 
   trying 2nd order, integrating factor of the form mu(x,y) 
   trying a symmetry of the form [xi=0, eta=F(x)] 
   checking if the LODE is missing y 
   -> Heun: Equivalence to the GHE or one of its 4 confluent cases under a power @ Moebius 
   -> trying a solution of the form r0(x) * Y + r1(x) * Y where Y = exp(int(r(x), dx)) * 2F1([a1, a2], [b1], f) 
   -> Trying changes of variables to rationalize or make the ODE simpler 
      trying a symmetry of the form [xi=0, eta=F(x)] 
      checking if the LODE is missing y 
      -> Heun: Equivalence to the GHE or one of its 4 confluent cases under a power @ Moebius 
      -> trying a solution of the form r0(x) * Y + r1(x) * Y where Y = exp(int(r(x), dx)) * 2F1([a1, a2], [b1], f) 
         trying a symmetry of the form [xi=0, eta=F(x)] 
         trying 2nd order exact linear 
         trying symmetries linear in x and y(x) 
         trying to convert to a linear ODE with constant coefficients 
   <- unable to find a useful change of variables 
      trying a symmetry of the form [xi=0, eta=F(x)] 
   trying to convert to an ODE of Bessel type 
   -> trying reduction of order to Riccati 
      trying Riccati sub-methods: 
         trying Riccati_symmetries 
         -> trying a symmetry pattern of the form [F(x)*G(y), 0] 
         -> trying a symmetry pattern of the form [0, F(x)*G(y)] 
         -> trying a symmetry pattern of the form [F(x),G(x)*y+H(x)] 
   -> trying with_periodic_functions in the coefficients 
      --- Trying Lie symmetry methods, 2nd order --- 
      `, `-> Computing symmetries using: way = 5`[0, y]
 

Solution by Maple

dsolve(diff(diff(y(x),x),x)+(11*WeierstrassP(x,a,b)*WeierstrassPPrime(x,a,b)-6*WeierstrassP(x,a,b)^2+1/2*a)*diff(y(x),x)/(WeierstrassPPrime(x,a,b)+WeierstrassP(x,a,b)^2)+(WeierstrassPPrime(x,a,b)^2-WeierstrassP(x,a,b)^2*WeierstrassPPrime(x,a,b)-WeierstrassP(x,a,b)*(6*WeierstrassP(x,a,b)^2-1/2*a))*y(x)/(WeierstrassPPrime(x,a,b)+WeierstrassP(x,a,b)^2)=0,y(x), singsol=all)
 

\[ \text {No solution found} \]

Solution by Mathematica

Time used: 0.0 (sec). Leaf size: 0

DSolve[((-(WeierstrassP[x, {a, b}]*(-1/2*a + 6*WeierstrassP[x, {a, b}]^2)) - WeierstrassP[x, {a, b}]^2*WeierstrassPPrime[x, {a, b}] + WeierstrassPPrime[x, {a, b}]^2)*y[x])/(WeierstrassP[x, {a, b}]^2 + WeierstrassPPrime[x, {a, b}]) + ((a/2 - 6*WeierstrassP[x, {a, b}]^2 + WeierstrassP[x, {a, b}]^3 - WeierstrassP[x, {a, b}]*WeierstrassPPrime[x, {a, b}])*y'[x])/(-WeierstrassP[x, {a, b}]^2 + WeierstrassPPrime[x, {a, b}]) + y''[x] == 0,y[x],x,IncludeSingularSolutions -> True]
 

Not solved