29.8.1 problem 206

Internal problem ID [4806]
Book : Ordinary differential equations and their solutions. By George Moseley Murphy. 1960
Section : Various 8
Problem number : 206
Date solved : Sunday, March 30, 2025 at 03:58:58 AM
CAS classification : [[_homogeneous, `class A`], _dAlembert]

\begin{align*} x y^{\prime }&=y+x \sec \left (\frac {y}{x}\right )^{2} \end{align*}

✓ Maple. Time used: 0.027 (sec). Leaf size: 31

ode:=x*diff(y(x),x) = y(x)+x*sec(y(x)/x)^2; 
dsolve(ode,y(x), singsol=all);
 

\[ \frac {x \sin \left (\frac {2 y}{x}\right )+2 y}{4 x}-\ln \left (x \right )-c_1 = 0 \]

✓ Mathematica. Time used: 0.283 (sec). Leaf size: 31

ode=x D[y[x],x]==y[x]+x Sec[y[x]/x]^2; 
ic={}; 
DSolve[{ode,ic},y[x],x,IncludeSingularSolutions->True]
 

\[ \text {Solve}\left [\frac {y(x)}{2 x}+\frac {1}{4} \sin \left (\frac {2 y(x)}{x}\right )=\log (x)+c_1,y(x)\right ] \]

✗ Sympy

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

TypeError : cannot determine truth value of Relational