link to University of Nebraska Lincoln
Department of Mathematics
203 Avery Hall
University of Nebraska Lincoln
Lincoln, NE 68588-0323
402-472-3731 (voice)
402-472-8466 (fax)

Math 489/Math 889
Stochastic Processes and
Advanced Mathematical Finance
Dunbar, Fall 2007


Solution of the Black Scholes Equation


Key Concepts

Key Concepts

  1. The Black-Scholes equation for the value of a European call option on a security can be solved by judicious changes of variables that reduce the equation to the heat equation, for which we can write down the solution.
  2. An example of how to choose and execute changes of variables to solve a partial differential equation.


Vocabulary

Vocabulary

  1. A differential equation with auxiliary initial conditions and boundary conditions, that is an initial value problem, is said to be well-posed if the solution exists, is unique, and small changes in the equation parameters, the initial conditions or the boundary conditions produce only small changes in the solutions.


Mathematical Ideas

Mathematical Ideas

Analytic Solution of the Black-Scholes Equation

This discussion is drawn from Section 4.2, pages 59-63; Section 4.3, pages 66-69; Section 5.3, pages 75-76; and Section 5.4, pages 77 - 81 of The Mathematics of Financial Derivatives: A Student Introduction by P. Wilmott, S. Howison, J. Dewynne, Cambridge University Press, Cambridge, 1995. Some ideas are also taken from Chapter 11 of Stochastic Calculus and Financial Applications by J. Michael Steele, Springer, New York, 2001.

We have to start somewhere, and to avoid the problem of deriving everything back to calculus, we will assert that the initial value problem for the heat equation on the infinite interval is well-posed, that is, the solution to the partial differential equation

 2 @u- = @--u - oo < x < oo , t > 0 @t @x2

with

u(x,0) = u0(x)

where the initial condition and the solution satisfy the following technical requirements:

  1. u0(x) has no more than a finite number of discontinuities of the jump kind,
  2. lim |x|--> oo u0(x)e-ax2 = 0 for any a > 0,
  3. lim |x|--> oo u(x,t)e-ax2 = 0 for any a > 0

exists for all time, is unique, and most importantly, can be represented as

 1 integral oo 2 u(x,t ) = - V~ --- u0(s)e- (x- s) /4t ds 2 pt - oo

Remark: This solution can derived in several different ways, I think the easiest way is to use Fourier transforms. The derivation of this solution representation is standard in any course or book on partial differential equations.

Remark: Mathematically, the conditions above are unnecessarily restrictive, and can be considerably weakened, however, they will be more than sufficient for all practical situations we encounter in mathematical finance.

Remark: The use of t for the time variable (instead of the more natural t) is to avoid a conflict of notation in the several changes of variables we will soon have to make.

The Black-Scholes terminal value problem for the value V (S,t) of a European call option on a security with price S at time t is

 2 @V--+ 1s2S2 @-V--+ rS @V--- rV = 0 @t 2 @S2 @S

with V (0,t) = 0, V (S,t) ~ S as S --> oo and

V (S,T ) = max(S - K, 0).

Note that this looks a little like the heat equation on the infinite interval in that it has a first derivative of the unknown with respect to time and the second derivative of the unknown with respect to the other (space) variable, but notice:

  1. Each time the unknown is differentiated with respect to S, it also multiplied by the independent variable S, so the equation is not a constant coefficient equation.
  2. There is a first derivative of V with respect to S in the equation.
  3. There is a zero-th order term V in the equation.
  4. The sign on the second derivative is the opposite of the heat equation form, so the equation is of backward parabolic form.
  5. The data of the problem is given at the final time T instead of the initial time 0, consistent with the backward parabolic form of the equation.
  6. There is a boundary condition V (0,t) = 0 specifying the value of the solution at one sensible boundary of the problem. The boundary is sensible since security values must only be zero or positive. This boundary condition says that any time the security value is 0, then the call value (with strike price K) is also worth 0.
  7. There is another boundary condition V (S,t) ~ S, as S -->  oo , but although this is financially sensible, (it says that for very large security prices, the call value with strike price K is approximately S) it is more in the nature of a technical condition, and we will ignore it without consequence.

We will eliminate each and every one of these objections with a suitable change of variables. The plan is to change variables to reduce the Black-Scholes terminal value problem to the heat equation, and then to use the known solution of the heat equation to represent the solution, and change variables back. This is a standard technique of solution in partial differential equations, and none of the transformations we are making are strange, unmotivated, or unknown.

First we will take t = T - t (1/2)s2 and S = Kex, and we will set

V(S, t) = Kv(x, t).

Remember, s is the volatility, r is the interest rate on a risk-free bond, and K is the strike price. In the changes of variables above, the choice for t reverses the sense of time, changing the problem from backward parabolic to forward parabolic. The choice for S is a well-known transformation based on experience with the Euler equidimensional equation in differential equations. In addition, the variables have been carefully scaled so as to make the transformed equation expressed in dimensionless quantities. All of these techniques are standard and are covered in most courses and books on partial differential equations and applied mathematics.

Some extremely wise advice adapted from page 186 of Stochastic Calculus and Financial Applications by J. Michael Steele. Springer, New York, 2001 is appropriate here.

“There is nothing particularly difficult about changing variables and transforming one equation to another, but there is an element of tedium and complexity that slows us down. There is no universal remedy for this molasses effect, but the calculations do seem to go more quickly if one follows a well-defined plan. If we know that V (S,t) satisfies an equation (like the Black-Scholes equation) we are guaranteed that we can make good use of the equation in the derivation of the equation for a new function v(x,t) defined in terms of the old if we write the old V as a function of the new v and write the new t and x as functions of the old t and S. This order of things puts everything in the direct line of fire of the chain rule; the partial derivatives V t, V S and V SS are easy to compute and at the end, the original equation stands ready for immediate use.”

Following our advice, we first have to write

 2 t = (1/2)s (T - t)

and

x = log(S/K)

so

@V @v dt @v - s2 ----= K ---.--- = K --- .----- @t @t dt @t 2

and

@V--= K @v-.dx- = K @v-.1- @S @x dS @x S

and

 2 ( ) @-V-- = @-- @V-- @S2 @S @S @ ( @v 1 ) = --- K ----- @S @x S ( ) @v- --1 @-- @v- -1 = K @x . S2 + K @S @x .S ( ) = K @v-.--1 + K @-- @v- .dx-.1-- @x S2 @x @x dS S @v - 1 @2v 1 = K ---.--- + K ----.--- @x S2 @x2 S2

and the terminal condition

V(S, T) = max(S - K, 0) = max(Kex - K, 0)

but V (S,T) = Kv(x, 0) so v(x, 0) = max(ex - 1, 0).

Now substitute all of the derivatives into the Black-Scholes equation to obtain:

 2 2 ( 2 ) ( ) K @v-. -s--+ s-S2 K @v-. -1-+ K @-v-.1-- + rS K @v-.1-- - rKv = 0 @t 2 2 @x S2 @x2 S2 @x S

Now the simplification begins:

  1. Factor the common factor K and cancel.
  2. Transpose the t-derivative to the other side, and divide through by s2/2
  3. Rename the remaining constant r/(s2/2) as k. k measures the ratio between the risk-free interest rate and the volatility.
  4. Cancel the S2 terms in the second derivative.
  5. Cancel the S terms in the first derivative.
  6. Gather up like order terms.

What remains is the rescaled, constant coefficient equation:

 2 @v-= @-v-+ (k - 1) @v-- kv. @t @x2 @x

We have made considerable progress, because

  1. Now there is only one dimensionless parameter k measuring the risk-free interest rate as a multiple of the volatility and another (more or less hidden) rescaled time to expiry (1/2)s2T , not the original 4 dimensioned quantities, namely K, T, s2 and r.
  2. The equation is defined on the interval - oo < x <  oo , since this x-interval defines 0 < S <  oo through the change of variables S = Kex.
  3. The equation is now constant coefficient. In principle, we could solve the equation directly now,

Instead, we will change the dependent variable scale yet again, by

 ax+bt v = e u(x,t).

where a and b are yet to be determined. Now using the product rule:

vt = beax+btu(x, t) + eax+btut

and

 ax+bt ax+bt vx = ae u(x,t) + e ux

and

v = a2eax+bt u(x,t ) + 2aeax+bt u + eax+btu xx x xx

Put these into our constant coefficient partial differential equation, cancel the common factor of eax+bt throughout and obtain:

bu + ut = a2u + 2aux + uxx + (k - 1)(au + ux)- ku

Now gather like terms:

 2 ut = uxx + [2a + (k - 1)]ux + [a + (k - 1)a - k - b]u.

Then choose a = -(k - 1)/2 so that the ux coefficient is 0, and then choose b = a2 + (k - 1)a - k = -(k + 1)2/4 so the u coefficient is likewise 0. With this choice, the equation is reduced to

ut = uxx.

We also need to transform the initial condition too. This would be

u(x, 0) = e-(-(k-1)/2)x-(-(k+1)2/4)0v(x,0) ((k- 1)/2)x x = e max(e - 1,0) = max(e((k+1)/2)x - e((k- 1)/2)x,0).

For future reference, we notice that this function is strictly positive when the argument x is strictly positive, that is u0(x) > 0 when x > 0, otherwise, u0(x) = 0 for x < 0.

Now we are in the final stage, since we are ready to apply the solution representation formula:

 1 integral oo 2 u(x,t) = - V~ --- u0(s)e-(x-s)/4t ds. 2 pt - oo

However, first we want to make a change of variable in the integration, by taking z = (s - x)/ V~ --- 2t, (and thereby dz = (-1/ V~ --- 2t) dx) so that the integration becomes:

 1 integral oo V~ --- 2 u(x, t) = V~ --- u0(z 2t + x)e-z /2 dz 2p - oo

We may as well only integrate over the domain where u0 > 0, that is for z > -x/ V~ --- 2t. On that domain, u0 = e((k+1)/2)(x+z V~ -- 2t) - e((k-1)/2)(x+z V~ - 2t) so we are down to:

 integral oo V~ -- integral oo V~ - V~ 1-- ek+21(x+z 2t)e-z2/2 dz - V~ 1-- e k-21 (x+z 2t)e-z2/2 dz 2p -x/ V~ 2t 2p -x/ V~ 2t

Call the two integrals I1 and I2 respectively.

We will evaluate I1 ( the one with the k + 1 term) first. This is easy, just completing the square in the exponent yields a standard, tabulated integral. The exponent is

 V~ --- V~ --- ((k + 1)/2)(x + z 2t)- z2/2 = (- 1/2)(z2 - 2t(k + 1)z) + ((k + 1)/2)x = (- 1/2)(z2 - V~ 2t-(k + 1)z + t(k + 1)2/2) + ((k + 1)/2)x + t(k + 1)2/4 V~ ---- = (- 1/2)(z - t /2(k + 1))2 + (k + 1)x/2 + t (k + 1)2/4

Therefore

 integral integral 1 oo k+1(x+zV ~ 2t)- z2/2 e(k+1)x/2+t(k+1)2/4 oo --1(z- V~ t/2(k+1))2 V~ --- V~ --e 2 e dz = ------ V~ --------- V~ -e 2 dz 2p - x/ 2t 2p -x/ 2t

Now, change variables again on the integral, choosing y = z - V~ ---- t /2(k + 1) so dy = dz, and all we need to change are the limits of integration:

 2 integral e(k+1)x/2+t(k+1)/4 oo (-1/2)y2 V~ --- V~ - V~ --- e dz 2p -x/ 2t- t/2(k+1)

The integral can be represented in terms of the cumulative distribution function of a normal random variable, typically denoted P. That is,

 V~ --- integral d 2 P(d) = (1/ 2p) e-y /2 dy - oo

so

 (k+1)x/2+t(k+1)2/4 I1 = e P(d1)

where d1 = x/ V~ --- 2t +  V~ ---- t/2(k + 1) (note the use of the symmetry of the integral!) The calculation of I2 is identical, except that (k + 1) is replaced by (k - 1) throughout.

That, is the solution of the heat equation is

u(x,t) = e(k+1)x/2+t(k+1)2/4P(d1) - e(k-1)x/2+t(k-1)2/4P(d2)

where d1 = x/ V~ 2t- +  V~ t-/2(k + 1) and d2 = x/ V~ 2t- +  V~ t-/2(k - 1).

Now, we must systematically unwind each of the changes of variables, from u, first v(x,t) = e(-1/2)(k-1)x-(1/4)(k+1)2tu(x,t). (Notice how many of the exponentials neatly combine and cancel!) Then put x = log(S/K), t = (1/2)s2(T - t) and V (S,t) = Kv(x,t).

The final solution is the Black-Scholes formula for the value of a European call option at time T with strike price K, if the current time is t and the underlying security price is S, the risk-free interest rate is r and the volatility is s:

 ( 2 ) ( 2 ) V (S,t) = SP log(S/K)--+V( ~ r +-s-/2)(T---t) - Ke -r(T -t)P log(S/K)--+V( ~ r---s-/2)(T---t)- s T - t s T - t

Note, usually one doesn’t see the solution as this full closed form solution. Instead, most versions of the solution write intermediate steps in small pieces, and then present the solution as an algorithm putting the pieces together to obtain the final answer.

Graphical Solution of the Black-Scholes Equation

Consider for purposes of graphical illustration the value of a call option with strike price K = 100. The risk-free interest rate per year, continuously compounded is 12%, so r = 0.12, the time to expiration is T = 1 measured in years, and the standard deviation per year on the return of the stock, or the volatility is s = 0.10 . The value of the call option at maturity plotted over a range of stock prices 70 < S < 130 surrounding the strike price is illustrated below:

Value of the call option at maturity

Now, we use the Black-Scholes formula above to compute the value of the option prior to expiration. With the same parameters as above the value of the call option is plotted over a range of stock prices 70 < S < 130 at time remaining to expiration t = 1 (red), t = 0.8, (orange), t = 0.6 (yellow), t = 0.4 (green), t = 0.2 (blue) and at expiration t = 0 (black).

Value of the call option at various times

Notice a couple of trends in the value from this graph:

  1. As the stock price increases, for a fixed time the option value increases,
  2. As the time to expiration decreases, for a fixed stock value price the value of the option decreases to the value at expiration.

We can also plot the solution of the Black-Scholes equation as a function of security price and the time to expiration as value surface:

Value surface from the Black-Scholes formula


Problems to Work for Understanding

  1. Explicitly evaluate the integral I2 in terms of the c.d.f. P and other elementary functions as was done for the integral I1.

    Solution

  2. What is the price of a European call option on a non-dividend-paying stock when the stock price is $52, the strike price is $50, the risk-free interest rate is 12% per annum (compounded continuously), the volatility is 30% per annum, and the time to maturity is 3 months?

    Solution

  3. What is the price of a European call option on a non-dividend paying stock when the stock price is $30, the exercise price is $29, the risk-free interest rate is 5%, the volatility is 25% per annum, and the time to maturity is 4 months?

    Solution

  4. Show that the Black-Scholes formula for the price of a call option tends to max(S - K, 0) as t --> T.

    Solution


Reading Suggestion:

  1. Section 4.2, pages 59-63; Section 4.3, pages 66-69; Section 5.3, pages 75-76; and Section 5.4, pages 77 - 81 of The Mathematics of Financial Derivatives: A Student Introduction by P. Wilmott, S. Howison, J. Dewynne, Cambridge University Press, Cambridge, 1995.
  2. Chapter 11 of Stochastic Calculus and Financial Applications by J. Michael Steele. Springer, New York, 2001.


Outside Readings and Links:

  1. Cornell University , Department of Computer Science, Prof. T. Coleman Rhodes and Prof. R. Jarrow, Numerical Solution of Black-Scholes Equation, Submitted by Chun Fan, Nov. 12, 2002.
  2. Monash University , Department of Mathematical Science, Eric. W. Chu Applied mathematicians transform the Black-Scholes Equation into a simpler forward form and solve it.It gives some examples and maple commands, Submitted by Chun Fan, Nov. 12, 2002.
  3. ¡/a¿ An applet for calculating the option value based on Black-Scholes model. Also contains tips on options, business news and literature on options. Submitted by Yogesh Makkar, November 19, 2003.
  4. ExcelEverywhere, a commercial application for spreadsheets on the Web A sample spreadsheet based calculator for calculating the option values, based on Black-Scholes model. Submitted by Yogesh Makkar, November 19,2003


I check all the information on each page for correctness and typographical errors. Nevertheless, some errors may occur and I would be grateful if you would alert me to such errors. I make every reasonable effort to present current and accurate information for public use, however I do not guarantee the accuracy or timeliness of information on this website. Your use of the information from this website is strictly voluntary and at your risk.

I have checked the links to external sites for usefulness. Links to external websites are provided as a convenience. I do not endorse, control, monitor, or guarantee the information contained in any external website. I don't guarantee that the links are active at all times. Use the links here with the same caution as you would all information on the Internet. This website reflects the thoughts, interests and opinions of its author. They do not explicitly represent official positions or policies of my employer.

Information on this website is subject to change without notice.

Steve Dunbar's Home Page, http://www.math.unl.edu/~sdunbar1
Email to Steve Dunbar, sdunbar1@unl.edu

Last modified: [an error occurred while processing this directive]