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?:
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.