2.1.91 Problem 89

Maple
Mathematica
Sympy

Internal problem ID [8802]
Book : Own collection of miscellaneous problems
Section : section 1.0
Problem number : 89
Date solved : Friday, April 25, 2025 at 05:11:34 PM
CAS classification : [[_2nd_order, _exact, _nonlinear], [_2nd_order, _reducible, _mu_x_y1], [_2nd_order, _reducible, _mu_y_y1], [_2nd_order, _reducible, _mu_xy]]

Maple. Time used: 0.053 (sec). Leaf size: 147
ode:=diff(diff(y(x),x),x)-y(x)*diff(y(x),x) = 2*x; 
dsolve(ode,y(x), singsol=all);
 
y=WhittakerM(ic128+1,14,i2x22)(ic12+6)+8c2WhittakerW(ic128+1,14,i2x22)+2(1i(x2c12)2)(c2WhittakerW(ic128,14,i2x22)+WhittakerM(ic128,14,i2x22))2x(c2WhittakerW(ic128,14,i2x22)+WhittakerM(ic128,14,i2x22))

Maple trace

Methods for second order ODEs: 
--- Trying classification methods --- 
trying 2nd order Liouville 
trying 2nd order WeierstrassP 
trying 2nd order JacobiSN 
differential order: 2; trying a linearization to 3rd order 
trying 2nd order ODE linearizable_by_differentiation 
trying 2nd order, 2 integrating factors of the form mu(x,y) 
trying differential order: 2; missing variables 
-> trying 2nd order, dynamical_symmetries, fully reducible to Abel through one \ 
integrating factor of the form G(x,y)/(1+H(x,y)*y)^2 
   --- trying a change of variables {x -> y(x), y(x) -> x} and re-entering meth\ 
ods for dynamical symmetries --- 
   -> trying 2nd order, dynamical_symmetries, fully reducible to Abel through o\ 
ne integrating factor of the form G(x,y)/(1+H(x,y)*y)^2 
trying 2nd order, integrating factors of the form mu(x,y)/(y)^n, only the sing\ 
ular cases 
trying symmetries linear in x and y(x) 
trying differential order: 2; exact nonlinear 
-> Calling odsolve with the ODE, diff(_b(_a),_a) = 1/2*_b(_a)^2+_a^2-_C1, _b(_a 
) 
   *** Sublevel 2 *** 
   Methods for first order ODEs: 
   --- Trying classification methods --- 
   trying a quadrature 
   trying 1st order linear 
   trying Bernoulli 
   trying separable 
   trying inverse linear 
   trying homogeneous types: 
   trying Chini 
   differential order: 1; looking for linear symmetries 
   trying exact 
   Looking for potential symmetries 
   trying Riccati 
   trying Riccati Special 
   trying Riccati sub-methods: 
      trying Riccati to 2nd Order 
      -> Calling odsolve with the ODE, diff(diff(y(x),x),x) = (-1/2*x^2+1/2*_C1 
)*y(x), y(x) 
         *** Sublevel 3 *** 
         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 
         <- No Liouvillian solutions exists 
         -> Trying a solution in terms of special functions: 
            -> Bessel 
            -> elliptic 
            -> Legendre 
            -> Whittaker 
               -> hyper3: Equivalence to 1F1 under a power @ Moebius 
               <- hyper3 successful: received ODE is equivalent to the 1F1 ODE 
            <- Whittaker successful 
         <- special function solution successful 
      <- Riccati to 2nd Order successful 
<- differential order: 2; exact nonlinear successful
 

Mathematica. Time used: 48.14 (sec). Leaf size: 318
ode=D[y[x],{x,2}]+D[y[x],x]*y[x]==2*x; 
ic={}; 
DSolve[{ode,ic},y[x],x,IncludeSingularSolutions->True]
 
y(x)24(24xParabolicCylinderD(14(2c12),i24x)+2iParabolicCylinderD(14(22c1),i24x)+c2(2ParabolicCylinderD(14(2c1+2),24x)24xParabolicCylinderD(14(2c12),24x)))ParabolicCylinderD(14(2c12),i24x)+c2ParabolicCylinderD(14(2c12),24x)y(x)2x224ParabolicCylinderD(14(2c1+2),24x)ParabolicCylinderD(14(2c12),24x)y(x)2x224ParabolicCylinderD(14(2c1+2),24x)ParabolicCylinderD(14(2c12),24x)
Sympy
from sympy import * 
x = symbols("x") 
y = Function("y") 
ode = Eq(-2*x - y(x)*Derivative(y(x), x) + Derivative(y(x), (x, 2)),0) 
ics = {} 
dsolve(ode,func=y(x),ics=ics)
 
NotImplementedError : The given ODE -(-2*x + Derivative(y(x), (x, 2)))/y(x) + Derivative(y(x), x) cannot be solved by the factorable group method