This is machine translation

Translated by Microsoft
Mouseover text to see original. Click the button below to return to the English verison of the page.

Note: This page has been translated by MathWorks. Please click here
To view all translated materals including this page, select Japan from the country navigator on the bottom of this page.


Determine option prices or sensitivities using Leisen-Reimer binomial tree model


PriceSens = optstocksensbylr(LRTree,OptSpec,Strike,Settle,ExerciseDates)
PriceSens = optstocksensbylr(LRTree,OptSpec,Strike,Settle,ExerciseDates,Name,Value)


PriceSens = optstocksensbylr(LRTree,OptSpec,Strike,Settle,ExerciseDates) calculates option prices or sensitivities using a Leisen-Reimer binomial tree model.

PriceSens = optstocksensbylr(LRTree,OptSpec,Strike,Settle,ExerciseDates,Name,Value) calculates option prices or sensitivities using a Leisen-Reimer binomial tree with additional options specified by one or more Name,Value pair arguments.

Input Arguments


Stock tree structure created by lrtree.


NINST-by-1 cell array of character vectors 'call' or 'put'.


NINST-by-1 (European/American) or NINST-by-NSTRIKES (Bermuda) matrix of strike price values. Each row is the schedule for one option. If an option has fewer than NSTRIKES exercise opportunities, the end of the row is padded with NaNs.


NINST-by-1 matrix of settlement or trade dates.


NINST-by-1(European/American) or NINST-by-NSTRIKEDATES (Bermuda) matrix of exercise dates. Each row is the schedule for one option. For a European option, there is only one ExerciseDate on the option expiry date. For the American type, the option can be exercised on any tree data between the ValuationDate and tree maturity. The last element of each row must be the same as the maturity of the tree.

Name-Value Pair Arguments

Specify optional comma-separated pairs of Name,Value arguments. Name is the argument name and Value is the corresponding value. Name must appear inside single quotes (' '). You can specify several name and value pair arguments in any order as Name1,Value1,...,NameN,ValueN.


NINST-by-1 flags with values of 0 (European/Bermuda) or 1 (American).

Default: 0


NOUT-by-1 or 1-by-NOUT cell array of character vectors indicating the nature and order of the outputs for the function. Possible values are Price, Delta, Gamma, Vega, Lambda, Rho, and All.

Default: Price

Output Arguments


NINST-by-1 expected prices or sensitivities values.


collapse all

This example shows how to compute option prices and sensitivities using a Leisen-Reimer binomial tree model. Consider European call and put options with an exercise price of $100 that expire on December 1, 2010. The underlying stock is trading at $100 on June 1, 2010 and has a volatility of 30% per annum. The annualized continuously compounded risk-free rate is 7% per annum. Using this data, compute the price, delta and gamma of the options using the Leisen-Reimer model with a tree of 25 time steps and the PP2 method.

AssetPrice  = 100;
Strike = 100;

ValuationDate = 'June-1-2010';
Maturity = 'December-1-2010'; 

% define StockSpec
Sigma = 0.3;

StockSpec = stockspec(Sigma, AssetPrice);

% define RateSpec
Rates = 0.07;
Settle = ValuationDate;
Basis = 1;
Compounding = -1;

RateSpec = intenvset('ValuationDate', ValuationDate, 'StartDates', Settle, ...
'EndDates', Maturity, 'Rates', Rates, 'Compounding', Compounding, 'Basis', Basis);

% build the Leisen-Reimer (LR) tree with 25 time steps
LRTimeSpec  = lrtimespec(ValuationDate, Maturity, 25); 

% use the PP2 method
LRMethod  = 'PP2';  

TreeLR = lrtree(StockSpec, RateSpec, LRTimeSpec, Strike, 'method', LRMethod);

% compute prices and sensitivities using the LR model:
OptSpec = {'call'; 'put'}; 
OutSpec = {'Price', 'Delta', 'Gamma'};

[Price, Delta, Gamma] = optstocksensbylr(TreeLR, OptSpec, Strike, Settle, ... 
Maturity, 'OutSpec', OutSpec)
Price = 


Delta = 


Gamma = 



Leisen D.P., M. Reimer. "Binomial Models for Option Valuation – Examining and Improving Convergence." Applied Mathematical Finance. Number 3, 1996, pp. 319–346.

Introduced in R2010b

Was this topic helpful?