nijso beishuizen wrote:
> Hi,
>
> The equation I try to solve is the final step in the solution of a
> nonlinear ode. The solution in the paper
> (of Cheb-Terrab and Kolokolnikov, it's on ArXiv) is given in implicit
> form as:
>
> sol : x + 1/(3*x^3*(y+x)^3)+x^(1+a)/(1+a) = %c;
>
Nijso,
Is this from E. S. Cheb-Terrab, T. Koloknikov, First Order ODEs,
Symmetries and Linear Transformations, European Journal of Applied
Mathematics, Vol. 14, No. 2, pp. 231-246 (2003)? I wrote an
implementation of this as part of the contrib_ode function. Look in
the share/contrib/diffequations/ directory. The code is ode1_lie.mac
and the examples and unit tests are in the tests subdirectory
The code works well for some cases but runs "forever" for other cases.
There are a number of test cases from Kamke and Murphy that are
commented out for this reason. I understand that the Maple routine
doesn't have the same problem. I haven't been able to fix this.
David