Код
dy/dx=a(y/x^2)-2
При этом ход решения:
Код
dy/dx=ay/x^2 => lny=-(a/x)+ln(C(x))
y=C(x)exp(-a/x)
=>
!! dc/dx=-2exp(a/x)
y=C(x)exp(-a/x)
=>
!! dc/dx=-2exp(a/x)
Решение для данного случая упирается в этот частный интеграл.
Код
c=-2(x*exp(a/x)-aEi(a/x))
a>0, x и y - больше либо раны нулю. При этом сначала x=max и y=0, с другой стороны предельное значение x=min, y=0