# European Option Pricing Formula in Risk-Aversive Markets

Authors: Shujin Wu; Shiyu Wang
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

In this study, using the method of discounting the terminal expectation value into its initial value, the pricing formulas for European options are obtained under the assumptions that the financial market is risk-aversive, the risk measure is standard deviation, and the price process of underlying asset follows a geometric Brownian motion. In particular, assuming the option writer does not need the risk compensation in a risk-neutral market, then the obtained results are degenerated into the famous Black–Scholes model (1973); furthermore, the obtained results need much weaker conditions than those of the Black–Scholes model. As a by-product, the obtained results show that the value of European option depends on the drift coefficientμ of its underlying asset, which does not display in the Black–Scholes model only because μ=r in a risk-neutral market according to the no-arbitrage opportunity principle. At last, empirical analyses on Shanghai 50 ETF options and S&P 500 options show that the fitting effect of obtained pricing formulas is superior to that of the Black–Scholes model.


---

## Body

## 1. Introduction

The option pricing theory began in 1900 when the French mathematician Louis Bachelier deduced an option pricing formula under the assumption that underlying asset prices follow a Brownian motion with zero drift. Since then, lots of researchers have contributed to the theory. Black and Scholes [1] present the very famous option pricing formula (i.e., Black–Scholes model) in a risk-neutral market and according to the no-arbitrage opportunity principle. Merton [2] shows the Black–Scholes-type model can be derived from weaker assumptions than in their original formulation and present some pricing methods for non-European options. Bakshi et al. [3] first derive an option pricing model that allows volatility, interest rates, and jumps to be stochastic. Gârleanu et al. [4] model demand-pressure effects on option prices. The model shows that demand pressure in one option contract increases its price by an amount proportional to the variance of the unhedgeable part of the option. Cai and Kou [5] propose a jump diffusion model for asset prices whose jump sizes have a mixed-exponential distribution, which is a weighted average of exponential distributions but with possibly negative weights, and then they extend the analytical tractability of the Black–Scholes model to alternative models. Bernarda and Czadob [6] investigate the pricing of basket options and more generally of complex exotic contracts depending on multiple indices. Their approach assumes that the underlying assets evolve as dependent GARCH1,1 processes. The dependence among the assets is modeled using a copula based on pair copula constructions. Bandi and Bertsimas [7] combine robust optimization and the idea of ε−arbitrage to propose a tractable approach to price a wide variety of options. Bao et al. [8] present a method that there is a possibility to get statistical arbitrage from Black–Scholes’s option price.In the last five years, there are still many researchers contributing to the theory of option pricing. Moretto et al. [9] study option pricing under deformed Gaussian distributions. Leippold and Scharer [10] develop a stochastic liquidity model, and they investigate discrete-time option pricing with stochastic liquidity. Hoka and Chanb [11] develop an option pricing method based on Legendre series expansion of the density function, and approximation formulas for pricing European type options are derived. Davison and Mamba [12] obtain a solution of the Black–Scholes equation with a nonsmooth boundary condition using symmetry methods. Willems [13] derives a series expansion for the price of a continuously sampled arithmetic Asian option in the Black–Scholes setting. The expansion is based on polynomials that are orthogonal with respect to the log-normal distribution. More literature studies can refer to Liu et al. [14], Friz et al. [15], Dubinsky et al. [16], Huh [17], Liu et al. [18], Siddiqi [19] and their studies.Although there are a huge number of literature studies on option pricing, they almost assume the financial markets are risk-neutral and complete, especially since Black and Scholes [1]. However, according to the theory and empirical analysis of risk, real financial markets are risk-aversive and incomplete. That is, investors need risk compensation for risky assets, and many risky assets cannot be duplicated by any portfolio constructed in real financial markets.In this study, using the method of discounting the terminal expectation value into its initial value, we obtain European option pricing formula under the assumptions that the financial market is risk-aversive, the risk measure is standard deviation, and the price process of underlying asset follows a geometric Brownian motion. In particular, if the option writer does not need the risk compensation in a risk-neutral market, then our obtained results are degenerated into the Black–Scholes model [1]; furthermore, our obtained results need much weaker conditions than those in the Black–Scholes model. At last, we take the Shanghai 50 ETF options, the first floor option in the Chinese financial market, and S&P 500 options as samples to compare the fitting effect. The empirical analyses show that the fitting effect of our pricing formulas is superior to that of the Black–Scholes model.


## 2. The Black–Scholes Formula

In this study, we will investigate European option pricing and compare our results with those of Black and Scholes [1]. Thus, we first retell the main results of Black & Scholes [1].Black and Scholes [1] present nine assumptions in the market for the security and for the option and then obtain their famous option pricing formula.Assumption 1.
Security price satisfies a geometric Brownian motion (GBM) model, where its drift coefficient and diffusion coefficient are constant through time. That is, the security price satisfies stochastic differential equation:(1)dStSt=μdt+σdBt,where μ and σ are constant and σ>0.Assumption 2.
The short-term interest rater is known and is constant through time. That is, the risk-free bond price satisfies ordinary differential equation:(2)dPtPt=rdt,where r is constant.Assumption 3.
The security pays no dividends or other distributions.Assumption 4.
There are no transaction costs in buying or selling the security or the option.Assumption 5.
The security can be continuously transacted.Assumption 6.
The amount of security can be arbitrarily divided.Assumption 7.
It is possible to borrow any fraction of the price of a security to buy it or to hold it, at the short-term interest rate.Assumption 8.
There are no penalties to short selling. A seller who does not own a security will simply accept the price of the security from a buyer and will agree to settle with the buyer on some future date by paying him an amount equal to the price of the security on that date.Assumption 9.
There is no-arbitrage opportunity.
When the above assumptions all hold, Black and Scholes [1] derived the pricing formula for European options, which is the Black–Scholes model.Theorem 1 (see [1]).
If Assumptions1 to 9 hold, then the values of European call option and European put option follow as(3)CS0,K,r,σ,τ=S0Φd2−Ke−rτΦd1,(4)PS0,K,r,σ,τ=Ke−rτΦ−d1−S0Φ−d2,where S0 is the initial price of underlying asset, K is the strike price of option, r is the short-term interest rate, σ is the diffusion coefficient of underlying asset, τ is the left expiration time of option, Φ⋅ is the cumulative density function of standard normal distribution, and(5)d1=1στlnS0K+r−12σ2τ,d2=d1+στ.


## 3. European Option Pricing in Risk-Aversive Markets

Black–Scholes model and its modified versions have some defects. In fact, because real financial markets are incomplete, an option may not be duplicated constantly, so its value deduced by the asset duplication method and no-arbitrage principle may lose the deductive basis. On the other hand, real financial markets are risk-aversive. Option seller undertakes the total risk and option buyer has no any risk, so option seller needs a reasonable risk compensation according to the theory of risk. In this section, we will deduce the option pricing formula in risk-aversive markets only under three assumptions, i.e., Assumptions1, 2, and 4 in Section 2, which is reasonably far more than the Black–Scholes model and its modified versions.In risk-aversive markets, assume the price process of some risky assetX by Xt,t≥0. Then the value of European call option at expiration time T with underlying asset X and strike price K follows as(6)XT−K+,here and in the sequel the operator ⋅+=max0,⋅. The value of European put option at expiration time T with underlying asset X and strike price K follows as(7)K−XT+.Note that no matter call option or put option, it is always its seller undertakes the total risk and its buyer has no any risk. According to risk theory, the seller reasonably requires some risk compensationλρX, where λ≥0 is the risk-compensation coefficient, ρ⋅ is the risk measure, and X is the risk size. After the seller has received the reasonable risk compensation, the seller takes risky asset as equivalent risk-free bond, so it follows from (6) that the value at time t of European call option with underlying asset X, strike price K, and expiration time T follows as(8)CtX,K,T,r,λ=e−rT−tEXT−K++λρXT−K+,∀t≤T,where r≥0 is the risk-free rate during t,T. Analogically, it follows from (7) that the value of European put option with underlying asset X, strike price K, and expiration time T follows as(9)PtX,K,T,r,λ=e−rT−tEK−XT++λρK−XT+,∀t≤T.In conclusion, we obtain the following proposition from (8) and (9).Proposition 1.
In risk-aversive market, assume that the risk measure isρ⋅ and the risk-compensation coefficient is λ≥0, and assuming a European option with underlying asset X, strike price K, and expiration time T, then its call-option value at time t follows as(10)CtX,K,T,r,λ=e−rT−tEXT−K++λρXT−K+,and its put-option value at time t follows as(11)PtX,K,T,r,λ=e−rT−tEK−XT++λρK−XT+,where r≥0 is the risk-free rate during t,T and t≤T.In order to obtain a closed-form solution to Proposition1, in the following, we always assume the price process of underlying asset follows some geometric Brownian motion model (1), and the risk measure is the standard deviation, i.e., ρZ=stdZ for any risk variable Z.Using Proposition1, we can deduce the value of European call option in risk-aversive markets and under Assumptions 1, 2, and 4.Theorem 2.
In a risk-aversive market, assuming that underlying asset follows a geometric Brownian motion with drift coefficientμ and diffusion coefficient σ>0, the current price of underlying asset is S0, risk-free interest rate is r through the time, the risk-compensation factor is λ≥0, and the risk measure is standard deviation, then the value of European call option with strike price K and left expiration time τ follows as(12)CS0,K,r,σ,τ,λ=S0eμ−rτΦd2−Ke−rτΦd1+λe−rτ⋅sqrtS02e2μτeσ2τΦd3−Φ2d2−KS0eμτΦd2Φ−d1+K2Φd1Φ−d1,where(13)d1=1στlnS0K+μ−12σ2τ,dm=d1+m−1στ,m=2,3.The Proofof Theorem 2 refers to Appendix A.If a financial market is risk-neutral, then investors treat expected return and deterministic return equally, so expectation yieldμ equals risk-free yield r, that is, μ=r. Otherwise, if μ≠r, there will exist arbitrage opportunity. Furthermore, in a risk-neutral financial market, the risk-compensation factor equals zero, that is, λ=0. Thus, it yields from Theorem 1 that we have the following corollary.Corollary 1.
In a risk-neutral market, assuming that underlying asset follows a geometric Brownian motion with drift coefficientμ and diffusion coefficient σ>0, the current price of underlying asset is S0, and risk-free yield is r through the time, then the value of European call option with strike price K and left expiration time τ follows as(14)CS0,K,r,σ,τ=S0Φd2−Ke−rτΦd1,where(15)d1=1στlnS0K+r−12σ2τ,d2=d1+στ.In the following, we will deduce the pricing formula of European put option. Using Proposition1, we will construct the pricing formula of European put option in risk-aversive markets and under Assumptions 1, 2, and 4.Theorem 3.
In a risk-aversive market, assuming that underlying asset follows a geometric Brownian motion with drift coefficientμ and diffusion coefficient σ>0, the current price of underlying asset is S0, risk-free yield is r through the time, the risk-compensation factor is λ≥0, and the risk measure is standard deviation, then the value of European put option with strike price K and left expiration time τ follows as(16)PS0,K,r,σ,τ,λ=Ke−rτΦ−d1−S0eμ−rτΦ−d2+λe−rτ⋅sqrtS02e2μτeσ2τΦ−d3−Φ2−d2−KS0eμτΦ−d2Φd1+K2Φ−d1Φd1,where(17)d1=1στlnS0K+μ−12σ2τ,dm=d1+m−1στ,m=2,3.The Proof of Theorem3 refers to Appendix B.According to the analysis before Corollary1, if a financial market is risk-neutral, then μ=r and λ=0. Furthermore, it yields from Theorem 3 that we have the following corollary.Corollary 2.
In a risk-neutral market, assuming that underlying asset follows a geometric Brownian motion with drift coefficientμ and diffusion coefficient σ>0, the current price of underlying asset is S0, and risk-free yield is r through the time, then the value of European put option with strike price K and left expiration time τ follows as(18)PS0,K,r,σ,τ=Ke−rτΦ−d1−S0Φ−d2,where(19)d1=1στlnS0K+r−12σ2τ,d2=d1+στ.Remark 1.
AlthoughCorollaries 1 and 2 obtain the same values as those in Theorem 1 for European call option and European put option, Corollaries 1 and 2 need much weaker conditions than those of Theorem 1. That is, Corollaries 1 and 2 improve Theorem 1.In fact, if Assumptions1, 2, and 4 in Section 2 hold, Corollaries 1 and 2 hold. However, the conditions that Theorem 1 (i.e., the Black–Scholes model) holds are Assumptions 1 to 9 in Section 2.


## 4. Empirical Analysis of Shanghai 50 ETF Options

In the section, we will present empirical analysis on Shanghai 50 ETF options, the first floor option in the Chinese financial market, and use the data of September 3 and 4, 2018, to compare the fitting effect of our pricing model and the Black–Scholes pricing model. All used data come from the CSMAR Database, which includes actual option priceC or P, trading date t, exercise date T, strike price K, and current price of underlying asset S0. In addition, the database also includes the historical volatility of Shanghai 50 ETF σ=0.1994 on September 3, 2018, and σ=0.2005 on September 4, 2018, and the 1-year deposit benchmark interest rate of Chinese Central Bank r=1.5%, which is chosen as the reference level of risk-free interest rate in the Chinese financial market. All data analyses in the following are worked out by the software MATLAB R2018b.


### 4.1. Parameter Estimation

There are 48 call options on September 3, 2018, and the current price of underlying assetS0=2.512. We take the annual average return rate of the last month as the drift coefficient of the underlying asset μ. According to simple computation, we obtain μ=2.557% and the left expiration time τ=T−t/365 (years), where t=2018/09/03, and then we work out d1, d2, and d3 by Theorem 2; see Table 1 for detailed data, where C is the actual closing price of call option, T is the expiration time of call option, K is the strike price of call option, S0 is the initial price of underlying asset, τ is the left expiration time of call option, and d1,d2, and d3 are the parameters in Theorem 2. Furthermore, we obtain the estimated value of the risk-compensation factor λ=0.0077963 by the least square method with R2=0.7058.Table 1
Parameter estimation ofd1, d2, and d3 for call options.
CTKS0τd1d2d30.00302018/09/262.752.5120.0630−1.8013−1.7512−1.70120.00732018/09/262.702.5120.0630−1.4347−1.3847−1.33460.01462018/09/262.652.5120.0630−1.0613−1.0112−0.96120.01502018/10/242.752.5120.1397−1.2038−1.1293−1.05470.01592018/12/262.952.5120.3123−1.4263−1.3149−1.20350.02202018/12/262.902.5120.3123−1.2729−1.1615−1.05010.02282018/10/242.702.5120.1397−0.9576−0.8831−0.80860.02752018/09/262.602.5120.0630−0.6807−0.6307−0.58060.02872018/12/262.852.5120.3123−1.1169−1.0054−0.89400.03502018/10/242.652.5120.1397−0.7068−0.6323−0.55780.03722018/12/262.802.5120.3123−0.9580−0.8466−0.73520.04512018/09/262.552.5120.0630−0.2928−0.2427−0.19270.04842018/12/262.752.5120.3123−0.7964−0.6849−0.57350.05002018/10/242.602.5120.1397−0.4513−0.3767−0.30220.06012018/12/262.702.5120.3123−0.6317−0.5203−0.40880.06952018/09/262.502.5120.06300.10280.15290.20290.06982018/10/242.552.5120.1397−0.1908−0.1162−0.04170.07422019/03/272.802.5120.5616−0.7049−0.5555−0.40610.07562018/12/262.652.5120.3123−0.4640−0.3525−0.24110.08782019/03/272.752.5120.5616−0.5844−0.4349−0.28550.09382018/12/262.602.5120.3123−0.2930−0.1816−0.07010.09442018/10/242.502.5120.13970.07490.14950.22400.09802018/09/262.452.5120.06300.50640.55650.60660.10262019/03/272.702.5120.5616−0.4616−0.3121−0.16270.11492018/12/262.552.5120.3123−0.1188−0.00730.10410.11902019/03/272.652.5120.5616−0.3365−0.1870−0.03760.12472018/10/242.452.5120.13970.34600.42050.49500.13612018/09/262.402.5120.06300.91840.96841.01850.13802019/03/272.602.5120.5616−0.2090−0.05960.08990.13902018/12/262.502.5120.31230.05890.17040.28180.15922018/10/242.402.5120.13970.62260.69710.77170.16232019/03/272.552.5120.5616−0.07910.07040.21980.16712018/12/262.452.5120.31230.24020.35170.46310.17702018/09/262.352.5120.06301.33901.38911.43910.18612019/03/272.502.5120.56160.05340.20290.35230.19682018/10/242.352.5120.13970.90510.97961.05410.19812018/12/262.402.5120.31230.42530.53670.64810.21102019/03/272.452.5120.56160.18860.33810.48750.22092018/09/262.302.5120.06301.76871.81871.86880.22972018/12/262.352.5120.31230.61420.72560.83710.23602018/10/242.302.5120.13971.19361.26811.34270.24212019/03/272.402.5120.56160.32660.47610.62550.26602018/12/262.302.5120.31230.80720.91861.03000.27042018/09/262.252.5120.06302.20782.25782.30790.27672019/03/272.352.5120.56160.46750.61690.76640.30742018/12/262.252.5120.31231.00441.11581.22730.31632018/09/262.202.5120.06302.65672.70682.75680.34502019/03/272.252.5120.56160.75850.90791.0574There are 60 put options on September 3, 2018, and the current price of underlying assetS0=2.512. Similarly, we obtain μ=2.93% and the left expiration time τ=T−t/365 (years), where t=2018/09/03, and then we work out d1, d2, and d3 by Theorem 3; see Table 2 for detailed data, where P is the actual closing price of put option, T is the expiration time of put option, K is the strike price of put option, S0 is the initial price of underlying asset, τ is the left expiration time of put option, and d1,d2, and d3 are the parameters in Theorem 3. Furthermore, we obtain the estimated value of the risk-compensation factor λ=0.013306 by the least square method with R2=0.6218.Table 2
Parameter estimation ofd1, d2, and d3 for put options.
PTKS0τd1d2d30.00152018/09/262.202.5120.06302.65672.70682.75680.00272018/09/262.252.5120.06302.20782.25782.30790.00552018/09/262.302.5120.06301.76871.81871.86880.01152018/09/262.352.5120.06301.33901.38911.43910.01652018/10/242.302.5120.13971.19361.26811.34270.02052018/09/262.402.5120.06300.91840.96841.01850.02472018/12/262.202.5120.31231.20611.31751.42890.02552018/10/242.352.5120.13970.90510.97961.05410.03242018/12/262.252.5120.31231.00441.11581.22730.03432018/09/262.452.5120.06300.50640.55650.60660.03712018/10/242.402.5120.13970.62260.69710.77170.04322018/12/262.302.5120.31230.80720.91861.03000.05282018/10/242.452.5120.13970.34600.42050.49500.05312018/09/262.502.5120.06300.10280.15290.20290.05652018/12/262.352.5120.31230.61420.72560.83710.06082019/03/272.252.5120.56160.75850.90791.05740.07302019/03/272.302.5120.56160.61140.76090.91030.07332018/10/242.502.5120.13970.07490.14950.22400.07342018/12/262.402.5120.31230.42530.53670.64810.07912018/09/262.552.5120.0630−0.2928−0.2427−0.19270.09022019/03/272.352.5120.56160.46750.61690.76640.09202018/12/262.452.5120.31230.24020.35170.46310.09892018/10/242.552.5120.1397−0.1908−0.1162−0.04170.10742019/03/272.402.5120.56160.32660.47610.62550.11172018/09/262.602.5120.0630−0.6807−0.6307−0.58060.11402018/12/262.502.5120.31230.05890.17040.28180.12442019/03/272.452.5120.56160.18860.33810.48750.12842018/10/242.602.5120.1397−0.4513−0.3767−0.30220.14012018/12/262.552.5120.3123−0.1188−0.00730.10410.14832019/03/272.502.5120.56160.05340.20290.35230.14982018/09/262.652.5120.0630−1.0613−1.0112−0.96120.16862018/12/262.602.5120.3123−0.2930−0.1816−0.07010.17502019/03/272.552.5120.5616−0.07910.07040.21980.19082018/09/262.702.5120.0630−1.4347−1.3847−1.33460.20032018/12/262.652.5120.3123−0.4640−0.3525−0.24110.20072018/10/242.702.5120.1397−0.9576−0.8831−0.80860.20092019/03/272.602.5120.5616−0.2090−0.05960.08990.22782018/12/262.702.5120.3123−0.6317−0.5203−0.40880.23122019/03/272.652.5120.5616−0.3365−0.1870−0.03760.23702018/09/262.752.5120.0630−1.8013−1.7512−1.70120.24402018/10/242.752.5120.1397−1.2038−1.1293−1.05470.26352019/03/272.702.5120.5616−0.4616−0.3121−0.16270.26572018/12/262.752.5120.3123−0.7964−0.6849−0.57350.28282018/09/262.802.5120.0630−2.1613−2.1112−2.06120.29762019/03/272.752.5120.5616−0.5844−0.4349−0.28550.30892018/12/262.802.5120.3123−0.9580−0.8466−0.73520.33232018/09/262.852.5120.0630−2.5149−2.4648−2.41480.33302019/03/272.802.5120.5616−0.7049−0.5555−0.40610.34772018/12/262.852.5120.3123−1.1169−1.0054−0.89400.38102018/09/262.902.5120.0630−2.8623−2.8123−2.76220.39392018/12/262.902.5120.3123−1.2729−1.1615−1.05010.43152018/09/262.952.5120.0630−3.2039−3.1538−3.10370.44052018/12/262.952.5120.3123−1.4263−1.3149−1.20350.48102018/09/263.002.5120.0630−3.5396−3.4896−3.43950.58132018/09/263.102.5120.0630−4.1947−4.1447−4.09460.68092018/09/263.202.5120.0630−4.8290−4.7789−4.72890.77982018/09/263.302.5120.0630−5.4438−5.3937−5.34370.87932018/09/263.402.5120.0630−6.0402−5.9901−5.94010.97902018/09/263.502.5120.0630−6.6193−6.5692−6.51921.08072018/09/263.602.5120.0630−7.1821−7.1320−7.0820In addition, we obtain the estimation of risk-free rater=1.768% by minimizing the mean square error of the Black–Scholes model on September 3, 2018. Note that r=1.768% is a little higher than the 1-year deposit benchmark interest rate of Chinese Central Bank 1.5%, so it is very reasonable that we take r=1.768% as the risk-free rate. Furthermore, it is far advantageous for the Black–Scholes model to improve its fitting effect.


### 4.2. Comparison of Pricing Effect

In the section, we will compare the fitting effect of our pricing formulas with the Black–Scholes model.


#### 4.2.1. Comparison of Pricing Effect on Call Options

There are 62 call options on September 4, 2018, and the current price of underlying assetS0=2.55. Based on the estimated parameters μ=2.557%, σ=0.2005, r=1.768%, λ=0.0077963 for call options, the call options of Shanghai 50 ETF on September 4, 2018, are priced by our obtained pricing formula in Theorem 3 and the Black–Scholes model, respectively; see Table 3 for detailed data, where T is the expiration time of call option, K is the strike price of call option, S0 is the initial price of underlying asset, τ is the left expiration time of call option, d1,d2, and d3 are the parameters in Theorem 3, C1 is the value of call option computed by the Black–Scholes model, C2 is the value of call option computed by Theorem 3, and C is the actual closing price of call option.Table 3
Pricing results by our pricing formula and the Black–Scholes model for call options.
TKS0τd1d2d3C1C2C2018/09/263.002.550.0603−3.2949−3.2457−3.19650.00050.00000.00052018/12/262.402.550.30960.55860.67020.78170.22390.20870.22002018/12/262.352.550.30960.74730.85890.97050.26010.24510.25552018/09/262.902.550.0603−2.6062−2.5570−2.50770.00170.00020.00072019/03/272.502.550.55890.15250.30240.45230.20790.18830.20362019/03/272.552.550.55890.02040.17030.32020.18160.16260.17732019/03/272.652.550.5589−0.2362−0.08630.06360.13640.11880.13532018/12/262.602.550.3096−0.1589−0.04730.06430.11070.09630.10672018/09/262.202.550.06033.00603.05523.10440.35480.35200.35102018/09/262.852.550.0603−2.2529−2.2036−2.15440.00300.00060.00122019/03/272.202.550.55891.00531.15521.30510.41600.39550.40702018/12/262.302.550.30960.94011.05171.16320.29910.28460.29252018/10/242.352.550.13701.11081.18501.25920.22640.21660.22552018/10/242.302.550.13701.40061.47481.54900.26980.26100.26532018/10/242.702.550.1370−0.7601−0.6859−0.61170.03560.02580.03002018/12/262.202.550.30961.33861.45011.56170.38400.37060.38172018/09/262.502.550.06030.40900.45820.50740.09050.08000.08792018/09/262.952.550.0603−2.9535−2.9042−2.85500.00100.00010.00052018/12/262.952.550.3096−1.2909−1.1794−1.06780.02290.01500.01972018/09/262.802.550.0603−1.8933−1.8441−1.79480.00520.00150.00192018/12/262.852.550.3096−0.9818−0.8703−0.75870.03710.02720.03412018/10/242.552.550.13700.01010.08430.15850.09010.07800.08542018/10/242.402.550.13700.82710.90130.97550.18620.17530.18372018/09/263.102.550.0603−3.9610−3.9118−3.86260.00010.00000.00042018/12/262.702.550.3096−0.4972−0.3856−0.27400.07330.06040.07002019/03/272.352.550.55890.56530.71520.86510.30160.28110.29852019/03/272.402.550.55890.42490.57480.72460.26790.24760.26282018/12/262.752.550.3096−0.6616−0.5501−0.43850.05880.04690.05502019/03/272.252.550.55890.85541.00531.15520.37580.35520.36792019/03/272.802.550.5589−0.6035−0.4536−0.30380.08560.07070.08102018/12/262.652.550.3096−0.3296−0.2181−0.10650.09050.07680.08712018/09/262.302.550.06032.10292.15212.20140.25710.25280.25382018/09/263.602.550.0603−6.9988−6.9496−6.90030.00000.00000.00042019/03/272.752.550.5589−0.4833−0.3334−0.18350.10040.08460.09882019/03/272.302.550.55890.70880.85871.00860.33760.31700.33162018/09/262.402.550.06031.23831.28751.33680.16580.15830.16412018/09/262.252.550.06032.54942.59862.64790.30550.30220.30402019/03/272.602.550.5589−0.10910.04080.19070.15780.13940.15492018/10/242.452.550.13700.54920.62340.69760.14970.13800.14612018/12/262.902.550.3096−1.1377−1.0262−0.91460.02920.02030.02662018/10/242.502.550.13700.27700.35120.42540.11750.10540.11352018/10/242.752.550.1370−1.0074−0.9332−0.85900.02520.01670.01952018/09/263.302.550.0603−5.2311−5.1819−5.13270.00000.00000.00032018/12/262.502.550.30960.19270.30430.41580.16070.14550.15722018/12/262.452.550.30960.37380.48540.59690.19070.17540.18702018/10/242.602.550.1370−0.2516−0.1774−0.10310.06750.05590.06232018/09/262.702.550.0603−1.1545−1.1053−1.05600.01530.00810.01032018/12/262.802.550.3096−0.8232−0.7116−0.60000.04690.03590.04452018/09/262.552.550.06030.00670.05590.10520.06210.05120.05902018/09/263.502.550.0603−6.4265−6.3773−6.32800.00000.00000.00032018/09/262.652.550.0603−0.7747−0.7255−0.67630.02530.01640.02062018/09/262.352.550.06031.66601.71521.76450.21010.20450.20772018/09/262.452.550.06030.81940.86860.91790.12540.11620.12322018/12/262.252.550.30961.13711.24871.36020.34050.32650.33302018/10/242.652.550.1370−0.5083−0.4340−0.35980.04950.03870.04422018/09/262.752.550.0603−1.5272−1.4780−1.42880.00900.00370.00482018/09/263.402.550.0603−5.8376−5.7884−5.73920.00000.00000.00032019/03/272.702.550.5589−0.3609−0.2110−0.06110.11730.10060.11622018/09/263.202.550.0603−4.6060−4.5568−4.50760.00000.00000.00032018/12/262.552.550.30960.01520.12680.23830.13410.11920.13112018/09/262.602.550.0603−0.3878−0.3386−0.28930.04050.03030.03772019/03/272.452.550.55890.28730.43720.58710.23670.21660.2319According to simple computing, the expectation and variance of absolute errors follow as(20)EC1−C^=0.0065574,VarC1−C^=0.0045044,EC2−C^=0.0032696,VarC2−C^=0.001929,where C is the actual closing price of call option, C1 is the value of call option computed by the Black–Scholes model, and C2 is the value of call option computed by Theorem 3. It is obvious that EC2−C^<EC1−C^. In the following, we will support the statement by the hypothesis test (i.e., t-test function of MATLAB R2018b).In fact, for the hypothesis test,(21)H0:EC2−C−EC1−C≥0,H1:EC2−C−EC1−C<0,Thet−statistics equals −5.9311 with degree of freedom 61, and its p value is 7.6261×10−8. Thus, we accept H1, i.e., EC2−C<EC1−C. That is, the prices of call options computed by Theorem 3 are far nearer to their actual prices than those computed by the Black–Scholes model.


#### 4.2.2. Comparison of Pricing Effect on Put Options

There are 62 put options on September 4, 2018, and the current price of underlying assetS0=2.55. Based on the estimated parameters μ=2.557%, σ=0.2005, r=1.768%, and λ=0.013306 for put options, the put options of Shanghai 50 ETF on September 4, 2018, are priced by our pricing formula in Theorem 2 and Black–Scholes model, respectively; see Table 4 for detailed data, where T is the expiration time of put option, K is the strike price of put option, S0 is the initial price of underlying asset, τ is the left expiration time of put option, d1,d2, and d3 are the parameters in Theorem 3, P1 is the value of put option computed by the Black–Scholes model, P2 is the value of put option computed by Theorem 3, and P is the actual closing price of put option.Table 4
Pricing results by our pricing formula and the Black–Scholes model for put options.
TKS0τd1d2d3P1P2P2018/12/262.502.550.30960.19270.30430.41580.09630.08390.09832019/03/272.302.550.55890.70880.85871.00860.05830.04780.06352019/03/272.252.550.55890.85541.00531.15520.04630.03640.05202018/12/262.202.550.30961.33861.45011.56170.01860.01040.02052018/12/262.902.550.3096−1.1377−1.0262−0.91460.36310.35680.36182018/09/262.302.550.06032.10292.15212.20140.00500.00070.00302018/09/263.202.550.0603−4.6060−4.5568−4.50760.64730.64710.64372018/10/242.702.550.1370v0.7601−0.6859−0.61170.18310.17030.17242018/09/262.502.550.06030.40900.45820.50740.04400.02780.03592018/09/262.652.550.0603−0.7747−0.7255−0.67630.12820.11400.11912019/03/272.552.550.55890.02040.17030.32020.15210.14130.15692018/10/242.602.550.1370−0.2516−0.1774−0.10310.11620.10050.10542018/09/262.202.550.06033.00603.05523.10440.00120.00000.00102019/03/272.402.550.55890.42490.57480.72460.08870.07750.09392019/03/272.602.550.5589−0.10910.04080.19070.17800.16770.18292018/09/262.852.550.0603−2.2529−2.2036−2.15440.30100.29800.29552018/12/262.402.550.30960.55860.67020.78170.06050.04750.06072019/03/272.452.550.55890.28730.43720.58710.10740.09620.11212019/03/272.652.550.5589−0.2362−0.08630.06360.20620.19670.20942019/03/272.802.550.5589−0.6035−0.4536−0.30380.30400.29730.30782018/09/262.452.550.06030.81940.86860.91790.02770.01400.02192018/12/262.952.550.3096−1.2909−1.1794−1.06780.40610.40130.40562018/12/262.252.550.30961.13711.24871.36020.02570.01610.02712018/09/262.802.550.0603−1.8933−1.8441−1.79480.25410.24900.24862018/10/242.652.550.1370−0.5083−0.4340−0.35980.14780.13320.13392018/12/262.752.550.3096−0.6616−0.5501−0.43850.24460.23410.23872018/12/262.852.550.3096−0.9818−0.8703−0.75870.32170.31400.32012018/09/262.702.550.0603−1.1545−1.1053−1.05600.16680.15570.15852018/09/263.502.550.0603−6.4265−6.3773−6.32800.94700.94680.94432019/03/272.352.550.55890.56530.71520.86510.07240.06140.07842018/10/242.552.550.13700.01010.08430.15850.08890.07280.07882018/12/262.702.550.3096−0.4972−0.3856−0.27400.20960.19790.20792018/10/242.752.550.1370−1.0074−0.9332−0.85900.22190.21100.21202018/12/262.452.550.30960.37380.48540.59690.07760.06400.07822018/10/242.402.550.13700.82710.90130.97550.03330.02040.02602018/09/263.102.550.0603−3.9610−3.9118−3.86260.54740.54720.54512018/10/242.502.550.13700.27700.35120.42540.06610.05030.05692018/09/263.302.550.0603−5.2311−5.1819−5.13270.74720.74700.74202018/10/242.302.550.13701.40061.47481.54900.01490.00630.01082018/12/262.352.550.30960.74730.85890.97050.04630.03430.04672018/09/262.402.550.06031.23831.28751.33680.01660.00620.01182018/09/263.402.550.0603−5.8376−5.7884−5.73920.84710.84690.84532018/09/262.602.550.0603−0.3878−0.3386−0.28930.09430.07790.08442018/09/262.252.550.06032.54942.59862.64790.00250.00020.00152018/10/242.352.550.13701.11081.18501.25920.02260.01180.01782019/03/272.702.550.5589−0.3609−0.2110−0.06110.23680.22800.24022018/09/262.352.550.06031.66601.71521.76450.00940.00230.00652018/09/262.952.550.0603−2.9535−2.9042−2.85500.39820.39740.39212018/10/242.452.550.13700.54920.62340.69760.04760.03300.03872018/12/262.652.550.3096−0.3296−0.2181−0.10650.17720.16450.17532018/12/262.552.550.30960.01520.12680.23830.12120.10730.12102019/03/272.752.550.5589−0.4833−0.3334−0.18350.26940.26160.27382018/09/262.752.550.0603−1.5272−1.4780−1.42880.20910.20120.20112018/09/263.602.550.0603−6.9988−6.9496−6.90031.04691.04671.04522019/03/272.202.550.55891.00531.15521.30510.03630.02720.04232018/12/262.802.550.3096−0.8232−0.7116−0.60000.28210.27290.28002018/12/262.302.550.30960.94011.05171.16320.03480.02390.03532018/09/262.902.550.0603−2.6062−2.5570−2.50770.34920.34760.34302018/09/262.552.550.06030.00670.05590.10520.06610.04890.05732019/03/272.502.550.55890.15250.30240.45230.12850.11740.13302018/12/262.602.550.3096−0.1589−0.04730.06430.14770.13420.14682018/09/263.002.550.0603−3.2949−3.2457−3.19650.44770.44730.4459According to simple computing, the expectation and variance of absolute errors follow as(22)EP1−P^=0.00836,VarP1−P^=0.0056946,EP2−P^=0.0044808,VarP2−P^=0.0029812,where P is the actual closing price of put option, P1 is the value of put option computed by the Black–Scholes model, and P2 is the value of put option computed by Theorem 3. It is obvious that EP2−P^<EP1−P^. In the following, we will support the statement by the hypothesis test (i.e., t-test function of MATLAB R2018b).In fact, for the hypothesis test,(23)H0:EP2−P−EP1−P≥0,H1:EP2−P−EP1−P<0,Thet−statistics equals −4.7567 with degree of freedom 61, and its p value is 6.2304×10−6. Thus, we accept H1, i.e., EP2−P<EP1−P. That is, the prices of put options computed by Theorem 3 are far nearer to their actual prices than those computed by the Black–Scholes model.Therefore, our pricing formulas in Theorem2 and Theorem 3 have less absolute errors than those of the Black–Scholes model for both call options and put options. That is, the fitting effect of our pricing formulas is superior to that of the Black–Scholes model.


## 4.1. Parameter Estimation

There are 48 call options on September 3, 2018, and the current price of underlying assetS0=2.512. We take the annual average return rate of the last month as the drift coefficient of the underlying asset μ. According to simple computation, we obtain μ=2.557% and the left expiration time τ=T−t/365 (years), where t=2018/09/03, and then we work out d1, d2, and d3 by Theorem 2; see Table 1 for detailed data, where C is the actual closing price of call option, T is the expiration time of call option, K is the strike price of call option, S0 is the initial price of underlying asset, τ is the left expiration time of call option, and d1,d2, and d3 are the parameters in Theorem 2. Furthermore, we obtain the estimated value of the risk-compensation factor λ=0.0077963 by the least square method with R2=0.7058.Table 1
Parameter estimation ofd1, d2, and d3 for call options.
CTKS0τd1d2d30.00302018/09/262.752.5120.0630−1.8013−1.7512−1.70120.00732018/09/262.702.5120.0630−1.4347−1.3847−1.33460.01462018/09/262.652.5120.0630−1.0613−1.0112−0.96120.01502018/10/242.752.5120.1397−1.2038−1.1293−1.05470.01592018/12/262.952.5120.3123−1.4263−1.3149−1.20350.02202018/12/262.902.5120.3123−1.2729−1.1615−1.05010.02282018/10/242.702.5120.1397−0.9576−0.8831−0.80860.02752018/09/262.602.5120.0630−0.6807−0.6307−0.58060.02872018/12/262.852.5120.3123−1.1169−1.0054−0.89400.03502018/10/242.652.5120.1397−0.7068−0.6323−0.55780.03722018/12/262.802.5120.3123−0.9580−0.8466−0.73520.04512018/09/262.552.5120.0630−0.2928−0.2427−0.19270.04842018/12/262.752.5120.3123−0.7964−0.6849−0.57350.05002018/10/242.602.5120.1397−0.4513−0.3767−0.30220.06012018/12/262.702.5120.3123−0.6317−0.5203−0.40880.06952018/09/262.502.5120.06300.10280.15290.20290.06982018/10/242.552.5120.1397−0.1908−0.1162−0.04170.07422019/03/272.802.5120.5616−0.7049−0.5555−0.40610.07562018/12/262.652.5120.3123−0.4640−0.3525−0.24110.08782019/03/272.752.5120.5616−0.5844−0.4349−0.28550.09382018/12/262.602.5120.3123−0.2930−0.1816−0.07010.09442018/10/242.502.5120.13970.07490.14950.22400.09802018/09/262.452.5120.06300.50640.55650.60660.10262019/03/272.702.5120.5616−0.4616−0.3121−0.16270.11492018/12/262.552.5120.3123−0.1188−0.00730.10410.11902019/03/272.652.5120.5616−0.3365−0.1870−0.03760.12472018/10/242.452.5120.13970.34600.42050.49500.13612018/09/262.402.5120.06300.91840.96841.01850.13802019/03/272.602.5120.5616−0.2090−0.05960.08990.13902018/12/262.502.5120.31230.05890.17040.28180.15922018/10/242.402.5120.13970.62260.69710.77170.16232019/03/272.552.5120.5616−0.07910.07040.21980.16712018/12/262.452.5120.31230.24020.35170.46310.17702018/09/262.352.5120.06301.33901.38911.43910.18612019/03/272.502.5120.56160.05340.20290.35230.19682018/10/242.352.5120.13970.90510.97961.05410.19812018/12/262.402.5120.31230.42530.53670.64810.21102019/03/272.452.5120.56160.18860.33810.48750.22092018/09/262.302.5120.06301.76871.81871.86880.22972018/12/262.352.5120.31230.61420.72560.83710.23602018/10/242.302.5120.13971.19361.26811.34270.24212019/03/272.402.5120.56160.32660.47610.62550.26602018/12/262.302.5120.31230.80720.91861.03000.27042018/09/262.252.5120.06302.20782.25782.30790.27672019/03/272.352.5120.56160.46750.61690.76640.30742018/12/262.252.5120.31231.00441.11581.22730.31632018/09/262.202.5120.06302.65672.70682.75680.34502019/03/272.252.5120.56160.75850.90791.0574There are 60 put options on September 3, 2018, and the current price of underlying assetS0=2.512. Similarly, we obtain μ=2.93% and the left expiration time τ=T−t/365 (years), where t=2018/09/03, and then we work out d1, d2, and d3 by Theorem 3; see Table 2 for detailed data, where P is the actual closing price of put option, T is the expiration time of put option, K is the strike price of put option, S0 is the initial price of underlying asset, τ is the left expiration time of put option, and d1,d2, and d3 are the parameters in Theorem 3. Furthermore, we obtain the estimated value of the risk-compensation factor λ=0.013306 by the least square method with R2=0.6218.Table 2
Parameter estimation ofd1, d2, and d3 for put options.
PTKS0τd1d2d30.00152018/09/262.202.5120.06302.65672.70682.75680.00272018/09/262.252.5120.06302.20782.25782.30790.00552018/09/262.302.5120.06301.76871.81871.86880.01152018/09/262.352.5120.06301.33901.38911.43910.01652018/10/242.302.5120.13971.19361.26811.34270.02052018/09/262.402.5120.06300.91840.96841.01850.02472018/12/262.202.5120.31231.20611.31751.42890.02552018/10/242.352.5120.13970.90510.97961.05410.03242018/12/262.252.5120.31231.00441.11581.22730.03432018/09/262.452.5120.06300.50640.55650.60660.03712018/10/242.402.5120.13970.62260.69710.77170.04322018/12/262.302.5120.31230.80720.91861.03000.05282018/10/242.452.5120.13970.34600.42050.49500.05312018/09/262.502.5120.06300.10280.15290.20290.05652018/12/262.352.5120.31230.61420.72560.83710.06082019/03/272.252.5120.56160.75850.90791.05740.07302019/03/272.302.5120.56160.61140.76090.91030.07332018/10/242.502.5120.13970.07490.14950.22400.07342018/12/262.402.5120.31230.42530.53670.64810.07912018/09/262.552.5120.0630−0.2928−0.2427−0.19270.09022019/03/272.352.5120.56160.46750.61690.76640.09202018/12/262.452.5120.31230.24020.35170.46310.09892018/10/242.552.5120.1397−0.1908−0.1162−0.04170.10742019/03/272.402.5120.56160.32660.47610.62550.11172018/09/262.602.5120.0630−0.6807−0.6307−0.58060.11402018/12/262.502.5120.31230.05890.17040.28180.12442019/03/272.452.5120.56160.18860.33810.48750.12842018/10/242.602.5120.1397−0.4513−0.3767−0.30220.14012018/12/262.552.5120.3123−0.1188−0.00730.10410.14832019/03/272.502.5120.56160.05340.20290.35230.14982018/09/262.652.5120.0630−1.0613−1.0112−0.96120.16862018/12/262.602.5120.3123−0.2930−0.1816−0.07010.17502019/03/272.552.5120.5616−0.07910.07040.21980.19082018/09/262.702.5120.0630−1.4347−1.3847−1.33460.20032018/12/262.652.5120.3123−0.4640−0.3525−0.24110.20072018/10/242.702.5120.1397−0.9576−0.8831−0.80860.20092019/03/272.602.5120.5616−0.2090−0.05960.08990.22782018/12/262.702.5120.3123−0.6317−0.5203−0.40880.23122019/03/272.652.5120.5616−0.3365−0.1870−0.03760.23702018/09/262.752.5120.0630−1.8013−1.7512−1.70120.24402018/10/242.752.5120.1397−1.2038−1.1293−1.05470.26352019/03/272.702.5120.5616−0.4616−0.3121−0.16270.26572018/12/262.752.5120.3123−0.7964−0.6849−0.57350.28282018/09/262.802.5120.0630−2.1613−2.1112−2.06120.29762019/03/272.752.5120.5616−0.5844−0.4349−0.28550.30892018/12/262.802.5120.3123−0.9580−0.8466−0.73520.33232018/09/262.852.5120.0630−2.5149−2.4648−2.41480.33302019/03/272.802.5120.5616−0.7049−0.5555−0.40610.34772018/12/262.852.5120.3123−1.1169−1.0054−0.89400.38102018/09/262.902.5120.0630−2.8623−2.8123−2.76220.39392018/12/262.902.5120.3123−1.2729−1.1615−1.05010.43152018/09/262.952.5120.0630−3.2039−3.1538−3.10370.44052018/12/262.952.5120.3123−1.4263−1.3149−1.20350.48102018/09/263.002.5120.0630−3.5396−3.4896−3.43950.58132018/09/263.102.5120.0630−4.1947−4.1447−4.09460.68092018/09/263.202.5120.0630−4.8290−4.7789−4.72890.77982018/09/263.302.5120.0630−5.4438−5.3937−5.34370.87932018/09/263.402.5120.0630−6.0402−5.9901−5.94010.97902018/09/263.502.5120.0630−6.6193−6.5692−6.51921.08072018/09/263.602.5120.0630−7.1821−7.1320−7.0820In addition, we obtain the estimation of risk-free rater=1.768% by minimizing the mean square error of the Black–Scholes model on September 3, 2018. Note that r=1.768% is a little higher than the 1-year deposit benchmark interest rate of Chinese Central Bank 1.5%, so it is very reasonable that we take r=1.768% as the risk-free rate. Furthermore, it is far advantageous for the Black–Scholes model to improve its fitting effect.


## 4.2. Comparison of Pricing Effect

In the section, we will compare the fitting effect of our pricing formulas with the Black–Scholes model.


### 4.2.1. Comparison of Pricing Effect on Call Options

There are 62 call options on September 4, 2018, and the current price of underlying assetS0=2.55. Based on the estimated parameters μ=2.557%, σ=0.2005, r=1.768%, λ=0.0077963 for call options, the call options of Shanghai 50 ETF on September 4, 2018, are priced by our obtained pricing formula in Theorem 3 and the Black–Scholes model, respectively; see Table 3 for detailed data, where T is the expiration time of call option, K is the strike price of call option, S0 is the initial price of underlying asset, τ is the left expiration time of call option, d1,d2, and d3 are the parameters in Theorem 3, C1 is the value of call option computed by the Black–Scholes model, C2 is the value of call option computed by Theorem 3, and C is the actual closing price of call option.Table 3
Pricing results by our pricing formula and the Black–Scholes model for call options.
TKS0τd1d2d3C1C2C2018/09/263.002.550.0603−3.2949−3.2457−3.19650.00050.00000.00052018/12/262.402.550.30960.55860.67020.78170.22390.20870.22002018/12/262.352.550.30960.74730.85890.97050.26010.24510.25552018/09/262.902.550.0603−2.6062−2.5570−2.50770.00170.00020.00072019/03/272.502.550.55890.15250.30240.45230.20790.18830.20362019/03/272.552.550.55890.02040.17030.32020.18160.16260.17732019/03/272.652.550.5589−0.2362−0.08630.06360.13640.11880.13532018/12/262.602.550.3096−0.1589−0.04730.06430.11070.09630.10672018/09/262.202.550.06033.00603.05523.10440.35480.35200.35102018/09/262.852.550.0603−2.2529−2.2036−2.15440.00300.00060.00122019/03/272.202.550.55891.00531.15521.30510.41600.39550.40702018/12/262.302.550.30960.94011.05171.16320.29910.28460.29252018/10/242.352.550.13701.11081.18501.25920.22640.21660.22552018/10/242.302.550.13701.40061.47481.54900.26980.26100.26532018/10/242.702.550.1370−0.7601−0.6859−0.61170.03560.02580.03002018/12/262.202.550.30961.33861.45011.56170.38400.37060.38172018/09/262.502.550.06030.40900.45820.50740.09050.08000.08792018/09/262.952.550.0603−2.9535−2.9042−2.85500.00100.00010.00052018/12/262.952.550.3096−1.2909−1.1794−1.06780.02290.01500.01972018/09/262.802.550.0603−1.8933−1.8441−1.79480.00520.00150.00192018/12/262.852.550.3096−0.9818−0.8703−0.75870.03710.02720.03412018/10/242.552.550.13700.01010.08430.15850.09010.07800.08542018/10/242.402.550.13700.82710.90130.97550.18620.17530.18372018/09/263.102.550.0603−3.9610−3.9118−3.86260.00010.00000.00042018/12/262.702.550.3096−0.4972−0.3856−0.27400.07330.06040.07002019/03/272.352.550.55890.56530.71520.86510.30160.28110.29852019/03/272.402.550.55890.42490.57480.72460.26790.24760.26282018/12/262.752.550.3096−0.6616−0.5501−0.43850.05880.04690.05502019/03/272.252.550.55890.85541.00531.15520.37580.35520.36792019/03/272.802.550.5589−0.6035−0.4536−0.30380.08560.07070.08102018/12/262.652.550.3096−0.3296−0.2181−0.10650.09050.07680.08712018/09/262.302.550.06032.10292.15212.20140.25710.25280.25382018/09/263.602.550.0603−6.9988−6.9496−6.90030.00000.00000.00042019/03/272.752.550.5589−0.4833−0.3334−0.18350.10040.08460.09882019/03/272.302.550.55890.70880.85871.00860.33760.31700.33162018/09/262.402.550.06031.23831.28751.33680.16580.15830.16412018/09/262.252.550.06032.54942.59862.64790.30550.30220.30402019/03/272.602.550.5589−0.10910.04080.19070.15780.13940.15492018/10/242.452.550.13700.54920.62340.69760.14970.13800.14612018/12/262.902.550.3096−1.1377−1.0262−0.91460.02920.02030.02662018/10/242.502.550.13700.27700.35120.42540.11750.10540.11352018/10/242.752.550.1370−1.0074−0.9332−0.85900.02520.01670.01952018/09/263.302.550.0603−5.2311−5.1819−5.13270.00000.00000.00032018/12/262.502.550.30960.19270.30430.41580.16070.14550.15722018/12/262.452.550.30960.37380.48540.59690.19070.17540.18702018/10/242.602.550.1370−0.2516−0.1774−0.10310.06750.05590.06232018/09/262.702.550.0603−1.1545−1.1053−1.05600.01530.00810.01032018/12/262.802.550.3096−0.8232−0.7116−0.60000.04690.03590.04452018/09/262.552.550.06030.00670.05590.10520.06210.05120.05902018/09/263.502.550.0603−6.4265−6.3773−6.32800.00000.00000.00032018/09/262.652.550.0603−0.7747−0.7255−0.67630.02530.01640.02062018/09/262.352.550.06031.66601.71521.76450.21010.20450.20772018/09/262.452.550.06030.81940.86860.91790.12540.11620.12322018/12/262.252.550.30961.13711.24871.36020.34050.32650.33302018/10/242.652.550.1370−0.5083−0.4340−0.35980.04950.03870.04422018/09/262.752.550.0603−1.5272−1.4780−1.42880.00900.00370.00482018/09/263.402.550.0603−5.8376−5.7884−5.73920.00000.00000.00032019/03/272.702.550.5589−0.3609−0.2110−0.06110.11730.10060.11622018/09/263.202.550.0603−4.6060−4.5568−4.50760.00000.00000.00032018/12/262.552.550.30960.01520.12680.23830.13410.11920.13112018/09/262.602.550.0603−0.3878−0.3386−0.28930.04050.03030.03772019/03/272.452.550.55890.28730.43720.58710.23670.21660.2319According to simple computing, the expectation and variance of absolute errors follow as(20)EC1−C^=0.0065574,VarC1−C^=0.0045044,EC2−C^=0.0032696,VarC2−C^=0.001929,where C is the actual closing price of call option, C1 is the value of call option computed by the Black–Scholes model, and C2 is the value of call option computed by Theorem 3. It is obvious that EC2−C^<EC1−C^. In the following, we will support the statement by the hypothesis test (i.e., t-test function of MATLAB R2018b).In fact, for the hypothesis test,(21)H0:EC2−C−EC1−C≥0,H1:EC2−C−EC1−C<0,Thet−statistics equals −5.9311 with degree of freedom 61, and its p value is 7.6261×10−8. Thus, we accept H1, i.e., EC2−C<EC1−C. That is, the prices of call options computed by Theorem 3 are far nearer to their actual prices than those computed by the Black–Scholes model.


### 4.2.2. Comparison of Pricing Effect on Put Options

There are 62 put options on September 4, 2018, and the current price of underlying assetS0=2.55. Based on the estimated parameters μ=2.557%, σ=0.2005, r=1.768%, and λ=0.013306 for put options, the put options of Shanghai 50 ETF on September 4, 2018, are priced by our pricing formula in Theorem 2 and Black–Scholes model, respectively; see Table 4 for detailed data, where T is the expiration time of put option, K is the strike price of put option, S0 is the initial price of underlying asset, τ is the left expiration time of put option, d1,d2, and d3 are the parameters in Theorem 3, P1 is the value of put option computed by the Black–Scholes model, P2 is the value of put option computed by Theorem 3, and P is the actual closing price of put option.Table 4
Pricing results by our pricing formula and the Black–Scholes model for put options.
TKS0τd1d2d3P1P2P2018/12/262.502.550.30960.19270.30430.41580.09630.08390.09832019/03/272.302.550.55890.70880.85871.00860.05830.04780.06352019/03/272.252.550.55890.85541.00531.15520.04630.03640.05202018/12/262.202.550.30961.33861.45011.56170.01860.01040.02052018/12/262.902.550.3096−1.1377−1.0262−0.91460.36310.35680.36182018/09/262.302.550.06032.10292.15212.20140.00500.00070.00302018/09/263.202.550.0603−4.6060−4.5568−4.50760.64730.64710.64372018/10/242.702.550.1370v0.7601−0.6859−0.61170.18310.17030.17242018/09/262.502.550.06030.40900.45820.50740.04400.02780.03592018/09/262.652.550.0603−0.7747−0.7255−0.67630.12820.11400.11912019/03/272.552.550.55890.02040.17030.32020.15210.14130.15692018/10/242.602.550.1370−0.2516−0.1774−0.10310.11620.10050.10542018/09/262.202.550.06033.00603.05523.10440.00120.00000.00102019/03/272.402.550.55890.42490.57480.72460.08870.07750.09392019/03/272.602.550.5589−0.10910.04080.19070.17800.16770.18292018/09/262.852.550.0603−2.2529−2.2036−2.15440.30100.29800.29552018/12/262.402.550.30960.55860.67020.78170.06050.04750.06072019/03/272.452.550.55890.28730.43720.58710.10740.09620.11212019/03/272.652.550.5589−0.2362−0.08630.06360.20620.19670.20942019/03/272.802.550.5589−0.6035−0.4536−0.30380.30400.29730.30782018/09/262.452.550.06030.81940.86860.91790.02770.01400.02192018/12/262.952.550.3096−1.2909−1.1794−1.06780.40610.40130.40562018/12/262.252.550.30961.13711.24871.36020.02570.01610.02712018/09/262.802.550.0603−1.8933−1.8441−1.79480.25410.24900.24862018/10/242.652.550.1370−0.5083−0.4340−0.35980.14780.13320.13392018/12/262.752.550.3096−0.6616−0.5501−0.43850.24460.23410.23872018/12/262.852.550.3096−0.9818−0.8703−0.75870.32170.31400.32012018/09/262.702.550.0603−1.1545−1.1053−1.05600.16680.15570.15852018/09/263.502.550.0603−6.4265−6.3773−6.32800.94700.94680.94432019/03/272.352.550.55890.56530.71520.86510.07240.06140.07842018/10/242.552.550.13700.01010.08430.15850.08890.07280.07882018/12/262.702.550.3096−0.4972−0.3856−0.27400.20960.19790.20792018/10/242.752.550.1370−1.0074−0.9332−0.85900.22190.21100.21202018/12/262.452.550.30960.37380.48540.59690.07760.06400.07822018/10/242.402.550.13700.82710.90130.97550.03330.02040.02602018/09/263.102.550.0603−3.9610−3.9118−3.86260.54740.54720.54512018/10/242.502.550.13700.27700.35120.42540.06610.05030.05692018/09/263.302.550.0603−5.2311−5.1819−5.13270.74720.74700.74202018/10/242.302.550.13701.40061.47481.54900.01490.00630.01082018/12/262.352.550.30960.74730.85890.97050.04630.03430.04672018/09/262.402.550.06031.23831.28751.33680.01660.00620.01182018/09/263.402.550.0603−5.8376−5.7884−5.73920.84710.84690.84532018/09/262.602.550.0603−0.3878−0.3386−0.28930.09430.07790.08442018/09/262.252.550.06032.54942.59862.64790.00250.00020.00152018/10/242.352.550.13701.11081.18501.25920.02260.01180.01782019/03/272.702.550.5589−0.3609−0.2110−0.06110.23680.22800.24022018/09/262.352.550.06031.66601.71521.76450.00940.00230.00652018/09/262.952.550.0603−2.9535−2.9042−2.85500.39820.39740.39212018/10/242.452.550.13700.54920.62340.69760.04760.03300.03872018/12/262.652.550.3096−0.3296−0.2181−0.10650.17720.16450.17532018/12/262.552.550.30960.01520.12680.23830.12120.10730.12102019/03/272.752.550.5589−0.4833−0.3334−0.18350.26940.26160.27382018/09/262.752.550.0603−1.5272−1.4780−1.42880.20910.20120.20112018/09/263.602.550.0603−6.9988−6.9496−6.90031.04691.04671.04522019/03/272.202.550.55891.00531.15521.30510.03630.02720.04232018/12/262.802.550.3096−0.8232−0.7116−0.60000.28210.27290.28002018/12/262.302.550.30960.94011.05171.16320.03480.02390.03532018/09/262.902.550.0603−2.6062−2.5570−2.50770.34920.34760.34302018/09/262.552.550.06030.00670.05590.10520.06610.04890.05732019/03/272.502.550.55890.15250.30240.45230.12850.11740.13302018/12/262.602.550.3096−0.1589−0.04730.06430.14770.13420.14682018/09/263.002.550.0603−3.2949−3.2457−3.19650.44770.44730.4459According to simple computing, the expectation and variance of absolute errors follow as(22)EP1−P^=0.00836,VarP1−P^=0.0056946,EP2−P^=0.0044808,VarP2−P^=0.0029812,where P is the actual closing price of put option, P1 is the value of put option computed by the Black–Scholes model, and P2 is the value of put option computed by Theorem 3. It is obvious that EP2−P^<EP1−P^. In the following, we will support the statement by the hypothesis test (i.e., t-test function of MATLAB R2018b).In fact, for the hypothesis test,(23)H0:EP2−P−EP1−P≥0,H1:EP2−P−EP1−P<0,Thet−statistics equals −4.7567 with degree of freedom 61, and its p value is 6.2304×10−6. Thus, we accept H1, i.e., EP2−P<EP1−P. That is, the prices of put options computed by Theorem 3 are far nearer to their actual prices than those computed by the Black–Scholes model.Therefore, our pricing formulas in Theorem2 and Theorem 3 have less absolute errors than those of the Black–Scholes model for both call options and put options. That is, the fitting effect of our pricing formulas is superior to that of the Black–Scholes model.


## 4.2.1. Comparison of Pricing Effect on Call Options

There are 62 call options on September 4, 2018, and the current price of underlying assetS0=2.55. Based on the estimated parameters μ=2.557%, σ=0.2005, r=1.768%, λ=0.0077963 for call options, the call options of Shanghai 50 ETF on September 4, 2018, are priced by our obtained pricing formula in Theorem 3 and the Black–Scholes model, respectively; see Table 3 for detailed data, where T is the expiration time of call option, K is the strike price of call option, S0 is the initial price of underlying asset, τ is the left expiration time of call option, d1,d2, and d3 are the parameters in Theorem 3, C1 is the value of call option computed by the Black–Scholes model, C2 is the value of call option computed by Theorem 3, and C is the actual closing price of call option.Table 3
Pricing results by our pricing formula and the Black–Scholes model for call options.
TKS0τd1d2d3C1C2C2018/09/263.002.550.0603−3.2949−3.2457−3.19650.00050.00000.00052018/12/262.402.550.30960.55860.67020.78170.22390.20870.22002018/12/262.352.550.30960.74730.85890.97050.26010.24510.25552018/09/262.902.550.0603−2.6062−2.5570−2.50770.00170.00020.00072019/03/272.502.550.55890.15250.30240.45230.20790.18830.20362019/03/272.552.550.55890.02040.17030.32020.18160.16260.17732019/03/272.652.550.5589−0.2362−0.08630.06360.13640.11880.13532018/12/262.602.550.3096−0.1589−0.04730.06430.11070.09630.10672018/09/262.202.550.06033.00603.05523.10440.35480.35200.35102018/09/262.852.550.0603−2.2529−2.2036−2.15440.00300.00060.00122019/03/272.202.550.55891.00531.15521.30510.41600.39550.40702018/12/262.302.550.30960.94011.05171.16320.29910.28460.29252018/10/242.352.550.13701.11081.18501.25920.22640.21660.22552018/10/242.302.550.13701.40061.47481.54900.26980.26100.26532018/10/242.702.550.1370−0.7601−0.6859−0.61170.03560.02580.03002018/12/262.202.550.30961.33861.45011.56170.38400.37060.38172018/09/262.502.550.06030.40900.45820.50740.09050.08000.08792018/09/262.952.550.0603−2.9535−2.9042−2.85500.00100.00010.00052018/12/262.952.550.3096−1.2909−1.1794−1.06780.02290.01500.01972018/09/262.802.550.0603−1.8933−1.8441−1.79480.00520.00150.00192018/12/262.852.550.3096−0.9818−0.8703−0.75870.03710.02720.03412018/10/242.552.550.13700.01010.08430.15850.09010.07800.08542018/10/242.402.550.13700.82710.90130.97550.18620.17530.18372018/09/263.102.550.0603−3.9610−3.9118−3.86260.00010.00000.00042018/12/262.702.550.3096−0.4972−0.3856−0.27400.07330.06040.07002019/03/272.352.550.55890.56530.71520.86510.30160.28110.29852019/03/272.402.550.55890.42490.57480.72460.26790.24760.26282018/12/262.752.550.3096−0.6616−0.5501−0.43850.05880.04690.05502019/03/272.252.550.55890.85541.00531.15520.37580.35520.36792019/03/272.802.550.5589−0.6035−0.4536−0.30380.08560.07070.08102018/12/262.652.550.3096−0.3296−0.2181−0.10650.09050.07680.08712018/09/262.302.550.06032.10292.15212.20140.25710.25280.25382018/09/263.602.550.0603−6.9988−6.9496−6.90030.00000.00000.00042019/03/272.752.550.5589−0.4833−0.3334−0.18350.10040.08460.09882019/03/272.302.550.55890.70880.85871.00860.33760.31700.33162018/09/262.402.550.06031.23831.28751.33680.16580.15830.16412018/09/262.252.550.06032.54942.59862.64790.30550.30220.30402019/03/272.602.550.5589−0.10910.04080.19070.15780.13940.15492018/10/242.452.550.13700.54920.62340.69760.14970.13800.14612018/12/262.902.550.3096−1.1377−1.0262−0.91460.02920.02030.02662018/10/242.502.550.13700.27700.35120.42540.11750.10540.11352018/10/242.752.550.1370−1.0074−0.9332−0.85900.02520.01670.01952018/09/263.302.550.0603−5.2311−5.1819−5.13270.00000.00000.00032018/12/262.502.550.30960.19270.30430.41580.16070.14550.15722018/12/262.452.550.30960.37380.48540.59690.19070.17540.18702018/10/242.602.550.1370−0.2516−0.1774−0.10310.06750.05590.06232018/09/262.702.550.0603−1.1545−1.1053−1.05600.01530.00810.01032018/12/262.802.550.3096−0.8232−0.7116−0.60000.04690.03590.04452018/09/262.552.550.06030.00670.05590.10520.06210.05120.05902018/09/263.502.550.0603−6.4265−6.3773−6.32800.00000.00000.00032018/09/262.652.550.0603−0.7747−0.7255−0.67630.02530.01640.02062018/09/262.352.550.06031.66601.71521.76450.21010.20450.20772018/09/262.452.550.06030.81940.86860.91790.12540.11620.12322018/12/262.252.550.30961.13711.24871.36020.34050.32650.33302018/10/242.652.550.1370−0.5083−0.4340−0.35980.04950.03870.04422018/09/262.752.550.0603−1.5272−1.4780−1.42880.00900.00370.00482018/09/263.402.550.0603−5.8376−5.7884−5.73920.00000.00000.00032019/03/272.702.550.5589−0.3609−0.2110−0.06110.11730.10060.11622018/09/263.202.550.0603−4.6060−4.5568−4.50760.00000.00000.00032018/12/262.552.550.30960.01520.12680.23830.13410.11920.13112018/09/262.602.550.0603−0.3878−0.3386−0.28930.04050.03030.03772019/03/272.452.550.55890.28730.43720.58710.23670.21660.2319According to simple computing, the expectation and variance of absolute errors follow as(20)EC1−C^=0.0065574,VarC1−C^=0.0045044,EC2−C^=0.0032696,VarC2−C^=0.001929,where C is the actual closing price of call option, C1 is the value of call option computed by the Black–Scholes model, and C2 is the value of call option computed by Theorem 3. It is obvious that EC2−C^<EC1−C^. In the following, we will support the statement by the hypothesis test (i.e., t-test function of MATLAB R2018b).In fact, for the hypothesis test,(21)H0:EC2−C−EC1−C≥0,H1:EC2−C−EC1−C<0,Thet−statistics equals −5.9311 with degree of freedom 61, and its p value is 7.6261×10−8. Thus, we accept H1, i.e., EC2−C<EC1−C. That is, the prices of call options computed by Theorem 3 are far nearer to their actual prices than those computed by the Black–Scholes model.


## 4.2.2. Comparison of Pricing Effect on Put Options

There are 62 put options on September 4, 2018, and the current price of underlying assetS0=2.55. Based on the estimated parameters μ=2.557%, σ=0.2005, r=1.768%, and λ=0.013306 for put options, the put options of Shanghai 50 ETF on September 4, 2018, are priced by our pricing formula in Theorem 2 and Black–Scholes model, respectively; see Table 4 for detailed data, where T is the expiration time of put option, K is the strike price of put option, S0 is the initial price of underlying asset, τ is the left expiration time of put option, d1,d2, and d3 are the parameters in Theorem 3, P1 is the value of put option computed by the Black–Scholes model, P2 is the value of put option computed by Theorem 3, and P is the actual closing price of put option.Table 4
Pricing results by our pricing formula and the Black–Scholes model for put options.
TKS0τd1d2d3P1P2P2018/12/262.502.550.30960.19270.30430.41580.09630.08390.09832019/03/272.302.550.55890.70880.85871.00860.05830.04780.06352019/03/272.252.550.55890.85541.00531.15520.04630.03640.05202018/12/262.202.550.30961.33861.45011.56170.01860.01040.02052018/12/262.902.550.3096−1.1377−1.0262−0.91460.36310.35680.36182018/09/262.302.550.06032.10292.15212.20140.00500.00070.00302018/09/263.202.550.0603−4.6060−4.5568−4.50760.64730.64710.64372018/10/242.702.550.1370v0.7601−0.6859−0.61170.18310.17030.17242018/09/262.502.550.06030.40900.45820.50740.04400.02780.03592018/09/262.652.550.0603−0.7747−0.7255−0.67630.12820.11400.11912019/03/272.552.550.55890.02040.17030.32020.15210.14130.15692018/10/242.602.550.1370−0.2516−0.1774−0.10310.11620.10050.10542018/09/262.202.550.06033.00603.05523.10440.00120.00000.00102019/03/272.402.550.55890.42490.57480.72460.08870.07750.09392019/03/272.602.550.5589−0.10910.04080.19070.17800.16770.18292018/09/262.852.550.0603−2.2529−2.2036−2.15440.30100.29800.29552018/12/262.402.550.30960.55860.67020.78170.06050.04750.06072019/03/272.452.550.55890.28730.43720.58710.10740.09620.11212019/03/272.652.550.5589−0.2362−0.08630.06360.20620.19670.20942019/03/272.802.550.5589−0.6035−0.4536−0.30380.30400.29730.30782018/09/262.452.550.06030.81940.86860.91790.02770.01400.02192018/12/262.952.550.3096−1.2909−1.1794−1.06780.40610.40130.40562018/12/262.252.550.30961.13711.24871.36020.02570.01610.02712018/09/262.802.550.0603−1.8933−1.8441−1.79480.25410.24900.24862018/10/242.652.550.1370−0.5083−0.4340−0.35980.14780.13320.13392018/12/262.752.550.3096−0.6616−0.5501−0.43850.24460.23410.23872018/12/262.852.550.3096−0.9818−0.8703−0.75870.32170.31400.32012018/09/262.702.550.0603−1.1545−1.1053−1.05600.16680.15570.15852018/09/263.502.550.0603−6.4265−6.3773−6.32800.94700.94680.94432019/03/272.352.550.55890.56530.71520.86510.07240.06140.07842018/10/242.552.550.13700.01010.08430.15850.08890.07280.07882018/12/262.702.550.3096−0.4972−0.3856−0.27400.20960.19790.20792018/10/242.752.550.1370−1.0074−0.9332−0.85900.22190.21100.21202018/12/262.452.550.30960.37380.48540.59690.07760.06400.07822018/10/242.402.550.13700.82710.90130.97550.03330.02040.02602018/09/263.102.550.0603−3.9610−3.9118−3.86260.54740.54720.54512018/10/242.502.550.13700.27700.35120.42540.06610.05030.05692018/09/263.302.550.0603−5.2311−5.1819−5.13270.74720.74700.74202018/10/242.302.550.13701.40061.47481.54900.01490.00630.01082018/12/262.352.550.30960.74730.85890.97050.04630.03430.04672018/09/262.402.550.06031.23831.28751.33680.01660.00620.01182018/09/263.402.550.0603−5.8376−5.7884−5.73920.84710.84690.84532018/09/262.602.550.0603−0.3878−0.3386−0.28930.09430.07790.08442018/09/262.252.550.06032.54942.59862.64790.00250.00020.00152018/10/242.352.550.13701.11081.18501.25920.02260.01180.01782019/03/272.702.550.5589−0.3609−0.2110−0.06110.23680.22800.24022018/09/262.352.550.06031.66601.71521.76450.00940.00230.00652018/09/262.952.550.0603−2.9535−2.9042−2.85500.39820.39740.39212018/10/242.452.550.13700.54920.62340.69760.04760.03300.03872018/12/262.652.550.3096−0.3296−0.2181−0.10650.17720.16450.17532018/12/262.552.550.30960.01520.12680.23830.12120.10730.12102019/03/272.752.550.5589−0.4833−0.3334−0.18350.26940.26160.27382018/09/262.752.550.0603−1.5272−1.4780−1.42880.20910.20120.20112018/09/263.602.550.0603−6.9988−6.9496−6.90031.04691.04671.04522019/03/272.202.550.55891.00531.15521.30510.03630.02720.04232018/12/262.802.550.3096−0.8232−0.7116−0.60000.28210.27290.28002018/12/262.302.550.30960.94011.05171.16320.03480.02390.03532018/09/262.902.550.0603−2.6062−2.5570−2.50770.34920.34760.34302018/09/262.552.550.06030.00670.05590.10520.06610.04890.05732019/03/272.502.550.55890.15250.30240.45230.12850.11740.13302018/12/262.602.550.3096−0.1589−0.04730.06430.14770.13420.14682018/09/263.002.550.0603−3.2949−3.2457−3.19650.44770.44730.4459According to simple computing, the expectation and variance of absolute errors follow as(22)EP1−P^=0.00836,VarP1−P^=0.0056946,EP2−P^=0.0044808,VarP2−P^=0.0029812,where P is the actual closing price of put option, P1 is the value of put option computed by the Black–Scholes model, and P2 is the value of put option computed by Theorem 3. It is obvious that EP2−P^<EP1−P^. In the following, we will support the statement by the hypothesis test (i.e., t-test function of MATLAB R2018b).In fact, for the hypothesis test,(23)H0:EP2−P−EP1−P≥0,H1:EP2−P−EP1−P<0,Thet−statistics equals −4.7567 with degree of freedom 61, and its p value is 6.2304×10−6. Thus, we accept H1, i.e., EP2−P<EP1−P. That is, the prices of put options computed by Theorem 3 are far nearer to their actual prices than those computed by the Black–Scholes model.Therefore, our pricing formulas in Theorem2 and Theorem 3 have less absolute errors than those of the Black–Scholes model for both call options and put options. That is, the fitting effect of our pricing formulas is superior to that of the Black–Scholes model.


## 5. Empirical Analysis of S&P 500 Options

In the section, we will present empirical analysis on S&P 500 options and use the data of April 1 and 2, 2019 to compare the fitting effect of our pricing model and the Black–Scholes pricing model. All used data come from the Chicago Board Options Exchange, which includes actual option priceC or P, trading date t, exercise date T, strike price K, and current price of underlying asset S0. We consider the out-of-money put and call options, which are more liquid and actively traded than in-the-money options. And observations with trading volume below average, prices less than $0.5 or left expiration time less than 10 days or longer than 360 days are discarded. The annualized historical volatility σ=0.1174 based on closing prices of the underlying asset over the last month. According to the Board of Governors of the Federal Reserve System (https://www.federalreserve.gov/releases/h15/data.htm), the annualized risk-free interest rate is 2.43% on April 1 and 2.42% on April 2, 2019.


### 5.1. Parameter Estimation

We consider the daily logarithmic returns of S&P 500 index closing prices from January 2015 to December 2018. The augmented Dickey–Fuller test shows the Dickey–Fuller statistic equals−10.283 with lag order 10 and p value is 0.01, which indicates that the time series of returns is stationary. Furthermore, the autocorrelogram shows the coefficients of autocorrelation mostly fall within double standard deviations (see Figure 1). Thus, we accept that the time series of daily logarithmic returns of S&P 500 index closing prices are stationary and independent. We take the annual average return rate based on 252 effective trading days as the drift coefficient of the underlying asset μ.Figure 1
The autocorrelogram of the daily logarithmic returns of S&P 500 index closing prices.There are 41 call options on April 1, 2019, and the current price of underlying assetS0=2867.19. According to simple computation, we obtain μ=4.932% and the left expiration time τ=T−t/365 (years), where t=2019/04/01, and then we work out d1, d2, and d3 by Theorem 2; see Table 5 for detailed data, where C is the actual closing price of call option, τ is the left expiration time of call option, K is the strike price of call option, and d1,d2, and d3 are the parameters in Theorem 2. Furthermore, we obtain the estimated value of the risk-compensation factor λ=−0.0246533 by the least square method with R2=0.9532.Table 5
Parameter estimation ofd1, d2, and d3 for call options.
CτKd1d2d324.400.04657528700.0393200.0646590.08999721.500.0465752875−0.029375−0.0040370.02130210.600.0465752900−0.371070−0.345732−0.3203934.500.0465752925−0.709832−0.684493−0.6591553.850.0465752930−0.777237−0.751898−0.7265593.140.0465752935−0.844526−0.819188−0.7938491.670.0465752950−1.045710−1.020372−0.9950330.650.0465752975−1.378754−1.353416−1.3280770.600.0465752980−1.445027−1.419689−1.39435046.150.12602728700.1047710.1464520.18813344.000.12602728750.0630100.1046910.14637131.000.1260272900−0.144713−0.103032−0.06135119.350.1260272925−0.350653−0.308972−0.2672917.000.1260272975−0.757303−0.715622−0.6739413.730.1260273000−0.958073−0.916392−0.8747113.300.1260273005−0.998026−0.956345−0.9146642.500.1260273015−1.077733−1.036052−0.99437166.700.22191828700.1525040.2078140.26312363.800.22191828750.1210330.1763430.23165248.600.2219182900−0.0355050.0198050.07511436.700.2219182925−0.190699−0.135390−0.08008028.320.2219182940−0.283181−0.227871−0.17256127.700.2219182950−0.344573−0.289263−0.23395413.030.2219183000−0.648446−0.593137−0.5378276.000.2219183050−0.947297−0.891987−0.8366783.100.2219183100−1.241288−1.185978−1.1306690.850.2219183200−1.815306−1.759996−1.70468628.300.2986302975−0.377841−0.313680−0.24951921.900.2986303000−0.508267−0.444106−0.37994517.000.3753423050−0.637911−0.565980−0.49404842.200.4712333000−0.313762−0.233165−0.15256715.100.4712333100−0.720596−0.639998−0.559401133.500.72054828750.2794200.3790830.47874669.700.7205483000−0.147614−0.0479500.05171333.360.7205483100−0.476619−0.376956−0.27729227.500.7205483125−0.557212−0.457549−0.35788515.700.7205483200−0.795178−0.695515−0.5958527.000.7205483300−1.103934−1.004271−0.90460775.300.7972603000−0.109289−0.0044550.10038018.700.7972603200−0.724912−0.620077−0.51524352.440.9698633100−0.319343−0.203716−0.088089There are 62 put options on April 1, 2019, and the current price of underlying assetS0=2867.19. Similarly, we obtain μ=4.932% and the left expiration time τ=T−t/365 (years), where t=2019/04/01, and then we work out d1, d2, and d3 by Theorem 3; see Table 6 for detailed data, where P is the actual closing price of put option, τ is the left expiration time of put option, K is the strike price of put option, and d1,d2, and d3 are the parameters in Theorem 3. Furthermore, we obtain the estimated value of the risk-compensation factor λ=0.0269767 by the least square method with R2=0.9561.Table 6
Parameter estimation ofd1, d2, and d3 for put options.
PτKd1d2d30.350.04657524007.0974587.1227977.1481360.420.04657524106.9333616.9586996.9840380.400.04657524206.7699426.7952816.8206190.400.04657524356.5260766.5514156.5767530.470.04657524406.4451216.4704606.4957990.460.04657524506.2837086.3090476.3343850.550.04657524556.2032486.2285876.2539260.540.04657524855.7239055.7492435.7745820.620.04657524905.6445775.6699155.6952540.560.04657525005.4863985.5117375.5370760.720.04657525155.2503135.2756525.3009900.740.04657525255.0937045.1190425.1443810.750.04657525504.7048784.7302174.7555550.800.04657526003.9385343.9638723.9892110.930.04657526053.8627113.8880503.9133891.380.04657526403.3359953.3613343.3866731.500.04657526503.1867873.2121263.2374641.800.04657526752.8162172.8415552.8668942.440.04657527002.4490932.4744322.4997702.550.04657527052.3760762.4014152.4267542.650.04657527102.3031952.3285332.3538723.030.04657527202.1578332.1831722.2085113.250.04657527252.0853532.1106922.1360304.500.04657527501.7249351.7502741.7756120.550.12602720008.7697708.8114518.8531321.150.12602722755.6788365.7205175.7621981.500.12602723504.9006554.9423364.9840171.950.12602724004.3955464.4372274.4789072.200.12602724254.1469244.1886054.2302853.200.12602725003.4161523.4578333.4995145.100.12602725752.7069832.7486642.7903456.050.12602726002.4751762.5168572.5585386.380.12602726102.3830772.4247582.4664398.550.12602726502.0181752.0598562.10153710.310.12602726751.7928981.8345791.87626012.350.12602727001.5697171.6113981.65307914.350.12602727251.3485931.3902741.43195518.070.12602727501.1294881.1711691.21285022.540.12602727750.9123660.9540470.99572827.700.12602728000.6971910.7388720.78055332.850.12602728250.4839290.5256100.56729141.250.12602728500.2725470.3142270.3559081.300.22191820006.6823806.7376906.7930005.390.22191823753.5753203.6306303.6859395.800.22191824003.3859993.4413083.4966189.260.22191825002.6479352.7032452.75855411.740.22191825502.2899032.3452122.40052212.900.22191825652.1838612.2391712.29448013.550.22191825752.1135112.1688202.22413015.320.22191826001.9388231.9941332.04944217.700.22191826251.7658071.8211171.87642619.600.22191826501.5944311.6497401.70505026.230.22191827001.2564761.3117861.36709528.700.22191827151.1563101.2116191.26692931.250.22191827251.0898391.1451481.20045835.950.22191827500.9247230.9800331.03534246.700.22191828000.5989480.6542570.70956759.320.22191828500.2789390.3342480.38955863.550.22191828600.2156110.2709200.32623064.700.22191828650.1840300.2393400.29464946.600.29863027500.8478740.9120350.97619657.500.29863028000.5670410.6312030.695364


### 5.2. Comparison of Pricing Effect

In the section, we will compare the fitting effect of our pricing formulas with the Black–Scholes model.


#### 5.2.1. Comparison of Pricing Effect on Call Options

There are 25 call options on April 2, 2019, and the current price of underlying assetS0=2867.24. Based on the estimated parameters μ=4.932%, σ=0.1174, r=2.42%, and λ=−0.0246533 for call options, the call options of S&P 500 on April 2, 2019, are priced by our obtained pricing formula in Theorem 1 and the Black–Scholes model, respectively; see Table 7 for detailed data, where K is the strike price of call option, τ is the left expiration time of call option, d1,d2, and d3 are the parameters in Theorem 3, C is the actual closing price of call option, C1 is the value of call option computed by the Black–Scholes model, and C2 is the value of call option computed by Theorem 1.Table 7
Pricing results by our pricing formula and the Black–Scholes model for call option.
Kτd1d2d3CC1C228700.0438360.0365110.0610930.08567521.0528.2630.7728750.043836−0.034298−0.0097160.01486618.1525.8526.2528800.043836−0.104985−0.080403−0.05582118.2023.5821.9928900.043836−0.245990−0.221408−0.19682613.7019.4514.2729000.043836−0.386509−0.361927−0.3373458.5015.877.6528700.1232880.1035320.1447580.18598345.0050.0654.7728750.1232880.0613100.1025350.14376040.8547.5550.9128950.123288−0.106850−0.065625−0.02439931.6838.3636.7329000.123288−0.148708−0.107483−0.06625829.5036.2733.4929250.123288−0.356924−0.315698−0.27447318.1027.0019.2629400.123288−0.481000−0.439775−0.39854914.1022.3512.2529500.123288−0.563367−0.522141−0.48091611.5019.618.1929550.123288−0.604445−0.563220−0.52199510.4518.346.3528700.2191780.1516570.2066240.26159164.0069.1475.3928750.2191780.1199900.1749570.22992463.8066.5871.7829000.219178−0.0375230.0174440.07241147.5054.7354.9829500.219178−0.348517−0.293550−0.23858326.7535.5827.8130000.219178−0.654284−0.599317−0.54435013.3021.948.9930250.219178−0.805261−0.750294−0.6953278.7516.882.4229000.2958900.0186620.0825280.14639462.7767.3269.1529000.3726030.0620390.1337080.20537676.3878.7481.7831000.372603−0.868516−0.796848−0.7251799.4018.993.5029000.4684930.1059470.1863100.26667391.4091.8796.1030000.468493−0.315908−0.235545−0.15518241.2652.2244.3831000.468493−0.723929−0.643566−0.56320315.0027.1312.45According to simple computing, the expectation and variance of absolute errors follow as(24)EC1−C^=6.902386,VarC1−C^=6.810772,EC2−C^=5.156575,VarC2−C^=9.726556,where C is the actual closing price of call option, C1 is the value of call option computed by the Black–Scholes model, and C2 is the value of call option computed by Theorem 2. It is obvious that EC2−C^<EC1−C^. In the following, we will support the statement by the hypothesis test (i.e., test function of MATLAB R2018b).In fact, for the hypothesis test,(25)H0:EC2−C−EC1−C≥0,H1:EC2−C−EC1−C<0,the t−statistics equals −2.1465 with degree of freedom 24, and its p value is 0.01854. Thus, we accept H1, i.e., EC2−C<EC1−C. That is, the prices of call options computed by Theorem 2 are far nearer to their actual prices than those computed by the Black–Scholes model.


#### 5.2.2. Comparison of Pricing Effect on Put Options

There are 39 put options on April 2, 2019, and the current price of underlying assetS0=2867.24. Based on the estimated parameters μ=4.932%, σ=0.1174, r=2.42%, and λ=0.0269767 for put options, the put options of S&P 500 on April 2, 2019, are priced by our obtained pricing formula in Theorem 3 and the Black–Scholes model, respectively; see Table 8 for detailed data, where K is the strike price of put option, τ is the left expiration time of put option, d1,d2, and d3 are the parameters in Theorem 3, P is the actual closing price of put option, P1 is the value of put option computed by the Black–Scholes model, and P2 is the value of put option computed by Theorem 3.Table 8
Pricing results by our pricing formula and the Black–Scholes model for put option.
Kτd1d2d3PP1P225500.0438364.8456584.8702404.8948220.590.000.0525750.0438364.4487764.4733584.4979400.700.000.1426000.0438364.0557294.0803114.1048930.810.000.3526500.0438363.2808463.3054283.3300101.200.011.5926600.0438363.1276253.1522073.1767901.400.022.0826750.0438362.8988712.9234532.9480351.640.043.0227000.0438362.5204492.5450312.5696131.950.145.2723000.1232885.4740865.5153115.5565361.300.000.0124000.1232884.4417204.4829464.5241711.950.000.1425000.1232883.4515043.4927303.5339553.260.011.1125900.1232882.5936062.6348322.6760575.400.204.6526000.1232882.5001312.5413562.5825816.000.275.3226150.1232882.3605892.4018142.4430406.200.416.4526250.1232882.2680052.3092302.3504567.100.537.2926700.1232881.8556961.8969221.9381479.201.6211.9826750.1232881.8103141.8515391.8927649.401.8212.6027000.1232881.5846661.6258921.66711711.503.1315.9327250.1232881.3610991.4023241.44354914.305.1619.6927750.1232880.9200510.9612771.00250220.5012.4128.3628000.1232880.7024990.7437240.78495025.8518.1833.2628250.1232880.4868810.5281060.56933131.5025.7138.6128500.1232880.2731620.3143870.35561239.0035.1944.5222000.2191784.9882695.0432375.0982042.600.000.0324000.2191783.4052983.4602653.5152325.800.021.1625000.2191782.6626362.7176032.7725709.000.224.1125400.2191782.3738572.4288242.48379111.000.546.2825500.2191782.3023732.3573402.41230711.950.676.9425750.2191782.1248822.1798492.23481613.201.108.8026000.2191781.9491062.0040732.05904014.801.7610.9926500.2191781.6025681.6575351.71250219.954.1316.4627250.2191781.0948321.1497991.20476629.9212.2327.7428250.2191780.4391680.4941350.54910351.1037.7849.9128500.2191780.2788790.3338460.38881458.0047.6857.1028650.2191780.1833790.2383460.29331465.2054.3661.8024250.2958902.8195132.8833792.94724510.760.163.1625000.2958902.3425892.4064552.47032114.750.706.5426750.2958901.2832051.3470711.41093732.509.9124.3728500.2958900.2909780.3548440.41871069.7955.3264.9128600.2958900.2361340.3000000.36386676.2259.7068.22According to simple computing, the expectation and variance of absolute errors follow as(26)EP1−P^=7.850124,VarP1−P^=29.945065,EP2−P^=3.492115,VarP2−P^=6.764690,where P is the actual closing price of put option, P1 is the value of put option computed by the Black–Scholes model, and P2 is the value of put option computed by Theorem 3. It is obvious that EP2−P^<EP1−P^. In the following, we will support the statement by the hypothesis test (i.e., t-test function of MATLAB R2018b).In fact, for the hypothesis test,(27)H0:EP2−P−EP1−P≥0,H1:EP2−P−EP1−P<0,Thet−statistics equals −4.4919 with degree of freedom 38, and its p value is 0.00019. Thus, we accept H1, i.e., EP2−P<EP1−P. That is, the prices of put options computed by Theorem 3 are far nearer to their actual prices than those computed by the Black–Scholes model.Therefore, our pricing formulas in Theorems2 and 3 have less absolute errors than those of the Black–Scholes model for both call options and put options. That is, the fitting effect of our pricing formulas is superior to that of the Black–Scholes model.


## 5.1. Parameter Estimation

We consider the daily logarithmic returns of S&P 500 index closing prices from January 2015 to December 2018. The augmented Dickey–Fuller test shows the Dickey–Fuller statistic equals−10.283 with lag order 10 and p value is 0.01, which indicates that the time series of returns is stationary. Furthermore, the autocorrelogram shows the coefficients of autocorrelation mostly fall within double standard deviations (see Figure 1). Thus, we accept that the time series of daily logarithmic returns of S&P 500 index closing prices are stationary and independent. We take the annual average return rate based on 252 effective trading days as the drift coefficient of the underlying asset μ.Figure 1
The autocorrelogram of the daily logarithmic returns of S&P 500 index closing prices.There are 41 call options on April 1, 2019, and the current price of underlying assetS0=2867.19. According to simple computation, we obtain μ=4.932% and the left expiration time τ=T−t/365 (years), where t=2019/04/01, and then we work out d1, d2, and d3 by Theorem 2; see Table 5 for detailed data, where C is the actual closing price of call option, τ is the left expiration time of call option, K is the strike price of call option, and d1,d2, and d3 are the parameters in Theorem 2. Furthermore, we obtain the estimated value of the risk-compensation factor λ=−0.0246533 by the least square method with R2=0.9532.Table 5
Parameter estimation ofd1, d2, and d3 for call options.
CτKd1d2d324.400.04657528700.0393200.0646590.08999721.500.0465752875−0.029375−0.0040370.02130210.600.0465752900−0.371070−0.345732−0.3203934.500.0465752925−0.709832−0.684493−0.6591553.850.0465752930−0.777237−0.751898−0.7265593.140.0465752935−0.844526−0.819188−0.7938491.670.0465752950−1.045710−1.020372−0.9950330.650.0465752975−1.378754−1.353416−1.3280770.600.0465752980−1.445027−1.419689−1.39435046.150.12602728700.1047710.1464520.18813344.000.12602728750.0630100.1046910.14637131.000.1260272900−0.144713−0.103032−0.06135119.350.1260272925−0.350653−0.308972−0.2672917.000.1260272975−0.757303−0.715622−0.6739413.730.1260273000−0.958073−0.916392−0.8747113.300.1260273005−0.998026−0.956345−0.9146642.500.1260273015−1.077733−1.036052−0.99437166.700.22191828700.1525040.2078140.26312363.800.22191828750.1210330.1763430.23165248.600.2219182900−0.0355050.0198050.07511436.700.2219182925−0.190699−0.135390−0.08008028.320.2219182940−0.283181−0.227871−0.17256127.700.2219182950−0.344573−0.289263−0.23395413.030.2219183000−0.648446−0.593137−0.5378276.000.2219183050−0.947297−0.891987−0.8366783.100.2219183100−1.241288−1.185978−1.1306690.850.2219183200−1.815306−1.759996−1.70468628.300.2986302975−0.377841−0.313680−0.24951921.900.2986303000−0.508267−0.444106−0.37994517.000.3753423050−0.637911−0.565980−0.49404842.200.4712333000−0.313762−0.233165−0.15256715.100.4712333100−0.720596−0.639998−0.559401133.500.72054828750.2794200.3790830.47874669.700.7205483000−0.147614−0.0479500.05171333.360.7205483100−0.476619−0.376956−0.27729227.500.7205483125−0.557212−0.457549−0.35788515.700.7205483200−0.795178−0.695515−0.5958527.000.7205483300−1.103934−1.004271−0.90460775.300.7972603000−0.109289−0.0044550.10038018.700.7972603200−0.724912−0.620077−0.51524352.440.9698633100−0.319343−0.203716−0.088089There are 62 put options on April 1, 2019, and the current price of underlying assetS0=2867.19. Similarly, we obtain μ=4.932% and the left expiration time τ=T−t/365 (years), where t=2019/04/01, and then we work out d1, d2, and d3 by Theorem 3; see Table 6 for detailed data, where P is the actual closing price of put option, τ is the left expiration time of put option, K is the strike price of put option, and d1,d2, and d3 are the parameters in Theorem 3. Furthermore, we obtain the estimated value of the risk-compensation factor λ=0.0269767 by the least square method with R2=0.9561.Table 6
Parameter estimation ofd1, d2, and d3 for put options.
PτKd1d2d30.350.04657524007.0974587.1227977.1481360.420.04657524106.9333616.9586996.9840380.400.04657524206.7699426.7952816.8206190.400.04657524356.5260766.5514156.5767530.470.04657524406.4451216.4704606.4957990.460.04657524506.2837086.3090476.3343850.550.04657524556.2032486.2285876.2539260.540.04657524855.7239055.7492435.7745820.620.04657524905.6445775.6699155.6952540.560.04657525005.4863985.5117375.5370760.720.04657525155.2503135.2756525.3009900.740.04657525255.0937045.1190425.1443810.750.04657525504.7048784.7302174.7555550.800.04657526003.9385343.9638723.9892110.930.04657526053.8627113.8880503.9133891.380.04657526403.3359953.3613343.3866731.500.04657526503.1867873.2121263.2374641.800.04657526752.8162172.8415552.8668942.440.04657527002.4490932.4744322.4997702.550.04657527052.3760762.4014152.4267542.650.04657527102.3031952.3285332.3538723.030.04657527202.1578332.1831722.2085113.250.04657527252.0853532.1106922.1360304.500.04657527501.7249351.7502741.7756120.550.12602720008.7697708.8114518.8531321.150.12602722755.6788365.7205175.7621981.500.12602723504.9006554.9423364.9840171.950.12602724004.3955464.4372274.4789072.200.12602724254.1469244.1886054.2302853.200.12602725003.4161523.4578333.4995145.100.12602725752.7069832.7486642.7903456.050.12602726002.4751762.5168572.5585386.380.12602726102.3830772.4247582.4664398.550.12602726502.0181752.0598562.10153710.310.12602726751.7928981.8345791.87626012.350.12602727001.5697171.6113981.65307914.350.12602727251.3485931.3902741.43195518.070.12602727501.1294881.1711691.21285022.540.12602727750.9123660.9540470.99572827.700.12602728000.6971910.7388720.78055332.850.12602728250.4839290.5256100.56729141.250.12602728500.2725470.3142270.3559081.300.22191820006.6823806.7376906.7930005.390.22191823753.5753203.6306303.6859395.800.22191824003.3859993.4413083.4966189.260.22191825002.6479352.7032452.75855411.740.22191825502.2899032.3452122.40052212.900.22191825652.1838612.2391712.29448013.550.22191825752.1135112.1688202.22413015.320.22191826001.9388231.9941332.04944217.700.22191826251.7658071.8211171.87642619.600.22191826501.5944311.6497401.70505026.230.22191827001.2564761.3117861.36709528.700.22191827151.1563101.2116191.26692931.250.22191827251.0898391.1451481.20045835.950.22191827500.9247230.9800331.03534246.700.22191828000.5989480.6542570.70956759.320.22191828500.2789390.3342480.38955863.550.22191828600.2156110.2709200.32623064.700.22191828650.1840300.2393400.29464946.600.29863027500.8478740.9120350.97619657.500.29863028000.5670410.6312030.695364


## 5.2. Comparison of Pricing Effect

In the section, we will compare the fitting effect of our pricing formulas with the Black–Scholes model.


### 5.2.1. Comparison of Pricing Effect on Call Options

There are 25 call options on April 2, 2019, and the current price of underlying assetS0=2867.24. Based on the estimated parameters μ=4.932%, σ=0.1174, r=2.42%, and λ=−0.0246533 for call options, the call options of S&P 500 on April 2, 2019, are priced by our obtained pricing formula in Theorem 1 and the Black–Scholes model, respectively; see Table 7 for detailed data, where K is the strike price of call option, τ is the left expiration time of call option, d1,d2, and d3 are the parameters in Theorem 3, C is the actual closing price of call option, C1 is the value of call option computed by the Black–Scholes model, and C2 is the value of call option computed by Theorem 1.Table 7
Pricing results by our pricing formula and the Black–Scholes model for call option.
Kτd1d2d3CC1C228700.0438360.0365110.0610930.08567521.0528.2630.7728750.043836−0.034298−0.0097160.01486618.1525.8526.2528800.043836−0.104985−0.080403−0.05582118.2023.5821.9928900.043836−0.245990−0.221408−0.19682613.7019.4514.2729000.043836−0.386509−0.361927−0.3373458.5015.877.6528700.1232880.1035320.1447580.18598345.0050.0654.7728750.1232880.0613100.1025350.14376040.8547.5550.9128950.123288−0.106850−0.065625−0.02439931.6838.3636.7329000.123288−0.148708−0.107483−0.06625829.5036.2733.4929250.123288−0.356924−0.315698−0.27447318.1027.0019.2629400.123288−0.481000−0.439775−0.39854914.1022.3512.2529500.123288−0.563367−0.522141−0.48091611.5019.618.1929550.123288−0.604445−0.563220−0.52199510.4518.346.3528700.2191780.1516570.2066240.26159164.0069.1475.3928750.2191780.1199900.1749570.22992463.8066.5871.7829000.219178−0.0375230.0174440.07241147.5054.7354.9829500.219178−0.348517−0.293550−0.23858326.7535.5827.8130000.219178−0.654284−0.599317−0.54435013.3021.948.9930250.219178−0.805261−0.750294−0.6953278.7516.882.4229000.2958900.0186620.0825280.14639462.7767.3269.1529000.3726030.0620390.1337080.20537676.3878.7481.7831000.372603−0.868516−0.796848−0.7251799.4018.993.5029000.4684930.1059470.1863100.26667391.4091.8796.1030000.468493−0.315908−0.235545−0.15518241.2652.2244.3831000.468493−0.723929−0.643566−0.56320315.0027.1312.45According to simple computing, the expectation and variance of absolute errors follow as(24)EC1−C^=6.902386,VarC1−C^=6.810772,EC2−C^=5.156575,VarC2−C^=9.726556,where C is the actual closing price of call option, C1 is the value of call option computed by the Black–Scholes model, and C2 is the value of call option computed by Theorem 2. It is obvious that EC2−C^<EC1−C^. In the following, we will support the statement by the hypothesis test (i.e., test function of MATLAB R2018b).In fact, for the hypothesis test,(25)H0:EC2−C−EC1−C≥0,H1:EC2−C−EC1−C<0,the t−statistics equals −2.1465 with degree of freedom 24, and its p value is 0.01854. Thus, we accept H1, i.e., EC2−C<EC1−C. That is, the prices of call options computed by Theorem 2 are far nearer to their actual prices than those computed by the Black–Scholes model.


### 5.2.2. Comparison of Pricing Effect on Put Options

There are 39 put options on April 2, 2019, and the current price of underlying assetS0=2867.24. Based on the estimated parameters μ=4.932%, σ=0.1174, r=2.42%, and λ=0.0269767 for put options, the put options of S&P 500 on April 2, 2019, are priced by our obtained pricing formula in Theorem 3 and the Black–Scholes model, respectively; see Table 8 for detailed data, where K is the strike price of put option, τ is the left expiration time of put option, d1,d2, and d3 are the parameters in Theorem 3, P is the actual closing price of put option, P1 is the value of put option computed by the Black–Scholes model, and P2 is the value of put option computed by Theorem 3.Table 8
Pricing results by our pricing formula and the Black–Scholes model for put option.
Kτd1d2d3PP1P225500.0438364.8456584.8702404.8948220.590.000.0525750.0438364.4487764.4733584.4979400.700.000.1426000.0438364.0557294.0803114.1048930.810.000.3526500.0438363.2808463.3054283.3300101.200.011.5926600.0438363.1276253.1522073.1767901.400.022.0826750.0438362.8988712.9234532.9480351.640.043.0227000.0438362.5204492.5450312.5696131.950.145.2723000.1232885.4740865.5153115.5565361.300.000.0124000.1232884.4417204.4829464.5241711.950.000.1425000.1232883.4515043.4927303.5339553.260.011.1125900.1232882.5936062.6348322.6760575.400.204.6526000.1232882.5001312.5413562.5825816.000.275.3226150.1232882.3605892.4018142.4430406.200.416.4526250.1232882.2680052.3092302.3504567.100.537.2926700.1232881.8556961.8969221.9381479.201.6211.9826750.1232881.8103141.8515391.8927649.401.8212.6027000.1232881.5846661.6258921.66711711.503.1315.9327250.1232881.3610991.4023241.44354914.305.1619.6927750.1232880.9200510.9612771.00250220.5012.4128.3628000.1232880.7024990.7437240.78495025.8518.1833.2628250.1232880.4868810.5281060.56933131.5025.7138.6128500.1232880.2731620.3143870.35561239.0035.1944.5222000.2191784.9882695.0432375.0982042.600.000.0324000.2191783.4052983.4602653.5152325.800.021.1625000.2191782.6626362.7176032.7725709.000.224.1125400.2191782.3738572.4288242.48379111.000.546.2825500.2191782.3023732.3573402.41230711.950.676.9425750.2191782.1248822.1798492.23481613.201.108.8026000.2191781.9491062.0040732.05904014.801.7610.9926500.2191781.6025681.6575351.71250219.954.1316.4627250.2191781.0948321.1497991.20476629.9212.2327.7428250.2191780.4391680.4941350.54910351.1037.7849.9128500.2191780.2788790.3338460.38881458.0047.6857.1028650.2191780.1833790.2383460.29331465.2054.3661.8024250.2958902.8195132.8833792.94724510.760.163.1625000.2958902.3425892.4064552.47032114.750.706.5426750.2958901.2832051.3470711.41093732.509.9124.3728500.2958900.2909780.3548440.41871069.7955.3264.9128600.2958900.2361340.3000000.36386676.2259.7068.22According to simple computing, the expectation and variance of absolute errors follow as(26)EP1−P^=7.850124,VarP1−P^=29.945065,EP2−P^=3.492115,VarP2−P^=6.764690,where P is the actual closing price of put option, P1 is the value of put option computed by the Black–Scholes model, and P2 is the value of put option computed by Theorem 3. It is obvious that EP2−P^<EP1−P^. In the following, we will support the statement by the hypothesis test (i.e., t-test function of MATLAB R2018b).In fact, for the hypothesis test,(27)H0:EP2−P−EP1−P≥0,H1:EP2−P−EP1−P<0,Thet−statistics equals −4.4919 with degree of freedom 38, and its p value is 0.00019. Thus, we accept H1, i.e., EP2−P<EP1−P. That is, the prices of put options computed by Theorem 3 are far nearer to their actual prices than those computed by the Black–Scholes model.Therefore, our pricing formulas in Theorems2 and 3 have less absolute errors than those of the Black–Scholes model for both call options and put options. That is, the fitting effect of our pricing formulas is superior to that of the Black–Scholes model.


## 5.2.1. Comparison of Pricing Effect on Call Options

There are 25 call options on April 2, 2019, and the current price of underlying assetS0=2867.24. Based on the estimated parameters μ=4.932%, σ=0.1174, r=2.42%, and λ=−0.0246533 for call options, the call options of S&P 500 on April 2, 2019, are priced by our obtained pricing formula in Theorem 1 and the Black–Scholes model, respectively; see Table 7 for detailed data, where K is the strike price of call option, τ is the left expiration time of call option, d1,d2, and d3 are the parameters in Theorem 3, C is the actual closing price of call option, C1 is the value of call option computed by the Black–Scholes model, and C2 is the value of call option computed by Theorem 1.Table 7
Pricing results by our pricing formula and the Black–Scholes model for call option.
Kτd1d2d3CC1C228700.0438360.0365110.0610930.08567521.0528.2630.7728750.043836−0.034298−0.0097160.01486618.1525.8526.2528800.043836−0.104985−0.080403−0.05582118.2023.5821.9928900.043836−0.245990−0.221408−0.19682613.7019.4514.2729000.043836−0.386509−0.361927−0.3373458.5015.877.6528700.1232880.1035320.1447580.18598345.0050.0654.7728750.1232880.0613100.1025350.14376040.8547.5550.9128950.123288−0.106850−0.065625−0.02439931.6838.3636.7329000.123288−0.148708−0.107483−0.06625829.5036.2733.4929250.123288−0.356924−0.315698−0.27447318.1027.0019.2629400.123288−0.481000−0.439775−0.39854914.1022.3512.2529500.123288−0.563367−0.522141−0.48091611.5019.618.1929550.123288−0.604445−0.563220−0.52199510.4518.346.3528700.2191780.1516570.2066240.26159164.0069.1475.3928750.2191780.1199900.1749570.22992463.8066.5871.7829000.219178−0.0375230.0174440.07241147.5054.7354.9829500.219178−0.348517−0.293550−0.23858326.7535.5827.8130000.219178−0.654284−0.599317−0.54435013.3021.948.9930250.219178−0.805261−0.750294−0.6953278.7516.882.4229000.2958900.0186620.0825280.14639462.7767.3269.1529000.3726030.0620390.1337080.20537676.3878.7481.7831000.372603−0.868516−0.796848−0.7251799.4018.993.5029000.4684930.1059470.1863100.26667391.4091.8796.1030000.468493−0.315908−0.235545−0.15518241.2652.2244.3831000.468493−0.723929−0.643566−0.56320315.0027.1312.45According to simple computing, the expectation and variance of absolute errors follow as(24)EC1−C^=6.902386,VarC1−C^=6.810772,EC2−C^=5.156575,VarC2−C^=9.726556,where C is the actual closing price of call option, C1 is the value of call option computed by the Black–Scholes model, and C2 is the value of call option computed by Theorem 2. It is obvious that EC2−C^<EC1−C^. In the following, we will support the statement by the hypothesis test (i.e., test function of MATLAB R2018b).In fact, for the hypothesis test,(25)H0:EC2−C−EC1−C≥0,H1:EC2−C−EC1−C<0,the t−statistics equals −2.1465 with degree of freedom 24, and its p value is 0.01854. Thus, we accept H1, i.e., EC2−C<EC1−C. That is, the prices of call options computed by Theorem 2 are far nearer to their actual prices than those computed by the Black–Scholes model.


## 5.2.2. Comparison of Pricing Effect on Put Options

There are 39 put options on April 2, 2019, and the current price of underlying assetS0=2867.24. Based on the estimated parameters μ=4.932%, σ=0.1174, r=2.42%, and λ=0.0269767 for put options, the put options of S&P 500 on April 2, 2019, are priced by our obtained pricing formula in Theorem 3 and the Black–Scholes model, respectively; see Table 8 for detailed data, where K is the strike price of put option, τ is the left expiration time of put option, d1,d2, and d3 are the parameters in Theorem 3, P is the actual closing price of put option, P1 is the value of put option computed by the Black–Scholes model, and P2 is the value of put option computed by Theorem 3.Table 8
Pricing results by our pricing formula and the Black–Scholes model for put option.
Kτd1d2d3PP1P225500.0438364.8456584.8702404.8948220.590.000.0525750.0438364.4487764.4733584.4979400.700.000.1426000.0438364.0557294.0803114.1048930.810.000.3526500.0438363.2808463.3054283.3300101.200.011.5926600.0438363.1276253.1522073.1767901.400.022.0826750.0438362.8988712.9234532.9480351.640.043.0227000.0438362.5204492.5450312.5696131.950.145.2723000.1232885.4740865.5153115.5565361.300.000.0124000.1232884.4417204.4829464.5241711.950.000.1425000.1232883.4515043.4927303.5339553.260.011.1125900.1232882.5936062.6348322.6760575.400.204.6526000.1232882.5001312.5413562.5825816.000.275.3226150.1232882.3605892.4018142.4430406.200.416.4526250.1232882.2680052.3092302.3504567.100.537.2926700.1232881.8556961.8969221.9381479.201.6211.9826750.1232881.8103141.8515391.8927649.401.8212.6027000.1232881.5846661.6258921.66711711.503.1315.9327250.1232881.3610991.4023241.44354914.305.1619.6927750.1232880.9200510.9612771.00250220.5012.4128.3628000.1232880.7024990.7437240.78495025.8518.1833.2628250.1232880.4868810.5281060.56933131.5025.7138.6128500.1232880.2731620.3143870.35561239.0035.1944.5222000.2191784.9882695.0432375.0982042.600.000.0324000.2191783.4052983.4602653.5152325.800.021.1625000.2191782.6626362.7176032.7725709.000.224.1125400.2191782.3738572.4288242.48379111.000.546.2825500.2191782.3023732.3573402.41230711.950.676.9425750.2191782.1248822.1798492.23481613.201.108.8026000.2191781.9491062.0040732.05904014.801.7610.9926500.2191781.6025681.6575351.71250219.954.1316.4627250.2191781.0948321.1497991.20476629.9212.2327.7428250.2191780.4391680.4941350.54910351.1037.7849.9128500.2191780.2788790.3338460.38881458.0047.6857.1028650.2191780.1833790.2383460.29331465.2054.3661.8024250.2958902.8195132.8833792.94724510.760.163.1625000.2958902.3425892.4064552.47032114.750.706.5426750.2958901.2832051.3470711.41093732.509.9124.3728500.2958900.2909780.3548440.41871069.7955.3264.9128600.2958900.2361340.3000000.36386676.2259.7068.22According to simple computing, the expectation and variance of absolute errors follow as(26)EP1−P^=7.850124,VarP1−P^=29.945065,EP2−P^=3.492115,VarP2−P^=6.764690,where P is the actual closing price of put option, P1 is the value of put option computed by the Black–Scholes model, and P2 is the value of put option computed by Theorem 3. It is obvious that EP2−P^<EP1−P^. In the following, we will support the statement by the hypothesis test (i.e., t-test function of MATLAB R2018b).In fact, for the hypothesis test,(27)H0:EP2−P−EP1−P≥0,H1:EP2−P−EP1−P<0,Thet−statistics equals −4.4919 with degree of freedom 38, and its p value is 0.00019. Thus, we accept H1, i.e., EP2−P<EP1−P. That is, the prices of put options computed by Theorem 3 are far nearer to their actual prices than those computed by the Black–Scholes model.Therefore, our pricing formulas in Theorems2 and 3 have less absolute errors than those of the Black–Scholes model for both call options and put options. That is, the fitting effect of our pricing formulas is superior to that of the Black–Scholes model.


## 6. Conclusion

In this study, we obtain the pricing formula of European options, including European call option and European put option, in a risk-aversive market. Corollaries of our obtained results improve the Black–Scholes model owning to its much weaker conditions. It follows from our obtained results that European option value depends on the drift coefficientμ of its underlying security, which does not display in the Black–Scholes model only because μ=r in a risk-neutral financial market according to the no-arbitrage opportunity principle. Empirical analyses show that the fitting effect of our pricing formulas is superior to that of the Black–Scholes model.



---

*Source: 9713521-2021-07-31.xml*