2.1.21 Problem 21

Internal problem ID [9005]
Book : First order enumerated odes
Section : section 1
Problem number : 21
Date solved : Sunday, March 30, 2025 at 01:58:38 PM
CAS classification : [_rational, _Riccati]

Unknown ode type.

✓ Maple. Time used: 0.008 (sec). Leaf size: 94

ode:=c*diff(y(x),x) = (a*x+b*y(x)^2)/r/x; 
dsolve(ode,y(x), singsol=all);
 

$$y=\frac{\sqrt{\frac{xab}{r^2{c}^2}}cr(\mathrm{BesselY}(1,2\sqrt{\frac{xab}{r^2{c}^2}})c_1+\mathrm{BesselJ}(1,2\sqrt{\frac{xab}{r^2{c}^2}}))}{b(c_1\mathrm{BesselY}(0,2\sqrt{\frac{xab}{r^2{c}^2}})+\mathrm{BesselJ}(0,2\sqrt{\frac{xab}{r^2{c}^2}}))}$$

Maple trace

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 sub-methods: 
      -> Trying a Liouvillian solution using Kovacics algorithm 
      <- No Liouvillian solutions exists 
   <- Abel AIR successful: ODE belongs to the 0F1 1-parameter (Bessel type) cla\ 
ss
 

Maple step by step

$$\begin{array}{lll}& & \text{Let’s solve}\\ & & c\left(\frac{d}{dx}y\left(x\right)\right)=\frac{ax+by{\left(x\right)}^2}{rx}\\ \bullet & & \text{Highest derivative means the order of the ODE is}1\\ & & \frac{d}{dx}y\left(x\right)\\ \bullet & & \text{Solve for the highest derivative}\\ & & \frac{d}{dx}y\left(x\right)=\frac{ax+by{\left(x\right)}^2}{rxc}\end{array}$$

✓ Mathematica. Time used: 0.291 (sec). Leaf size: 207

ode=c*D[y[x],x]==(a*x+b*y[x]^2)/(r*x); 
ic={}; 
DSolve[{ode,ic},y[x],x,IncludeSingularSolutions->True]
 

$$\begin{array}{r}y(x)\to \frac{\sqrt{a}\sqrt{x}(2\mathrm{BesselY}(1,\frac{2\sqrt{a}\sqrt{b}\sqrt{x}}{cr})+c_1\mathrm{BesselJ}(1,\frac{2\sqrt{a}\sqrt{b}\sqrt{x}}{cr}))}{\sqrt{b}(2\mathrm{BesselY}(0,\frac{2\sqrt{a}\sqrt{b}\sqrt{x}}{cr})+c_1\mathrm{BesselJ}(0,\frac{2\sqrt{a}\sqrt{b}\sqrt{x}}{cr}))}\\ y(x)\to \frac{\sqrt{a}\sqrt{x}\mathrm{BesselJ}(1,\frac{2\sqrt{a}\sqrt{b}\sqrt{x}}{cr})}{\sqrt{b}\mathrm{BesselJ}(0,\frac{2\sqrt{a}\sqrt{b}\sqrt{x}}{cr})}\end{array}$$

✗ Sympy

from sympy import * 
x = symbols("x") 
a = symbols("a") 
b = symbols("b") 
c = symbols("c") 
r = symbols("r") 
y = Function("y") 
ode = Eq(c*Derivative(y(x), x) - (a*x + b*y(x)**2)/(r*x),0) 
ics = {} 
dsolve(ode,func=y(x),ics=ics)
 

RecursionError : maximum recursion depth exceeded