Black-Scholes Model - Maple Programming Help

Online Help

All Products    Maple    MapleSim

Home : Support : Online Help : Applications and Example Worksheets : Applications : applications/BlackScholes

Black-Scholes Model

In this application, we compute the option price using three different methods. The first method is to derive the analytical solution to the option price based on the classical Black-Scholes model. Next, we compute the option price through Monte Carlo simulation based on the Black-Scholes model for stock price estimation.  Finally, we use the Black-Scholes differential equation model to estimate the option price.

Model Overview

Analytic Solution

Monte Carlo Simulation

Differential Equations

Overview of the Model

We consider the classical Black-Scholes model with single risky asset that follows a geometric Brownian motion:




where (Wt,t0) is a standard Brownian motion, σ0 is the constant volatility, r0 is the constant risk-free rate and S00 is the initial asset price.


Under these conditions, for any t0, the stock price St is given by the following formula:



We consider a security with time to maturity T and the payoff function:



Payoff of the form PS={1S>K0SK corresponds to a digital call options with strike price, K.

We will consider several methods for computing the price of this security.


r   =  

σ  =  

S0 =  

T  =  

K  =  

Analytic Solution


UStatistics:-RandomVariableLogNormalrσ22T, σT:

ST can be represented in the form:


where Φ is a lognormal random variable with parameters Trσ22 and σT.


The price of this option can be computed as the discounted expected payoff of the option:




P_analyticPST , K



V_analytic&ExponentialE;rTStatistics:-ExpectedValueP_analytic assuming r &gt; 0&comma; &sigma; &gt; 0&comma; S0&gt;0&comma; K&gt;0&comma; T&gt;0



Analytic Price


We can use the analytic result to study the various market sensitivities. For example, we can symbolically compute the delta of our option.

&Delta;diffV_analytic&comma; S0



Here is a formula for the Gamma:

local Γ  factordiffΔ&comma; S0



We can also use the symbolic formula to plot the option price as a function of the parameters.

Monte Carlo Simulation

Alternatively, we can estimate the expectation using Monte Carlo simulation to compute the option price.


The discrete-time version of the model is:



where h&equals;1N, and &Phi; is drawn from the lognormal distribution with parameters rs22h and σh.


We can use this expression to generate a sample path for the price of our risky asset.


Simulating Stock Prices

*please be patient, it may take a few seconds to generate a sample

Number of Replications

Number of Updates




Note: We know the distribution of the final stock price. To compute the option price, we need only to simulate the final stock price, and not the whole stock path.


We can verify the above analytic result using Monte Carlo simulation.


Monte Carlo Simulation

Number of Replications

Option Price

Standard Error


Differential Equations

Finally, we can use the Black-Scholes differential equations to compute the option price.





The key boundary condition is:





Another obvious condition:





Finally, if StK for some t<T, then it holds with a high probability that STK. Our option will thus be exercised and produce cash flow P&lpar;STK&rpar;&approx;PST




Numeric PDE Solver

Space Step

Time Step

Option Price