29.11.12 problem 303

Internal problem ID [4903]
Book : Ordinary differential equations and their solutions. By George Moseley Murphy. 1960
Section : Various 11
Problem number : 303
Date solved : Sunday, March 30, 2025 at 04:10:07 AM
CAS classification : [_linear]

\begin{align*} \left (x^{2}+1\right ) y^{\prime }&=1+x^{2}-y \,\operatorname {arccot}\left (x \right ) \end{align*}

✓ Maple. Time used: 0.002 (sec). Leaf size: 23

ode:=(x^2+1)*diff(y(x),x) = 1+x^2-y(x)*arccot(x); 
dsolve(ode,y(x), singsol=all);
 

\[ y = \left (\int {\mathrm e}^{-\frac {\operatorname {arccot}\left (x \right )^{2}}{2}}d x +c_1 \right ) {\mathrm e}^{\frac {\operatorname {arccot}\left (x \right )^{2}}{2}} \]

✓ Mathematica. Time used: 3.489 (sec). Leaf size: 37

ode=(1+x^2)D[y[x],x]==(1+x^2)-y[x] ArcCot[x]; 
ic={}; 
DSolve[{ode,ic},y[x],x,IncludeSingularSolutions->True]
 

\[ y(x)\to e^{\frac {1}{2} \cot ^{-1}(x)^2} \left (\int _1^xe^{-\frac {1}{2} \cot ^{-1}(K[1])^2}dK[1]+c_1\right ) \]

✗ Sympy

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

Timed Out