MathGroup Archive 2001

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

Search the Archive

Re: right StoppingTest

  • To: mathgroup at smc.vnet.net
  • Subject: [mg31734] Re: right StoppingTest
  • From: bghiggins at ucdavis.edu (Brian Higgins)
  • Date: Thu, 29 Nov 2001 00:14:04 -0500 (EST)
  • References: <9tvgln$fql$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

Borut,
Here is a simple example that should help you get going:

eqns = {y'[t] == 4y[t], y[0] == 1}
sol = NDSolve[eqns, y[t], {t, 0, 10}, 
    StoppingTest :> If[y[t] > 200, lastT = t; True, False]]
Plot[y[t] /. sol, {t, 0, lastT}]

Cheers,

Brian




"Borut L" <borut at email.si> wrote in message news:<9tvgln$fql$1 at smc.vnet.net>...
> Hello,
> 
> I am solving a coupled system
> 
> o'[t] == a d[t] o[t] - b o[t]
> d'[t] == -a d[t] o[t]
> i'[t] == b o[t]
> 
> with NDSolve for o, d and i
> 
> There is a maximum for o[t] and I would like to stop the integration
> when it is reached and be given the value of o and time t when this happens.
> 
> I can't make the option StoppingTest to work it out properly.
> 
> Please suggest any help.
> 
> 
> Thanks,
> 
> Borut


  • Prev by Date: Compiling functions
  • Next by Date: Re: RE: A question about Plot
  • Previous by thread: Re: right StoppingTest
  • Next by thread: Re: right StoppingTest