Difficult ODEs!!

[ Follow Ups ] [ Post Followup ] [ Netlib Discussion Forum ] [ FAQ ]

Posted by Lautaro Vergara on December 16, 1997 at 09:53:25:

Could somebody give some hints about how could I solve
the following stiff ODE?:

(x*y'(x)-3*y(x)=-1/sqrt(2+2*y'+4*x*y'') (1)

with BCs


where X is a suitably chosen 'asymptotic' value
(that could well be 1 or 2) and A is a parameter
that must be numerically found, such that the solution
does not diverge within the interval [0,X].

This is just an example for a more involved set of
coupled nonlinear ODEs of the form


where F is a simple function of their arguments,
like the LHS in (1), the integrals are all finite
(as in the RHS of (1)), G is more complicated but
well behaved function. The BCs that involve a set of
unknown parameters, like A above, that must be
found such that the solutions are finite. The BCs
are such that it there exists at most a discrete
set of acceptable solutions (this is known).

I have tried to implement relaxation, as in
Numerical Recipes, but I have found serious
I have also converted eq. (1) above as a DAE, but
I have found no code that could sove such system,
when the BCs are NOT consistent (COLDAE does this,
but only for consistent BCs).

I shall appreciate very much some hints about
this or information about a program that could
serve as a basis for solving such kind of problems.

Thank you very much in advance.


Lautaro Vergara

Follow Ups: