Simple Monte Carlo Options Pricer In Python

Today we will be pricing a vanilla call option using a monte carlo simulation in Python. Monte Carlo models are used by quantitative analysts to determine accurate and fair prices for securities. Typically, these models are implemented in a fast low level language such as C++. However, for the sake of ease, we’ll be using Python.

Pre-Requisites:

Below is a list of pre-requisite knowledge to get the most out of this tutorial.

Required:

• Calculus
• Probability and Statistics
• Very basic programming

Recommended:

• Stochastic Processes
• Stochastic Calculus or an introductory asset pricing class

Understanding The Math

Before diving into the code, we’ll cover some of the basic financial mathematics necessary to understand the model.

Assumptions

In our model we will make the following assumptions.

• The price of our stock will generally increase with respect to time
• The expected return is a fixed rate of the current share price
• The stock follows a random-walk behavior

The Stock Price Evolution Model

We will be using a stochastic differential equation as the model of our stock price evolution. It will consist of a component representing the expected return on the stock Sₜ over an infinitesimal period of time dt, represented by:

Where Sₜ is the price of our share with respect to time, and μ is our expected rate of return.

To make our model representative of real share price behavior we must also add a stochastic element represented by a fixed fraction of our price σ, our share price Sₜ and a random walk process dWₜ to construct

Altogether, we get the equation:

Black-Scholes Option Pricing Model

The Black-Scholes option pricing model tells us the the price of a vanilla option with a compounding rate r , expiration T and payoff function f is:

Solving for Expectation

In order to finish deriving an expression for our fair options price, we must solve the expression within the expected payoff function. To begin with, we will use a risk-neutral stochastic differential equation (see stock price evolution model) for the expectation of Sₜ:

By taking the log of both sides of the equation and using Ito’s lemma, a solution to the stochastic process:

substituting a constant coefficient for d, (a convenient model simplification) we will obtain the final expression:

Wₜ is a brownian motion process with a mean 0 and variance T and normally distributed random variable N(0,1):

Our final expression for Sₜ is:

Raising the e to the equation, we can simplify to:

We will continuously sample random numbers for N(0,1) into the equation while maintaining a running sum:

Where K is the strike price. The average of this expression where n is the number of random samples we take will be

Final Expression

Our final expression for the fair price of the option will be the above expression multiplied by e⁻ʳᵗ to complete our original black Scholes model.

The Code

import math
import random

class SimpleMCPricer():
    def __init__(self, expiry, strike, spot, vol, r, paths):
        #The sigma value on the left side of the exponent
        self.variance = vol**2 * expiry
        #The sigma value on the right side of the e exponent
        self.root_Variance = math.sqrt(self.variance)
        #Corresponds to the (-1/2 * sigma^2)
        self.itoCorr = -0.5*self.variance
        ##Corresponds to S0e^(rT - 1/2 sigma^2T)
        self.movedSpot = spot*math.exp(r*expiry + self.itoCorr)
        self.runningSum = 0
        ##Simulate for all paths
        for i in range(0,paths):
            thisGauss = random.randrange(0,1000,1)
            thisGauss = thisGauss/1000
            ##Our rootVariance already has been multiplied by the expiry
            thisSpot = self.movedSpot*math.exp(self.root_Variance*thisGauss)
            #Determine payoff of this specific path
            thisPayoff = thisSpot - strike
            #Value of option is zero is our price is less than the strike
            thisPayoff = thisPayoff if thisPayoff > 0 else 0
            self.runningSum+=thisPayoff
        
        self.mean = self.runningSum/paths
        self.mean*= math.exp(-r * expiry)
    
    def getMean(self):
        return round(self.mean,2)

Run the Model!

Let us run the model on an option with expiration in 2 years, with a strike price of 32 dollars, a current price of 30 dollars, a 10% volatility parameter, and a 3% rate of return. We will simulate 1,000,000 paths and determine the fair price.

model = SimpleMCPricer(2,32,30,.1,0.03,1000000)
model.getMean()

As you can see, the calculated fair price of the option is 1.79 dollars.

1.79

Sources and Further Reading

[1] C++ Design Patterns and Derivatives Pricing, Mark S. Joshi

[2] The Concepts and Practices of Mathematical Finance, Mark S. Joshi

[3] Stock price modelling: Theory and Practice, Abdelmoula Dmouj

Originally published at https://www.pranav.ai.