MathGroup Archive 2006

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Relatively simple, but problematic, non-linear ODE


<ab_def at prontomail.com> wrote in message news:elbi31$k16$1 at smc.vnet.net...

> f[rho_?NonNegative] := Module[{y, z},
>  y = y /. First@ NDSolve[
>    {y'[z] == 1/2*(rho*y[z] + Sqrt[(rho^2 - 1)*y[z]^2 + 4]),
>     y[0] == 0},
>    y, {z, 0, Pi},
>    StoppingTest -> Re[y[z]] - 2];
>  Function @@ {z,
>    Piecewise[{{y[z], z < y[Domain[]][[1, 2]]}}, 2]}
> ]

 Thanks a lot, Maxim -- I see now, using your suggestion,
 that the StoppingTest option is the key here.

regards,
alan 


  • Prev by Date: Re: Problem with text rendering on Linux.
  • Next by Date: Re: On Reduce
  • Previous by thread: Re: Relatively simple, but problematic, non-linear ODE
  • Next by thread: Re: Points sampled by N[Derivative[]]