April 21st, 2008, 5:36 pm
Technically, S=infty is a natural boundary of the underlying diffusion, so no boundary conditions are necessary or allowed.However, for numerical PDE truncations, it is nice to put in something reasonable.If you put in the wrong behavior V(t,Smax) = Smax^2, you will still get the correct answer with a PDE solverfor any finite S as Smax ->infty, IMO.The right behavior is really (*) V(t,S) ~ S^2 exp{(r + sig^2) t)}, as S >> 1To derive the right behavior for any payoff, use Paul's idea, or simply write the payoff as phi(y), using y = Log S(T). Then, the exact answer isV(t,S) = v(t,x), wherev(t,x) = exp(-r t) int exp{-(y - x - mu t)^2/(2 sig^2 t)} phi(y) dy/Sqrt[2 Pi sig^2 t] . Now x = log S(0) = log S. Let x >> 1, and the integral will reduce to something easy: you can simply take phi(y) ~ exp(2 y),which is how I got (*). Same thing works for any (integrable) payoff. You couldanalyze payoff's that grow as S^p for some power once and for all by this method and be done with it.regards,
Last edited by
Alan on April 20th, 2008, 10:00 pm, edited 1 time in total.