Abstract
We state existence and localisation results for a fully nonlinear boundary value problem using the upper and lower solutions method. With this study we aim to contribute to a better understanding of some analytical features of a problem arising in financial modelling related to the introduction of transaction costs in the classical BlackScholes model. Our result concerns stationary solutions which become interesting in finance when the time does not play a relevant role such as, for instance, in perpetual options.
1 Introduction
In 1973, Fisher Black and Myron Scholes suggested a model that became fundamental for the valuation of financial derivatives in a complete frictionless market. Along with the noarbitrage possibilities, the classical BlackScholes model assumes that in order to replicate exactly the returns of a certain derivative, the hedging portfolio is continuously adjusted by transactioning the underlying asset of the derivative. This fact can only happen if no transaction costs exist when buying or selling financial assets. Otherwise, a continuous adjustment would imply that those costs, such as taxes or fees, would become infinitely large.
Hence the introduction of transaction costs in the model is a problem that has been motivating the work of several authors and has led to the study of new models that generalise the classical BlackScholes model.
In this paper, we aim to give a contribution to better understanding of some analytical features of the problem. We are concerned with the existence and localisation results for the nonlinear secondorder Dirichlet boundary problem
where and p, q are positive constants (p and q can be assumed nonnegative, but the assumption is neither interesting from the mathematical viewpoint nor reasonable for applications to finance). We also consider that . This assumption turns out to be quite natural in some financial settings, for instance, if we are dealing with call options.
This problem is related to the study of stationary solutions of a nonlinear parabolic equation that models the valuation of a call option in presence of transaction costs. These stationary solutions give the option value V as a function of the stock price, which can be interesting when dealing with a model where the time does not play a relevant role such as, for instance, in perpetual options.
Since with , the equation of (1) is of the following type:
However, we will consider it in the form presented in (1) in order to keep some similarity to the financial setting problem, as we will see later.
The arguments used rely on the upper and lower solutions method (see [1,2]).
The results we present were motivated by the work of Amster et al. (see [3]). There, the authors studied an analogous problem but without obtaining a localisation result. Then they proceeded their study using a limit of a nonincreasing (nondecreasing) sequence of upper (respectively lower) solutions. Those results were proved under some Lipschitz condition, which has some financial implications. We do not assume such a condition and, moreover, we obtain localisation information on the solution.
In this paper, we state some existence results for the problem (1), following a previous study contained in [4], and improving the localisation result. Relating to the financial model, this implies more accuracy in the information of the option regarding the behaviour of the stock value, in a model where transaction costs are considered. In Section 2, we consider an auxiliary problem. Using the results proved there, we state in Section 3 an existence and localisation theorem for (1). In Section 4, we make some comments on the financial problem that suggested our analytical study. With this section we aim to illustrate the interplay between finance and mathematics which, in fact, is very challenging and often stimulates the use of innovative mathematical and computational techniques.
2 Auxiliary result
Observe that the equation of (1) is from the algebraic point of view a secondorder equation in the variable . So, solving it algebraically in order of , we obtain the equivalent form
which leads us to consider the equation
where
This fact suggests the study of the auxiliary problem
where
The main argument to prove Theorem 2 relies on the method of upper and lower solutions.
We recall that is a lower solution of (2) if
Similarly, an upper solution of (2) is defined by reversing the inequalities in (3). A solution of (2) is a function u which is simultaneously a lower and an upper solution. A function f is said to satisfy the Nagumo condition on some given subset if there exists a positive continuous function , , such that
and
Lemma 1Suppose that
Then:
1. The functionis a lower solution of the problem (2).
2. Ifis small enough, then the function
is a lower solution of the problem (2).
3. The function
is an upper solution of the problem (2).
Proof 1. If we plug in the first member of the equation in (2), we have that
and by (5)
2. Consider now the function
Observe first that
and
On the other hand, we have that
So, if we plug the function in the first member of the equation in (2), we obtain
which, for k small enough, is nonnegative.
3. Consider now the function
Then, in an analogous way, we plug in the first member of the equation in (2). We observe that
and
So, the three assertions of the lemma are proved. □
Remark 1 The value is a suitable upper bound for the possible values that can take in the above assertion 2. In fact, for and since , easy computations show that
which is used in the last step of the proof of assertion 2.
Now we state an existence and localisation result for the auxiliary problem (2).
1. The problem (2) has a solutionVsuch that
2. Ifis small enough, the problem (2) has a solutionVsuch that
3. If, the functionis a solution of the problem (2).
Proof By the previous lemma, we know already that there are lower and upper solutions for the problem (2). It is also clear that they are well ordered, that is,
So, if the function H satisfies the Nagumo condition, the thesis will follow by [1] or [2].
In order to prove assertions 1 and 2, we consider the sets, respectively,
and
Function H satisfies the Nagumo condition in and in . In fact, we have that
Thus, for some positive constants and , we have that
and
So, by what was said above, the first two assertions of the thesis hold.
The third assertion follows directly from obvious computations. □
Proposition 1Consider the problem (2) and the solutionVgiven by Theorem 1. ThenVis convex and satisfiesin.
Proof The convexity of the solution V follows easily by the equation of the problem (2)
Noting that
then
But, using Theorem 1, assertion 1, we have that . In fact, since , it follows that , for . Letting , we obtain that .
Then we have
□
3 Existence and localisation result
We return to the original problem (1) and state an existence and uniqueness result. We also provide information on the localisation of the solution.
Theorem 2Consider the nonlinear Dirichlet boundary value problem (1). The following assertions hold:
1. The functionis a solution of the problem (1) if and only if.
2. If, then the problem (1) has a convex solutionVsuch that
3. Ifandis small enough, then the problem (1) has a convex solutionVsuch that
4. Moreover, Vis the unique convex solution of (1) in any of the above cases.
Proof 1. It is easily verified that satisfies the equation of (1) and , and that the boundary condition results if and only if .
2. Let denote a convex solution of the auxiliary problem (2) given by Theorem 1. Then, using (2) and Proposition 1, we have that
therefore is a solution of the equation
Since and , we conclude that is a solution of the Dirichlet boundary value problem (1). Statement (6) is obtained from Theorem 1.
3. The proof is the same as for the previous case.
4. The uniqueness result follows immediately from Theorem 2.1 of [3]. □
Remark 2 We remark that Theorem 2.1 of [3] guarantees the existence and uniqueness of a solution but does not give information on localisation, due to the method used. However, using that uniqueness result, from the above theorem, assertions 2 and 3, we can derive the following corollary.
Corollary 1Let. Ifis small enough, then the nonlinear Dirichlet boundary value problem (1) has a convex solutionVsuch that
Figure 1 illustrates the localisation for of the convex solution for the problem (1) whose existence is given by assertion 2 of Theorem 2.
Figure 1. Localisation of the convex solution for the problem (1) stated by Theorem 2, assertion 2.
Assertion 3 of Theorem 2 improves the localisation result contained in assertion 2 as long as k is small enough. In fact, more accuracy is obtained near c.
Figure 2 is an illustration of the new localisation of the convex solution in this case.
Figure 2. Localisation of the convex solution for the problem (1) stated by Corollary 1.
4 Application to a generalised BlackScholes model
For better understanding how that above equation appears from financial option pricing, we present in this section some of the financial and mathematical features that lie behind the model. Namely, we will refer to the introduction of transaction costs in the classical BlackScholes (BS) model. The classical BS model concerns the price of a call or a put option on an underlying asset when the price of this asset itself is modelled as a geometric Brownian motion. We consider the BS equation
with the final condition, which represents the payoff function,
where V is the option price, S is the stock price at time t, T is the maturity time, r is the interest short rate and σ is the volatility of the stock returns.
When approaching the question of introducing transaction costs in the classical BlackScholes model, one of the first main references is the work presented by Leland in [5], that suggests a market with proportional transaction costs. That is, considering ν the number of shares (it is positive if the agent buys or negative if the agent sells) and S the price of the asset at time t, the costs of the transaction of ν shares at time t are given by
where the constant depends on the parts involved in the transaction.
Leland [5] suggested a new strategy using transaction costs proportional to the monetary value of any buy or sale of the asset. Leland’s replication strategy consisted in using the common BlackScholes formulae in periodical revisions of the portfolio but with an appropriately enlarged volatility. This is a model widely accepted in the financial industry. However, there are some mathematical problems with this approach, as referred by Kabanov and Safarian [6]: the terminal value of the replicating portfolio does not converge to the terminal payoff of the derivative if the transaction costs do not depend on the number of revisions (tending to infinity), limiting discrepancy that can be calculated explicitly.
In this work we use a slightly different structure for the transaction costs that was presented by Amster et al.[3]. We assume that the individual cost of the transaction of each share diminishes as the number of shares transactioned increases, which is represented by considering the cost as the percentage given by
where ν is the number of shares traded ( in a buy or in a sale) and are constants depending on the individual investor.
We consider the hedging portfolio consisting in an option of value V (short position) and Δ shares of the underlying asset of price S (long position), and consider a strategy where the portfolio is reviewed every δt, with δt a finite, fixed time step.
The inclusion of transaction costs in the classical BlackScoles model leads us to the following equation: (see [4])
where μ is the drift coefficient and σ is the volatility of the underlying asset.
Considering a small enough such that
the following nonlinear equation is obtained:
The equation obtained (10) is clearly an extension of the classical BlackScholes equation. The introduction of our particular model for the transaction costs in the option pricing market led us to a partial differential equation that contains the BlackScholes terms with an additional nonlinear term modelling the presence of transaction costs. We also use an adjusted volatility in the model , not the real volatility σ.
When concerned with the existence and localisation of stationary solutions of the above equation, in order to consider models where the influence of time does not play an important role, we are led to the following equation:
Dividing by the positive constant and writing and , we obtain
This equation is precisely the equation of the problem (1) studied in the previous section. Contributions to better understanding of the analytical features of the model can be useful for future developments of the financial model itself. In fact, research is full of many examples that illustrate the fruitful interplay between finance and mathematics with reciprocal advantages.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
Both authors collaborated in all the steps concerning the research, discussion and achievements presented in the final manuscript.
Acknowledgements
This research was supported by the Fundação para a Ciência e Tecnologia through the project PEstOE/EGE/UI0491/2011 and by the European Union in the FP7PEOPLE2012ITN Program under Grant Agreement Number 304617 (FP7 Marie Curie Action, Project MultiITN STRIKE  Novel Methods in Computational Finance. We thank the anonymous referees for their useful comments.
References

Mawhin, J: Points fixes, points critiques et problèmes aux limites, Presses University Montreal, Montreal (1985)

De Coster, C, Habets, P: TwoPoint Boundary Value Problems Lower and Upper Solutions, Elsevier, Amsterdam (2006)

Amster, P, Averbuj, CG, Mariani, MC, Rial, D: A BlackScholes option pricing model with transaction costs. J. Math. Anal. Appl.. 303, 688–695 (2005). Publisher Full Text

Grossinho, MR, Morais, E: A note on a stationary problem for a BlackScholes equation with transaction costs. Int. J. Pure Appl. Math.. 51, 579–587 (2009). PubMed Abstract

Leland, HE: Option pricing and replication with transaction costs. J. Finance. 40, 1283–1301 (1985). Publisher Full Text

Kabanov, Y, Safarian, M: Markets with Transaction Costs  Mathematical Theory, Springer, Berlin (2009)