4.1.41 \(f(x)^2+y'(x)=f'(x)+y(x)^2\)

ODE
\[ f(x)^2+y'(x)=f'(x)+y(x)^2 \] ODE Classification

odeadvisor timed out

Book solution method
Riccati ODE, Generalized ODE

Mathematica
cpu = 15.1413 (sec), leaf count = 0 , could not solve

DSolve[f[x]^2 + Derivative[1][y][x] == y[x]^2 + Derivative[1][f][x], y[x], x]

Maple
cpu = 0.308 (sec), leaf count = 30

\[\left [y \left (x \right ) = f \left (x \right )+\frac {{\mathrm e}^{\int 2 f \left (x \right )d x}}{\textit {\_C1} -\left (\int {\mathrm e}^{\int 2 f \left (x \right )d x}d x \right )}\right ]\] Mathematica raw input

DSolve[f[x]^2 + y'[x] == y[x]^2 + f'[x],y[x],x]

Mathematica raw output

DSolve[f[x]^2 + Derivative[1][y][x] == y[x]^2 + Derivative[1][f][x], y[x], x]

Maple raw input

dsolve(diff(y(x),x)+f(x)^2 = diff(f(x),x)+y(x)^2, y(x))

Maple raw output

[y(x) = f(x)+exp(Int(2*f(x),x))/(_C1-Int(exp(Int(2*f(x),x)),x))]