MathGroup Archive 2006

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

Search the Archive

Re: Discrete Event handling in Mathematica

  • To: mathgroup at smc.vnet.net
  • Subject: [mg68034] Re: Discrete Event handling in Mathematica
  • From: "ben" <benjamin.friedrich at gmail.com>
  • Date: Fri, 21 Jul 2006 05:37:21 -0400 (EDT)
  • References: <e9nl3f$a0k$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

Hi Raj,

Does this do the job?

NDSolve[{\[Theta]'[t] == Sin[\[Theta][t]], \[Theta][0] == 0.5},
\[Theta], {t,
    0, 100}, Method -> {EventLocator, "Event" -> \[Theta]l (t - tl) -
0.01,
      "EventAction" :> Print["tl=", tl, ", \[Theta]=", \[Theta][t]]}]

Bye
Ben


Raj schrieb:

> Hi!
>
> Could somebody tell me if Mathematica can handle discrete events in
> Ordinary Differential Equations.
>
> For example in the following code
>
> Clear[\[Theta], \[Theta]l, tl]
> \[Theta]l = Sin[0.5]
> tl = 0;
> Reap[NDSolve[{\[Theta]'[t] == Sin[\[Theta][t]], \[Theta][0] ==
>         0.5}, \[Theta], {t, 0, 100},
>     Method -> {EventLocator, "Event" -> \[Theta]l (t - tl) - 0.01,
>         "EventAction" :> Sow[{tl = t, \[Theta]l = Sin[\[Theta][t]]}]}]]
>
>  I want Mathematica to print out tl and \[Theta]l everytime the "Event"
> condition is met.
> 
> Thanks,
> 
> Raj


  • Prev by Date: Re: MapThread and If
  • Next by Date: Re: MapThread and If
  • Previous by thread: Re: Discrete Event handling in Mathematica
  • Next by thread: Re: Discrete Event handling in Mathematica