The Barone-Adesi Whaley Formula to Price American Options Revisited

This paper presents a method to solve the American option pricing problem in the Black Scholes framework that generalizes the Barone-Adesi, Whaley method [1]. An auxiliary parameter is introduced in the American option pricing problem. Power series expansions in this parameter of the option price and of the corresponding free boundary are derived. These series expansions have the Baroni-Adesi, Whaley solution of the American option pricing problem as zero-th order term. The coefficients of the option price series are explicit formulae. The partial sums of the free boundary series are determined solving numerically nonlinear equations that depend from the time variable as a parameter. Numerical experiments suggest that the series expansions derived are convergent. The evaluation of the truncated series expansions on a grid of values of the independent variables is easily parallelizable. The cost of computing the n-th order truncated series expansions is approximately proportional to n as n goes to infinity. The results obtained on a set of test problems with the first and second order approximations deduced from the previous series expansions outperform in accuracy and/or in computational cost the results obtained with several alternative methods to solve the American option pricing problem [1]-[3]. For example when we consider options with maturity time between three and ten years and positive cost of carrying parameter (i.e. when the continuous dividend yield is smaller than the risk free interest rate) the second order approximation of the free boundary obtained truncating the series expansions improves substantially the Barone-Adesi, Whaley free boundary [1]. The website: http://www.econ.univpm.it/recchioni/finance/w20 contains material including animations, an interactive application and an app that helps the understanding of the paper. A general reference to the work of the authors and of their coauthors in mathematical finance is the website: http://www.econ.univpm.it/recchioni/finance.


Introduction
American call and put options are one of the most traded products in financial markets.They are traded either standing alone or embedded in a variety of financial contracts such as, for example, convertible bonds, mortgages or life insurance policies.The fast and accurate evaluation of American option prices and of the corresponding free boundaries is an important problem in mathematical finance.Let us restrict our attention to the American option pricing problem in the Black Scholes framework.Many methods have been suggested to solve this problem.In particular several hybrid methods have been suggested.These methods combine analytical and numerical approximations.For example let us mention the hybrid methods proposed by Geske, Johnson (1984) [4], Barone-Adesi, Whaley (1987) [1], Kim (1990) [5], Bunch, Johnson, (1992) [6], Bjerksund, Stensland (1993) [7], Ju, Zhong (1999) [2], Barone-Adesi (2005) [8] and Zhu (2006) [9].
In [1] Barone-Adesi and Whaley write the American option price as the sum of the price of the corresponding European option and of a quantity called early exercise premium.The European option price is given by the Black Scholes formula and the early exercise premium is approximated with the solution of a free boundary value problem for an ordinary differential equation.This ordinary differential equation is obtained dropping the time derivative term in the partial differential equation satisfied by the early exercise premium.Barone-Adesi and Whaley [1] give a simple formula for the solution of this free boundary value problem for an ordinary differential equation.Moreover they determine an approximation of the free boundary solving numerically a nonlinear equation that depends from the time variable as a parameter.This approximate solution of the American option pricing problem is called Barone-Adesi, Whaley formula and is widely used in the financial markets by practitioners.An exhaustive review of the methods used to solve the American option pricing problem and of the developments of the Barone-Adesi, Whaley method during the period 1987-2005 can be found in Barone-Adesi (2005) [8].For example in 1999 Ju, Zhong [2] reconsidered the Barone-Adesi, Whaley formula of the early exercise premium.The Ju, Zhong formula [2] introduces a correction to the Barone-Adesi, Whaley approximation of the early exercise premium.This correction consists in writing the early exercise premium as the product of the Barone-Adesi, Whaley early exercise premium times a time-independent function determined solving an ordinary differential equation.When long dated options are considered, the Ju, Zhong formula improves the approximate option price obtained with the Barone-Adesi, Whaley formula.
Given a positive integer n, Geske and Johnson [4] approximate the price of an American put option using an n-fold compound option.They assume that exercise decisions are taken only at some known time values.These time values are a set of n points.In [4] Geske and Johnson deduce a formula to approximate the American put option price with a piecewise solution of the Black Scholes partial differential equation subject to boundary conditions imposed at the decision times.Moreover, using Richardson extrapolation, they show how to approximate the Geske, Johnson formula with a simple polynomial expression.Bunch and Johnson [6] refine the results obtained in [4] determining the n exercise times that maximize the accuracy of the option prices obtained.
In [7] Bjerksund and Stensland approximate the solution of the American option pricing problem assuming a flat early exercise boundary and using a trigger price.Bjerksund and Stensland reduce the evaluation of an American call option with exercise price E and maturity time T to the evaluation of a European call up-and-out barrier option with knock-out barrier X, strike price E and maturity time T. A rebate given by X E − is received by the holder of the option at the knock-out time when the option is exercised prior to maturity time.The barrier X is the flat boundary that approximates the free boundary of the American option pricing problem.In [7] the problem of choosing X is studied.In [10] the approximation of the free boundary used in [7] is refined.In fact in [10] the time interval where the problem is studied is divided in two disjoint subintervals and a flat early exercise boundary is used in each subinterval.Zhu (2006) [9] considers the American put option pricing problem and derives an explicit formula of the American put option price associated to a numerically approximated free boundary.This formula is a Taylor's series expansion with infinitely many terms.Each term of this Taylor's expansion considered contains several integrals that must be evaluated numerically.In [11] I. J. Kim, Jang, K. T. Kim show that the numerical evaluation of Zhu's formula is cumbersome and suggest a method to approximate the free boundary of the American option pricing problem.This method consists in the numerical solution of the integral equation satisfied by the free boundary deduced in [3] by Little, Pant, Hou.The solution of the American put option pricing problem suggested in [11] combines the integral formula of the option price obtained by I. J. Kim [5] with the approximation of the corresponding free boundary obtained solving numerically the integral equation presented in [3].Note that in the option price formula contained in [5] there are several integrals that must be evaluated numerically.
To solve the American option pricing problem instead of using hybrid methods it is possible to use only numerical methods.For example the finite differences method (see [12]), the Monte Carlo method (see [13]- [18]), and the regression method (see [19] [20]) can be used to solve the American option pricing problem.
Usually hybrid methods are computationally cheaper than numerical methods.However in many circumstances numerical methods provide approximate solutions of the American option pricing problem that are more accurate than those obtained with hybrid methods.In fact, at least in principle, the solutions provided by numerical methods can be made arbitrarily accurate choosing appropriately the values of the parameters that define the approximation computed.Instead many hybrid methods have a certain accuracy that depends from the problem under consideration and this accuracy cannot be changed choosing parameter values.That is most of the solutions found with hybrid methods do not converge to the exact solution of the American option pricing problem when a suitable limit is taken.Moreover most hybrid methods give satisfactory results when pricing problems with short maturity times are considered.The results obtained with these methods deteriorate when problems with medium or long maturity times are considered.
This paper presents a hybrid method to solve the American option pricing problem.We introduce an auxiliary parameter in the American option pricing problem and we deduce power series expansions in this parameter of the option price and of the corresponding free boundary.Explicit formulae (depending from the free boundary) are given for the coefficients of the option price series.The partial sums of the free boundary series are determined solving numerically nonlinear equations that depend from the time variable as a parameter.These series expansions are a formal solution of the American option pricing problem.Numerical experiments suggest that the series obtained are convergent.The zero-th order term of the series expansions is the Barone-Adesi, Whaley solution of the American option pricing problem [1] (i.e. the Barone-Adesi, Whaley formula).The first order approximation of the option price deduced from the expansions developed here has some similarities with the early exercise premium formula suggested by Ju, Zhong [2].
Test problems taken from [1] [2] and [3] are studied.The behaviour of the truncated series expansions on these test problems is studied.In particular in the numerical experiments presented we use the n-th order approximate solutions deduced from the expansions when n = 0, 1, 2 to solve the test problems considered.These experiments show that each approximation order of the solution deduced from the expansions adds roughly one correct significant digit to the results obtained.Moreover for n = 0, 1, ••• the computation of the n-th order approximation deduced from the expansions of the solution of the American option pricing problem on a grid of values of the independent variables is easily parallelizable and its computational cost is "substantially" linear in n as n goes to infinity.In particular the numerical experiments show that when we consider options with intermediate maturity times (i.e.: maturity times ranging in the interval 3 -10 years) the first and the second order approximations of the solution obtained from the series expansions improve substantially the approximate solution obtained using the Barone-Adesi, Whaley formula (see in Section 4, Table 1, Table 3, Table 4 and Figure 2).For example the improvement obtained with the higher order terms of the expansions is significant when we compare the approximations of the free boundary of the American option pricing problem obtained using the Barone-Adesi, Whaley formula with those obtained using the n-th order truncated power series expansions, n = 1, 2 (see Section 4, Table 1, Table 4 and Figure 2).Note that the Barone-Adesi, Whaley formula gives excellent results when we consider options with short or with long maturity times and that in these circumstances there is no room for improvements of practical value (see [1]).
The website: http://www.econ.univpm.it/recchioni/finance/w20contains material including animations, an interactive application and an app that helps the understanding of the paper.More general references to the work of the authors and of their coauthors in mathematical finance are available in the website: http://www.econ.univpm.it/recchioni/finance.
The paper is organized as follows.In Section 2 we formulate the American call option pricing problem in the Black Scholes framework and we introduce the auxiliary parameter that is used to solve it.In Section 3 we deduce the perturbation expansions in this auxiliary parameter of the American call option price and of the corresponding free boundary.The analysis of Sections 2 and 3 can be easily extended from the case of the American call option pricing problem to the case of the American put option pricing problem.This extension is omitted for simplicity.In Section 4 we present the results obtained with the method developed in Sections 2 and 3 on a set of test problems involving American call and put options.These results are compared with those discussed in the scientific literature obtained with some alternative methods to solve the American option pricing problem.

The American Option Pricing Problem in the Black Scholes Framework
We follow Barone-Adesi and Whaley [1] and we consider the problem of pricing American call and put options on commodities in the Black Scholes framework.
Let t be a real variable that denotes time and S t , t > 0, be a real stochastic process that models the commodity price, that is for t > 0 the random variable S t represents the commodity price at time t.We assume that under the risk neutral measure the commodity price satisfies the following stochastic differential equation: where b, σ are real parameters, z t , t > 0, is the standard Wiener process such that 0 0 z = and dz t is its sto- chastic differential.Equation ( 1) is known as Black Scholes asset price equation.The parameter 0 σ > is the volatility or instantaneous standard deviation and b is the cost of carrying parameter.In the most common situations we have b r d = − where r > 0 is the risk free interest rate and d > 0 is the continuous dividend yield, see [1].When needed Equation ( 1) is equipped with an initial condition.
To keep the exposition simple we study only the American call option pricing problem.The American put option pricing problem can be studied analogously.However in the test problems presented in Section 4 the method developed here to solve the American option pricing problem is used to evaluate both call and put options.
Let t = 0 be the current time, consider the problem of pricing an American call option having exercise price E > 0 and maturity time T > 0 written on a commodity whose price S t , t > 0, satisfies (1).The price ( ) , 0 t T < < , of this option and the corresponding free boundary ( ) f S t , 0 t T < < , solve the following problem [21]: 0, 0 , 0 , 2 with boundary conditions: , , 0 , and final condition: Problem ( 2), ( 3), ( 4), ( 5), ( 6) is the American call option pricing problem in the Black Scholes framework.It is a free boundary value problem for the partial differential Equation (2) whose unknowns are: the option price ( ) , 0 t T < < , and the free boundary ( ) f S t , 0 t T < < , see [22].Note that the boundary condition (3) can be omitted.In fact condition (3) follows from the degeneracy in 0 S = of (2) and from the fact that in mathematical finance only bounded solutions of (2) are meaningful (see [21] Chapter 3, p. 48-49).However to make easier the understanding of some choices made later (see Section 3) we prefer to state (3) explicitly.
Let us consider the change of variable: T t τ = − , 0 t T ≤ ≤ , the variable τ is called time to maturity, and let us define: 3), (4), ( 5), (6) rewritten in the variables S, τ for the unknowns C A , S * becomes: ( ) with boundary conditions: and initial condition: ) , be the Black Scholes price of the European call option having strike price E, maturity time T and parameters r, b, σ .The price ( ) , has a simple expression given by the Black Scholes formula [21].Note that C E is defined for 0 S > , 0 T τ < < .As done in [1] we seek a solution of problem ( 7), ( 8), ( 9), (10), (11) given by the sum of ( ) Substituting ( 12) in ( 7), ( 8), ( 9), ( 10), (11) and using the fact that C E satisfies the Black Scholes partial differential Equation (7), the boundary condition (8) and the initial condition (11) we obtain the following problem: with boundary conditions: , , , 0 , and initial condition: Problem ( 13), ( 14), ( 15), ( 16), ( 17) is a free boundary value problem for the partial differential Equation (13) in the unknowns ( ) We assume that the early exercise premium e has the following form (see [1]): where , is a sufficiently regular function that will be chosen later and , is an auxiliary unknown that must be determined solving ( 13), ( 14), ( 15), ( 16), ( 17), (18).Substituting (18) in (13) we have: where Equation ( 18) and the previous choice of K imply that the boundary conditions ( 14), ( 15), ( 16) can be rewritten respectively as: Note that in the formulation of problem ( 20), ( 21), ( 22), ( 23) we use the two variables τ and K, and recall that these two variables are linked by the condition This implies that when the function f is well behaved in 0 τ = the function e defined in (18) satisfies the initial condition (17).
Let us recall that in [24] a similar approach has been used in the study of barrier options.In fact in [24] it is considered the problem of pricing (put up-and-out) barrier options with time-dependent parameters in the Black Scholes framework.An auxiliary parameter is introduced in the barrier option pricing problem and a perturbation expansion in this parameter of the barrier option price is deduced.Note that the perturbation problem studied in [24] is a regular perturbation problem, while the perturbation problem considered here when 0 + →  is a singular perturbation problem.

A Series Expansion of the Solution of the American Option Pricing Problem
Let us drop the initial condition (28) from problem ( 24), ( 25), ( 26), ( 27), (28), (29).That is let us consider the equation: with the boundary conditions: Recall that once determined f  , S *  as solution of (30), (31), (32), (33) we will recover the function e using (29) and that when f  is well behaved in 0 τ = the function e determined in this way will satisfy the initial condition (28) as a consequence of the choice Let us as- sume that the following expansions in powers of  of the functions f  , S *  , 0 1 ≤ ≤  , hold: where the functions ˆj f , ˆj S * , 0,1, j =  are independent of  .For later convenience we define the partial sums ,  , of the series (34), (35), that is: , Note that for 0,1, j =  the problems that follow define the function ˆj f for ( ) To give a meaning to the sums contained in (34) and (36) we extend with zero the function ˆj f when ( ) With abuse of notation ˆj f denotes both the original function and the extended function 0,1, j =  .We impose (30), (31), ( 32 .This is what has been done by Barone-Adesi and Whaley [1] in the solution of the zero-th order problem.We simply extend the idea of Barone-Adesi and Whaley from the zero-th order problem to the n-th order problem, n = 1, 2, •••.Proceeding in this way we obtain: For n = 0 the zero-th order problem is: with boundary conditions: for n = 1, 2, ••• the n-th order problem is: with boundary conditions: The problems (38), (39), (40), ( 41) and ( 42), ( 43), (44), (45) are respectively free boundary value problems for the ordinary differential Equations ( 38) and (42).These problems depend from the parameter τ , 0 T τ < < .The boundary conditions (40), ( 41) and ( 44), (45) are the boundary conditions imposed in 1 =  derived from (32), (33).Note that we have: 0 where f 0 , 0 S * is the Barone-Adesi, Whaley solution of the American option pricing problem deduced in [1].Numerical experiments have shown that the approximations of the first few orders of f and S * deduced from the series expansions (34), (35) with the boundary conditions (40), ( 41) and ( 44), (45) imposed to the partial sums (37) evaluated in 1 =  are of higher quality than the corresponding approximations of the same order in  obtained imposing the boundary conditions (32), (33) order by order in powers of  in 0 =  .In particular these numerical experiments show that the approximations of , obtained solving problem (42), ( 43), ( 44), (45) are more accurate than the corresponding approximations obtained evaluating in 1 =  the partial sums , where ˆi S * , 1, 2 i = are the terms obtained imposing in 0 =  the boundary conditions (32), (33) at the first and at the second order in powers of  .
From now on instead of using the notation 0 f , 0, 1 S * =  introduced previously we denote the solution of the zero-th order problem (38), (39), ( 40), (41) with the notation f 0 , 0 S * .This choice emphasizes the fact that the zero-th order term of the expansions (34), (35) determined solving the previous problems is the Barone-Adesi, Whaley solution of the American option pricing problem.As done in [1] we seek a solution f 0 of (38), (39), ( 40), (41) of the following form: ( where in (46) the functions 0,0 A and q are auxiliary unknowns that must be determined.Substituting (46) in (38) we have: Equation ( 47) is satisfied if we impose that: the quadratic Equation (48) in the unknown q is easily solved, and one of its solutions is: From (49) it follows that ( ) q K is positive when K is positive.This means that when ( ) q k is given by (49) the function given by ( 46) satisfies (39).The second solution of (48) is negative when K is positive and must be discarded since the corresponding function f 0 given by (46) does not satisfy (39).
Substituting the formulae (46), (49) in the Equations ( 40), (41) we obtain respectively: ( ) that depends from the parameter τ .Given τ this last equation is easily trans- formed in a fixed point problem and solved numerically using Banach iteration, 0 T τ ≤ ≤ .The approximate solution of (51) determined in this way is substituted in (50) to obtain 0,0 A .The zero-th order term f 0 given by ( 46), ( 49), (50) with the numerical approximation of 0 S * substituted in 0,0 A multiplied by K (see (18)) is the Baroni-Adesi, Whaley formula of the early exercise premium (see [1] formula (20)).The numerical approximation of 0 S * obtained solving numerically (51) is the Barone-Adesi, Whaley approximation of the free boundary (see [1], formula (19)).Note that as already said with abuse of notation in the previous formulae ( ) 0 S τ * , 0 T τ < < , denotes both the unknown of the nonlinear Equation (51) and its numerical approximation.This ambiguity is reflected in the functions 0,0 A and f 0 .Let 1, 2, n =  we seek a solution ˆn f of the n-th order problem (42), ( 43), ( 44), (45) of the following form: , where the functions  , are auxiliary unknowns that must be determined impos- ing (42), ( 43), ( 44), (45).Substituting formula (52) in Equation (42) it is easy to see that in order to satisfy (42) it is sufficient to impose that the functions  , satisfy the following systems of linear equations: , .
To keep the notation simple in (53), ( 54), (55) we have omitted the dependence from K of the functions , are determined solving the n-th linear system of linear equations contained in (53), (54), (55).In fact for 1, 2, n =  the n-th system of li- near equations contained in (53), (54), (55) is a system of 2n linear equations in the 2n unknowns , that can be solved by backward substitution starting from the n-th Equation (53).Note that due to its special form the computational cost of solving these linear systems obtained in (53), ( 54 , . To keep the notation simple in (56), (57) we have omitted the dependence from τ of , 1 and of the solution of the n-th linear system contained in (53), ( 54), (55), 0 T τ < < .Equation ( 57) is a nonlinear equation in the unknown ( ) , depending from the parameter τ , that defines implicitly ( ) In the numerical experiments of Section 4 this equation is transformed in a fixed point problem and solved numerically using Banach iteration.The numerical approximation of ( )  , obtained from the solution of (57) is substituted in (56) and together with the solution of the n-th linear system contained in (53), ( 54), ( 55) determines (56).Note that when 1, 2, n =  , for simplicity with abuse of notation we denote with ( ) , both the unknown of (57) and its numerical approximation.A careful inspection of formulae ( 46), ( 51) and ( 52), (57) shows that for 0,1, n =  the numerical evaluation on a grid of values of the S and τ variables of the n-th order approximations ˆn f , of the solution of the American call option pricing problem is easily parallelized.

Numerical Results
Let us discuss the numerical results obtained on a set of test problems with the solution method of the American option pricing problem developed in Sections 2 and 3.
We use the trinomial tree method [17] with n T = 1000 time steps to compute the "true value" of the option prices considered in our experiments.The choice n T = 1000 guarantees four correct significant digits in the option prices computed in this Section.The "true value" of the corresponding free boundaries of the American call options considered in our experiments is computed solving numerically the following integral equation (see [25], [3] and the reference therein): where in (58), ( 59), (60) T is the maturity time and E is the strike price of the American call option considered.The integral operator contained in (58) is approximated with the composite rectangular rule with time step 0 τ ∆ > and ( ) be the solution of (58) it is easy to see that we have: When we consider an American call option and we have b r < the algorithm that solves the discretized version of (58) starts from 0 Recall that when b r ≥ (i.e. when the continuous dividend yield d is smaller or equal to zero) in [1] page 307 it is shown that the American call option price reduces to the corresponding European call option price.In fact in this case the free boundary is "at infinity" and the American call option must be exercised at maturity time, that is the value of the early exercise premium is identically zero.
The integral Equation (58) must be modified to deal with American put options (see [3]).Moreover in the case of put options the algorithm that solves the corresponding discretized integral equation starts from 0 In the numerical experiments discussed below the discretized versions of Equation (58) for call options and of the analogous equation for put options (see [3]) are solved iteratively at the points . That is we consider the unknowns ( ) , , and 0, a S τ ∆ is assigned as specified above.Of course when we consider put options Equation ( 58) must be substituted with a different integral equation, see [3], and as a consequence Equation (62) must be modified coherently.Equations ( 61), (62) (and their analogous for put options) are solved using Banach iteration.For 1, 2, , m ν =  and 0,1, ∆ be the j-th element of the sequence generated by Banach iteration applied to (61), (62), 1, 2, , m ν =  .The Banach iteration associated to is stopped at the smallest value of the index j that satisfies the condition: .
Note that in general the stopping value of the index j defined by (63) depends from ν , 1, 2, , m ν =  .The "true values" of the option prices and of the corresponding free boundaries defined previously are used as benchmarks to test the approximate solutions of the American option pricing problem computed with the method developed in Sections 2 and 3 and with some alternative methods taken from the scientific literature.
We begin our numerical experiments studying some test problems taken from [1] Section C.4.These test problems consider options on long-term U.S. Treasury bonds (time to maturity up to three years) and long term care insurance inflation options (time to maturity up to ten years and beyond).
In the first experiment we use the values of the Black Scholes parameters of Table V in [1].That is we consider the following three sets of parameter values: Table 1 shows several free boundary approximations of the American call option pricing problems specified above when 3 T τ = = , 5, 10.From left to right Table 1 shows the free boundary

S T *
is the "true free boundary" used as benchmark) and the approximations ( ) , discussed in Sections 2 and 3. Recall that 0, 1 0

S S * * = =
 is the Barone-Adesi Whaley approximation of the free boundary.Table 1 shows that for 0,1 n = going from the n-th order approximation to ( ) 1 n + -th order approximation of the free boundary roughly adds one correct significant digit to the approximation of the free boundary found.This effect is particularly evident when 0 b > , in fact in this case the Barone-Adesi, Whaley approximation of the free boundary is poor and has no correct significant digits.Note that positive values of b correspond to values of the continuous dividend yield d r b = − smaller than the risk free interest rate r.Recall that for American call options when 0 b > and 0 r b − > the integral Equation (58) becomes singular as 0 τ + → , see [25].When T τ = (i.e. when 0 t = ) Table 2 shows the asset price S, the European call option price C E obtained evaluating the Black Scholes formula, the approximations of the American call option price obtained using the trinomial tree method C T (i.e.C T is the "true value" of the option price used as benchmark), the Barone-Adesi, Whaley option price C BW obtained from ( 12), (18) when  12), (18) and the series expansions presented in Sections 2 and 3 truncated after . Furthermore Table 2 shows the relative errors: . Table 1 and Table 2 suggest that increasing the approximation order of the solution of the American call option pricing problem that has been deduced from the expansions in powers of  developed in Sections 2 and 3 (that is increasing n) it is possible to improve substantially the results obtained with the Barone-Adesi, Whaley formula (i.e. the result obtained when 0 n = ).Figure 1 shows the "true" free boundaries of the American call option pricing problem as a function of τ that have been obtained solving numerically the discretized version of the integral Equation ( 58) when E = 100,          , and b = −0.04,or b = 0.0, or b = 0.04, and 0 T τ < < .Figure 2 shows the "true" free boundary, the Barone-Adesi, Whaley free boundary, the first and the second order approximations of the free boundary obtained from the expansions of Sections 2 and 3 as a function of the time to maturity τ , 0 T τ < < , when E = 100, T = 10, r = 0.08, 0.2 σ = , b = 0.04.In particular Figure 2 shows that the first order approximation of the free boundary improves significantly the zero-th order approximation of the free boundary (i.e. the Barone-Adesi, Whaley free boundary) and that the second order approximation of the free boundary refines the result obtained with the first order approximation.The approximated option prices are obtained using (12), (18) and the zero-th, first and second order approximations of the sum the series expansions developed in Sections 2 and 3. Figure 2 and Figure 3 suggest that the second order corrections of the option price and of the free boundary are necessary to obtain satisfactory approximations of the solution of the American call option pricing problem when we consider  Figure 4 shows the American call option price C when T τ = (i.e. when t = 0), T = 10, E = 100, r = 0.08, 0.2 σ = , b 0.04 as a function of the asset price S. Note that the option price approximation obtained using the zero-th order term (i.e. the Barone-Adesi, Whaley price) is not accurate (see Figure 4 dotted line).This is a consequence of the fact that the zero-th order approximation of the free boundary is unsatisfactory (see Figure 2, Figure 4).Note that in Figure 4 the payoff function (solid line) is a lower bound for the "true" American call and the values of the following Maximum Absolute Errors (MAE): where in (68), (69) i takes the values 1, 2, ,5  and j takes values 1, 2, , 4  .Table 5 of [2] compares on the test problems considered the accuracy of the put option prices computed with some well known methods used to solve the American put option pricing problem including the methods of Geske and Johnson [4], Bunch and Johnson [6], Brodie and Detemple [19] (see [2] for further details).
shows that the root mean square errors and the maximum absolute errors of the put option prices obtained using the second order approximation of the solution of the American put option pricing problem deduced from the analogous for put options of ( 12), ( 18) and of the expansions introduced in Sections 2 and 3 outperform those shown in Table 5 of [2].In particular note that the second order approximation of the put option price obtained with the method developed in Sections 2 and 3 outperforms in accuracy the best approximation of the American put option price obtained with the Brodie and Detemple method shown in Table 5 of [2].Moreover Table 3 shows that the behaviour of the series expansions developed in Sections 2 and 3 applied to the put option pricing problem is similar to their behaviour in the case of the call option pricing problem shown in Table 2.In fact in Table 3 when n = 0, 1 going from the n-th order approximation to the ( ) 1 n + -th order ap- proximation of the American put option price we gain roughly one correct significant digit in the approximate solution found.A similar behaviour has already been observed in the call option case in Table 2.
The last test problem studied is taken from Little, Pant, Hou [3].We consider the American put option pricing problem when T = 1, E = 100, r = 0.07, 0.25 = − = (i.e.b = 0.04).We compare the free boundaries obtained using the approximations shown in Table 5 and Table 6 of [3] and the first three order approximations, ( ) S τ * , 0 T τ < < , denotes the "true" free boundary of the American put option pricing problem obtained solving numerically the integral equation for the put free boundary presented in [3].Note that the second order approximation of the free boundary obtained with the method developed in Sections 2 and 3 (adapted for put options) (shown in Table 4) is roughly of the same quality of the numerically computed free boundary obtained in [3] solving the integral equation analogous to (58) satisfied by the free boundary of the American put option pricing problem (see Table 4, Table 5 of [3]).
Finally let us compare the computational times needed to obtain the approximations of the solution of the American option pricing problem that have been considered in this Section.Let us consider the American call option pricing problem defined by S = 70,  5 shows the time (in seconds) required to compute for n = 0, 1, 2 the coefficients required to compute the same option price and the corresponding free boundary with the trinomial tree method when n T = 1000 steps and solving numerically the integral equation satisfied by the free boundary.These computational times have been obtained using an Intel CORE 3i processor.Table 5 shows that the time required to compute the first three order approximations of the solution of the American call option pricing problem (n = 0, 1, 2) deduced from ( 12), (18) and the series expansions of Sections 2, 3 is negligible when compared to the time required to compute the same quantities (price and free boundary) with the trinomial tree method and the numerical solution of the integral equation.Furthermore, as already said, it is easy to see that for 0,1, n =  the computation on a grid of values of the S and τ variables of the n-th order approximation of the option price and of the corresponding free boundary derived from ( 12), (18) and the series expansions introduced in Sections 2, 3 can be easily parallelized.In fact the computation of the free boundary at a given order on a grid of maturity times can be done in parallel since at each order in  the condition that defines the free boundary in the method developed Sections 2 and 3 is a local condition.Furthermore for any fixed maturity time the evaluation of the two series (called inner and outer expansions).Let ( ) O ⋅ be the Landau symbol, roughly speaking in the K variable the inner expansion of the option price holds when ), (33) to the series expansions (34), (35), order by order in powers of  .For n = 0, 1, ••• the unknowns of the n-th order problem are the functions ˆn f , that as unknown of the n-th order problem, we use S * , n = 0, 1, •••.This is due to the fact that since we are interested only in the solution of problem (30), (31), (32), (33) when 1 =  we impose the boundary conditions (32), (33) in 1 =  instead of imposing them in 0 =  Equation (50) defines 0,0 A as a function of ( ) 0 S τ * and Equation (51) is a nonlinear equation in the unknown ( ) 0 S τ * ), (55) grows linearly in n when n goes to infinity.Finally for 1boundary conditions (44), (45), that is imposing respectively:


is the unknown of the integral Equation (58), moreover in (58) the function is the standard normal cumulative distribution and the functions d and d ξ are given by: approximated on the set of evenly spaced nodes of step 0 τ ∆ > used to discretize the integral operator.This discretized version of the integral Equation (58) is solved as a fixed point problem using Banach iteration.Let     denote the integer part of •, for 0 of the discretized version of the integral Equation (58) at the node ν set of Black Scholes parameters we consider three call options with strike price 100 E = and maturity time 3 T = , or 5 T = , or 10 T = .The maturity time T is expressed in years.
the discretized version of the integral Equation (58) (i.e.

Figure 3 (
a) and Figure 3(b) show the "true" and the approximated prices of the American call option corresponding to the free boundaries shown in Figure 2 as a function of the time to maturity τ , 0 T τ < < , T = 10, when the asset price S takes the values S = 90 (Figure 3(a)) and S = 110 (Figure 3(b)).

Figure 2 .
Figure 2. American call option free boundary S * as a function of the time to maturity τ when E = 100, T = 10, r = 0.08, σ = 0.2, b = 0.04: "true" free boundary (solid line), Barone-Adesi, Whaley free boundary (dotted line), first order approximation of the free boundary (dashed line), second order approximation of the free boundary (dash-dotted line).

Figure 3 .
Figure 3.American call option price C as a function of the time to maturity τ for two values of the asset price S = 90 (a), S = 110 (b), when E = 100, T = 10, r = 0.08, σ = 0.2, b = 0.04."True" price C T (solid line), Barone-Adesi, Whaley price C BW (dotted line), first order approximation ,1, 1 A C =  of the price (dashed line), second order approximation ,2, 1 A C =  of the price of the free boundary obtained with the method developed in Sections 2 and 3 (adapted for put options) shown in Table4.In Table4


, and the time (in seconds)

Table 4 .
Approximations of the free boundary of an American put option with maturity T = 1 and strike price E = 100.

Table 1 .
Approximations of the free boundary of an American call option with intermediate maturity T and strike price E = 100.

Table 2 .
Approximations of the price of an American call option with intermediate maturity T, strike price E = 100 and negative cost of carrying.

Table 3 .
Approximations of the price of an American put option with T = 3, E = 100, r = 0.08, σ = 0.2 and relative and absolute errors committed.