Quantnotes.com
Help Contact
Fundamentals
Publications
Software & Data
Book Reviews
Job Listings
Event Listings
Forums
Edutainment
Useful Links
About Us

 

QuantFinanceJobs.com

:: Collaborations ::

:: QF Journal ::

 

:: Solving the Black-Scholes Equation ::

We have the Black-Scholes equation for European non-dividend-paying call option

[Graphics:Images/solvingbs_gr_1.gif]

with boundary conditions

[Graphics:Images/solvingbs_gr_2.gif]

and

[Graphics:Images/solvingbs_gr_3.gif]

These conditions  for a European call option with expiry T and exercise price E arise from: an option being worthless if the price of the underlying security is zero itself; the price of an option C(t) cannot be greater than the price of the underlying security S(t) otherwise an aribtrage can achieved by buying the security and selling the call option.

In fact, C(t) <= S(t), i.e. the maximum value of the call option is the difference between the spot price and the exercise price.

The Black-Scholes equation resembles a linear parabolic equation similar to the diffusion equation.  A general linear parabolic equation is of the form:

[Graphics:Images/solvingbs_gr_4.gif]

where a and b are constants.  This can always be reduced to a diffusion equation

[Graphics:Images/solvingbs_gr_5.gif]

by choosing the  substitution

[Graphics:Images/solvingbs_gr_6.gif]

We find that this condition can be satisfied by setting

[Graphics:Images/solvingbs_gr_7.gif]

and

[Graphics:Images/solvingbs_gr_8.gif]

If we assume r and [Graphics:Images/solvingbs_gr_9.gif] to be constant, the Black-Scholes equation can be reduced to the diffusion equation.  To make this possible we set

[Graphics:Images/solvingbs_gr_10.gif]

The result is

[Graphics:Images/solvingbs_gr_11.gif]

where

[Graphics:Images/solvingbs_gr_12.gif]

Having used this substitution, our initial condition becomes

[Graphics:Images/solvingbs_gr_13.gif]

Now our equation looks more like the diffusion equation and we can obtain the similarity solution using

[Graphics:Images/solvingbs_gr_14.gif]

By substituting this expression for v into the parabolic equation above and choosing values of [Graphics:Images/solvingbs_gr_15.gif] and [Graphics:Images/solvingbs_gr_16.gif] such that all terms involving u(x, [Graphics:Images/solvingbs_gr_17.gif]) are eliminated, we obtain

[Graphics:Images/solvingbs_gr_18.gif]

with

[Graphics:Images/solvingbs_gr_19.gif]

For now it is just a matter of solving the diffusion equation by methods of Fourier Transform.

The final solution is therefore

[Graphics:Images/solvingbs_gr_20.gif]

where

[Graphics:Images/solvingbs_gr_21.gif]

and

[Graphics:Images/solvingbs_gr_22.gif]

Here N(.) is the Cumulative Normal Distribution Function, the exact form is given by:

[Graphics:Images/solvingbs_gr_23.gif]

A simpler polynomial approximation up to 6 decimal places accurate can be used to represent this function instead:

[Graphics:Images/solvingbs_gr_24.gif]

and

[Graphics:Images/solvingbs_gr_25.gif]

Here

[Graphics:Images/solvingbs_gr_26.gif]

This is the exact solution for calculating the value of an European non-dividend paying call option. For put option of the same expiry T and exercise price E, the put-call parity formula can be used if the value of the call option is first calculated

[Graphics:Images/solvingbs_gr_27.gif]

This gives value of P as

[Graphics:Images/solvingbs_gr_28.gif]

 

Written by Alessio Farhadi.

Back for more

 

Copyright © 2001-04 Quantnotes.com. All rights reserved. Legal Notice | Privacy Notice