✓ Maple : cpu = 0.037 (sec), leaf count = 28
This is Bernoulli first order non-linear ODE.
Let
Integrating factor is
Hence
Let
Verification
restart; eq:=diff(y(x),x)+f(x)*y(x)^2+g(x)*y(x) = 0; beta:=exp(-Int(g(x),x)): my_sol:=beta/(Int(f(x)*beta,x)+_C1); odetest(y(x)=my_sol,eq); 0