MathGroup Archive 1999

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

Search the Archive

Re: Stochastic Differentail equations

  • To: mathgroup at
  • Subject: [mg19527] Re: Stochastic Differentail equations
  • From: Jens-Peer Kuska <kuska at>
  • Date: Sun, 29 Aug 1999 17:21:25 -0400
  • Organization: Universitaet Leipzig
  • References: <7q9ghb$>
  • Sender: owner-wri-mathgroup at

Hi Dimitris,

stochasitc differential equations need very special numerical methods.

The typical way in physics ist to make several assumptions about the
propertiers of the noise (spectral function ..) and obtain a partial
diffential equation. The Fokker-Planck equation is the most popular

An ordinary differntial equation solver make several assumptions about
derivatives of the rhight hand side. None of these assumptions are valid
for equations with a random variable inside.

I remember some Physical Review articles about the numerical simulation
of ode's with a stochasic rhs. Mail me if I should search for on or two

NDSolve[] uses high order methods to solve the differential equations
it needs that the derivatives of orders >3 exist.

Hope that helps

Dimitris Kugiumtzis wrote:
> I am not a regular reader of this newsgroup, so no flames if you
> find the question too simple or often repeated.
> I could not find in Mathematica how one can insert random components
> in the differential equations, when one uses NDSolve (giving rise
> to stochastic equations or as often referred to dynamic noise in the
> generated data). I am thinking something like for example
> x'[t]== -(y[t] + z[t])+Random
> where, the derivatives of y and z are defined accordingly (imagine
> for example the Lorenz equations in Mathematica with noise).
> If anyone can add something I would appreciate.
> Dimitris

  • Prev by Date: Re: Freeing memory in Mathematica
  • Next by Date: Re: Problem with conditional definitions
  • Previous by thread: Stochastic Differentail equations
  • Next by thread: GaussianQuadrature and truncated normal