MathGroup Archive 1999

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

Search the Archive

help with pde

  • To: mathgroup at
  • Subject: [mg16110] help with pde
  • From: "Will Holt" <will at>
  • Date: Thu, 25 Feb 1999 08:25:03 -0500
  • Organization: Customer of Planet Online
  • Sender: owner-wri-mathgroup at

Hi everyone,

I am trying to solve the Black-Scholes pde in mathematica subject to the
usual boundary conditions.


V[0,T]=0,V[S,T]=Max[S-K,0], for some particular S[0], K and T.

The Mathematica command "NDSolve" requires that the first argument must have
both an equation and an initial condition.
The problem is that the law motion for S is


where dZ is a Wiener process that can be substituted by epsilon*SQRt[dt],
with epsilon a random drawing from a standardised normal distribution.

I thought that I did not have to include this law motion for S in NDSolve,
but in case I do have to include it, how can I do this?

Also, on a different issue, how can I generate three series of normal random
numbers that are correlated amongst each other: rho12, rho23, rho13? e.g.
dX1=rho12dX2+SQRT[1-rho12^2]*de12, with de12 a standard Wiener process
independent of dX2, and so on for dX2 and dX3.

Any help will be greatly appreciated.


  • Prev by Date: Re: Evaluate/Module
  • Next by Date: model boat hull design methods
  • Previous by thread: Re: MathLink & Metrowerks CodeWarrior
  • Next by thread: Re: help with pde