|
:: Solving the Black-Scholes Equation ::
We have the Black-Scholes equation for European non-dividend-paying
call option
![[Graphics:Images/solvingbs_gr_1.gif]](Images/solvingbs_gr_1.gif)
with boundary conditions
![[Graphics:Images/solvingbs_gr_2.gif]](Images/solvingbs_gr_2.gif)
and
![[Graphics:Images/solvingbs_gr_3.gif]](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]](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]](Images/solvingbs_gr_5.gif)
by choosing the substitution
![[Graphics:Images/solvingbs_gr_6.gif]](Images/solvingbs_gr_6.gif)
We find that this condition can be satisfied by setting
![[Graphics:Images/solvingbs_gr_7.gif]](Images/solvingbs_gr_7.gif)
and
![[Graphics:Images/solvingbs_gr_8.gif]](Images/solvingbs_gr_8.gif)
If we assume r and
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]](Images/solvingbs_gr_10.gif)
The result is
![[Graphics:Images/solvingbs_gr_11.gif]](Images/solvingbs_gr_11.gif)
where
![[Graphics:Images/solvingbs_gr_12.gif]](Images/solvingbs_gr_12.gif)
Having used this substitution, our initial condition becomes
![[Graphics:Images/solvingbs_gr_13.gif]](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]](Images/solvingbs_gr_14.gif)
By substituting this expression for v into the parabolic equation above
and choosing values of
and
such that all terms involving u(x, )
are eliminated, we obtain
![[Graphics:Images/solvingbs_gr_18.gif]](Images/solvingbs_gr_18.gif)
with
![[Graphics:Images/solvingbs_gr_19.gif]](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]](Images/solvingbs_gr_20.gif)
where
![[Graphics:Images/solvingbs_gr_21.gif]](Images/solvingbs_gr_21.gif)
and
![[Graphics:Images/solvingbs_gr_22.gif]](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]](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]](Images/solvingbs_gr_24.gif)
and
![[Graphics:Images/solvingbs_gr_25.gif]](Images/solvingbs_gr_25.gif)
Here
![[Graphics:Images/solvingbs_gr_26.gif]](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]](Images/solvingbs_gr_27.gif)
This gives value of P as
![[Graphics:Images/solvingbs_gr_28.gif]](Images/solvingbs_gr_28.gif)
Written by Alessio Farhadi.
Back for more 
|