Research Article Empirical Pricing of Chinese Defaultable Corporate Bonds Based on the Incomplete Information Model

Similar documents
Structural Models of Credit Risk and Some Applications

Research Article The Volatility of the Index of Shanghai Stock Market Research Based on ARCH and Its Extended Forms

Introduction Credit risk

Pricing Convertible Bonds under the First-Passage Credit Risk Model

A No-Arbitrage Theorem for Uncertain Stock Model

Theoretical Problems in Credit Portfolio Modeling 2

Time-changed Brownian motion and option pricing

Research Article The European Vulnerable Option Pricing with Jumps Based on a Mixed Model

A Simple Model of Credit Spreads with Incomplete Information

No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate

NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours

Models for Credit Risk in a Network Economy

MODELING DEFAULTABLE BONDS WITH MEAN-REVERTING LOG-NORMAL SPREAD: A QUASI CLOSED-FORM SOLUTION

American Option Pricing Formula for Uncertain Financial Market

Research Article Risk Measurement for Portfolio Credit Risk Based on a Mixed Poisson Model

Modeling Credit Risk with Partial Information

Option Pricing under Delay Geometric Brownian Motion with Regime Switching

Asset Pricing Models with Underlying Time-varying Lévy Processes

Valuing Coupon Bond Linked to Variable Interest Rate

Research Article A Novel Machine Learning Strategy Based on Two-Dimensional Numerical Models in Financial Engineering

Option Pricing Formula for Fuzzy Financial Market

Credit Derivatives and Risk Aversion

Dynamic Portfolio Choice II

Research Article A Mathematical Model of Communication with Reputational Concerns

Credit Risk and Underlying Asset Risk *

A NOVEL BINOMIAL TREE APPROACH TO CALCULATE COLLATERAL AMOUNT FOR AN OPTION WITH CREDIT RISK

Lecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r.

Chapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets

Financial Risk Management

A Comparison of Credit Risk Models

M5MF6. Advanced Methods in Derivatives Pricing

The term structure model of corporate bond yields

arxiv: v2 [q-fin.pr] 23 Nov 2017

Asset-based Estimates for Default Probabilities for Commercial Banks

Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.

Barrier Options Pricing in Uncertain Financial Market

Distortion operator of uncertainty claim pricing using weibull distortion operator

Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing

Stochastic Volatility (Working Draft I)

Pricing of Futures Contracts by Considering Stochastic Exponential Jump Domain of Spot Price

Economathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t

Optimal stopping problems for a Brownian motion with a disorder on a finite interval

Self-Exciting Corporate Defaults: Contagion or Frailty?

Pricing Dynamic Guaranteed Funds Under a Double Exponential. Jump Diffusion Process. Chuang-Chang Chang, Ya-Hui Lien and Min-Hung Tsay

Credit Risk Models with Filtered Market Information

Research Article Welfare Comparison of Leader-Follower Models in a Mixed Duopoly

Modelling the Term Structure of Hong Kong Inter-Bank Offered Rates (HIBOR)

Slides for Risk Management Credit Risk

An overview of some financial models using BSDE with enlarged filtrations

Stochastic Processes and Stochastic Calculus - 9 Complete and Incomplete Market Models

Decomposing swap spreads

Valuation of Defaultable Bonds Using Signaling Process An Extension

Modelling Credit Spread Behaviour. FIRST Credit, Insurance and Risk. Angelo Arvanitis, Jon Gregory, Jean-Paul Laurent

Pricing of a European Call Option Under a Local Volatility Interbank Offered Rate Model

1 The continuous time limit

Research on the Determinants of China s Corporate Bond Credit Spreads

AMH4 - ADVANCED OPTION PRICING. Contents

Credit Risk: Modeling, Valuation and Hedging

The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations

A Note about the Black-Scholes Option Pricing Model under Time-Varying Conditions Yi-rong YING and Meng-meng BAI

MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS.

TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING

Youngrok Lee and Jaesung Lee

The Black-Scholes Model

The Actuary Pricing of an Innovative Housing Mortgage Insurance

The Double Skorohod Map and Real-Time Queues

A r b i t r a g e C r a s h e s i n t h e C o n v e rt i b l e B o n d M a r k e t

EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS

Drawdowns Preceding Rallies in the Brownian Motion Model

Research Article Robust Stability Analysis for the New Type Rural Social Endowment Insurance System with Minor Fluctuations in China

AN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL

25857 Interest Rate Modelling

Estimation of Default Risk in CIR++ model simulation

Portfolio optimization problem with default risk

Credit Risk modelling

Contagion models with interacting default intensity processes

Unified Credit-Equity Modeling

The Black-Scholes Model

Valuing Early Stage Investments with Market Related Timing Risk

Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.

BIRKBECK (University of London) MSc EXAMINATION FOR INTERNAL STUDENTS MSc FINANCIAL ENGINEERING DEPARTMENT OF ECONOMICS, MATHEMATICS AND STATIS- TICS

DOI: /s Springer. This version available at:

Whether Cash Dividend Policy of Chinese

CMBS Default: A First Passage Time Approach

Hedging with Life and General Insurance Products

CREDIT RATINGS. Rating Agencies: Moody s and S&P Creditworthiness of corporate bonds

The stochastic calculus

Exact Sampling of Jump-Diffusion Processes

the Role of Contingent Convertible Bond in Capital Structure Decisions Weikeng Chen [401248]

MORNING SESSION. Date: Wednesday, April 30, 2014 Time: 8:30 a.m. 11:45 a.m. INSTRUCTIONS TO CANDIDATES

Credit Risk : Firm Value Model

STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL

CDS Pricing Formula in the Fuzzy Credit Risk Market

arxiv:cond-mat/ v2 [cond-mat.str-el] 5 Nov 2002

Research Article Portfolio Selection with Subsistence Consumption Constraints and CARA Utility

Financial Derivatives Section 5

Systemic Influences on Optimal Investment

Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models

Risk Minimization Control for Beating the Market Strategies

Credit Portfolio Risk

Transcription:

Mathematical Problems in Engineering, Article ID 286739, 5 pages http://dx.doi.org/10.1155/2014/286739 Research Article Empirical Pricing of Chinese Defaultable Corporate Bonds Based on the Incomplete Information Model Li Ping and Wang Xiaoxu Department of Finance, Beihang University, Beijing 100191, China Correspondence should be addressed to Li Ping; lipingxx@126.com Received 29 December 2013; Accepted 9 February 2014; Published 24 March 2014 Academic Editor: Chuangxia Huang Copyright 2014 L. Ping and W. Xiaoxu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The default of Suntech Power made the year 2013 in China the first year of default of bond markets. People are also clearly aware of the default risk of corporate bonds and find that fair pricing for defaultable corporate bonds is very important. In this paper we first give the pricing model based on incomplete information, then empirically price the Chinese corporate bond 11 super JGBS from Merton s model, reduced-form model, and incomplete information model, respectively, and then compare the obtained prices with the real prices. Results show that all the three models can reflect the trend of bond prices, but the incomplete information model fits the real prices best. In addition, the default probability obtained from the incomplete information model can discriminate the credit quality of listed companies. 1. Introduction Because of the subprime crisis and the European sovereign debt crisis, global economy has been seriously hit. Subsequently, Standard and Poor downgraded the US credit rating. In China, the no default phenomenon in bond markets was broken by the default of the Suntech Power in 2013. So the year 2013 in China became the first year of default of bond markets. The year 2012 was of great significance for Chinese bond market. The United Credit indicated that there were 4 defaults in Chinese interbank bond market during 2012. The overall default probability of corporate bonds was 0.17 if in terms of the bond maturity and 0.12 if in terms of bond issuance. This means that the year 2012 was a new starting point of Chinese bond market. The company Chaori Solar Energy has attracted public attention since the issuance. And due to the high default risk, it was once predicted to be the first default of domestic bonds. Because of the higher and higher default risk for Chinese corporate bond, it is very important to study the issues of pricing and risk management for defaultable bonds. There are three main models for defaultable bond pricing: the traditional structural model [1] and reduced-form model [2] and the emerging incomplete information model. The traditional models are carried out under the assumption of perfect information and that default can be measured by observable changes in the asset value. However, in the real world, the volatility and growth rate of a company s asset value cannot be directly observed, which makes the default probability based on the traditional model biased with the actual situation and also leads to the result that the expected absolute level is low. Therefore, Duffie and Lando [3] relaxed the complete information assumptions in traditional models and proposed the incomplete information model. Then Giesecke and Goldberg [4] and Liu[5] further investigated the incomplete information model under the assumption of incomplete information and gave different types of incomplete information model. In this paper, we give a new pricing formula for defaultable corporate bonds by combining the incomplete information model with the credit risk premium based on default probability. Then we calculate the default probability and price for the Chinese corporate bond 11 super JGBS from the three models: structural model, reduced-form model, and incomplete information model. Results show that all the three models can reflect the trend of bond prices, but the incomplete information model can fit the real prices best and

2 Mathematical Problems in Engineering can discriminate the credit quality of listed companies more effectively. The rest of the paper is organized as follows. In Section 2, we introduce the traditional pricing models for defaultable corporate bonds. Then in Section 3 we propose a new pricing model based on the incomplete information model and the credit risk premium. An empirical calculation is given in Section4,and then Section5 concludes the paper. 2. Traditional Pricing Models for Defaultable Bonds In this section, we will introduce the two traditional pricing models: Merton s structural model and reduced-form model. 2.1. Structural Model. The structural model was first put forward by Merton [1] and then extended by Black and Cox [6], in which the company s debt is regarded as a contingent claim of the company s assets. When the company s asset value is lower than its liabilities, the company will go bankrupt. The value of equity can be obtained by the Black-Scholes option pricing model [7], and the debt value is the company s total value minus the equity value. In this paper, we assume that stocks and bonds constitute the company s assets structure. In Merton s model, if the asset valueislowerthantheliabilityvalueatthematurityofdebt, thecorporatedefaults.heassumedthatthevalueofcorporate bonds basically depends on three factors: the rate of return of a risk-free bond, default probability, and various regulations and restrictions in the bond terms such as maturity date, coupon rate, seniority of default event, and sinking fund. 2.2. Reduced-Form Model. The reduced-form model assumes that the default is exogenous and will occur with no warnings. It directly defines the default time according to default intensity, which avoids to model the unobservable corporate value. The representative work includes Jarrow and Turnbull [2] and Duffie and Singleton[8]. They treated the default as a jump process, and the default process can depend on exogenous macrostate variables. The reduced-form model has the following characteristics: (1) it assumes that the market was complete which is accordant with the no-arbitrage assumption; (2) the credit risk can be obtained by the default probability; (3) the default is a random process; (4) the default recovery rate is an exogenous variable. The reduced-form model does not consider the relationship between the default and corporate value, so some researchers made some extension. In the reduced-form model, default time is decided by the first jump of the exogenous jump process and the parameters of default intensity are given by the market data, and information can be observed in the market. However, the real market is asymmetric, and we can hardly avoid the information distortion. Therefore, in addition to the part information of the company s assets and liabilities, other information such as the default boundaries is difficult to observe. 3. The Incomplete Information Model Combined with Credit Risk Premium Different from the traditional models, the incomplete information model is to analyze credit risk based on incomplete information. In reality, default occurs when there is no omen, and the information about company s asset value and default boundary cannot be completely observed by investors. That is, the information is incomplete. We cannot predict when defaultoccursandcannotpricethebondaccordingtothe exact asset value and default boundaries. Under the condition that the information about company s asset value and boundaries is not complete, investors can only observe the default according to Giesecke [9]. When companies issue bonds, they only announce their initial debt value V 0 to investors. Then after they issue the bonds, investors receive the incomplete information about corporate value such as the company s accounting statements with noise. We use a complete probability space (Ω, I, P)to represent the uncertainty and use I t to denote the information filtration. From the information incompleteness we know that corporate value V t I t.wethenassumethatv t follows a geometric Brownian motion with volatility σ and drift μ = m (1/2)σ 2. M t is used to denote the historical lowest value of the company; then the distribution of M t is Ψ (t, x) =P(M t x) =Φ( x μt (1) )+exp (2μx σ t σ 2 )Φ(x+μt σ t ), where, Φ is the standard normal distribution. Then the derivative of Ψ(t, x) with respect to t is φ (t, x) = 1 2σ [( μ t x 2 x+μt t 3 )e2μx/σ φ( σ t ) ( μ t + x (2) x )φ(μt t3 σ t )], where φ is the standard normal density function. Duffie and Lando [3] showed that when the default boundary D is known, the price trend A t satisfies t A (t) = λ s ds, (3) 0 where λ t = φ(t, d)/(1 Ψ(t, d)) is the default density, d= ln(d). We use τ to denote the default time; then its conditional distribution is T p (t, T) =P(τ T I t )=1 E[exp ( λ s ds I t )] t = Ψ (T, d) Ψ(t, d). 1 Ψ(t, d) For t 0, the conditional distribution of M t with respect to I t is assumed to be H(t, x). From the information (4)

Mathematical Problems in Engineering 3 Table 1: Basic information of the bond. Bond code Bond name Bond issuer Issue date Maturity date Credit rating Coupon rate 112061 11 super JGBS Shanghai Chaori Solar Energy Science 07-03-2012 07-03-2017 A 8.98% Table 2: Estimates of related parameters. Parameters σ e σ V μ e t T Default probability Bond price Estimation 0.61 0.36 0.73 1.32 5 37.37% 102.36 incompleteness we know that H(t, x) is continuous and increasing. Then from Duffie and Lando [3], the conditional survival probability L t is V 0 L t =1 H (t, x) dg (x), (5) where G(x) = e x, x<0.ifthepricetrenda t is continuous, then V 0 A t = log (1 H (t, x) dg (x)). (6) When investors can only observe the default we have Assuming that V 0 =0,then H (t, x) =Ψ(t, x). (7) 0 A t = log [1 Ψ (t, x) e x dx]. (8) Together with (3)weobtain p (s, V) =Φ( V μs ) e V+sm V υs Φ( σ s σ s ) + 1 γ e(1 γ)v Φ( μs V σ s ) 1 γ ev+sβ Φ( δs V σ s ), where υ=μ+σ 2, γ = 1 + (2μ/σ 2 ), δ=μ γσ 2, β= μγ+ γ 2 σ 2 /2. For simplicity we set Γ= t/σ;then From (2)and(9)we get p (t, 0) =Φ( μγ)+e tm Φ ( VΓ) + 1 γ Φ(μΓ) 1 γ etβ Φ (δγ). A t = log (1 p(s,0)), λ t = Then the default probability is p (t, T) =P(τ T ς t )=1 e A t A T =1 e T λ t ds s = (9) (10) p (t, 0) 1 p(t, 0). (11) p (T, 0) p(t, 0). 1 p(t, 0) (12) From the above formulas we can see that, for different hypothesis about the information of the corporate value and default boundary, the default probability and credit risk premium are different. Let r be the risk-free interest rate; let R t and p t be the yield and default probability of the corporate bond at time t, respectively. Then for a single period bond, (1 + R t )(1 p t )= 1+r;thatis, R t =r+ p t (1+r). (13) 1 p t We put s= (p/(1 p))(1+r);thens is the credit risk premium of the corporate bond. For multiperiod corporate bond, we assume that r i (i = 1,2,...,n)is the risk-free interest rate for maturity i, s i is the credit risk premium, and C i isthecashflow;thentheexpected price of a corporate bond is n P= i=1 (1 + r i +s i ). (14) This is the formula for pricing a corporate bond considering the credit risk premium. In the next section we will use this formula to empirically calculate the price for a Chinese corporatebondandcomparethepricesobtainedfrom(14) with those obtained from the traditional structural model and reduced-form model. 4. Empirical Pricing In this section, we choose the 11 super JGBS bond with high default risk for empirical analysis. The 11 super JGBS bond was issued by Shanghai Chaori Solar Energy Science Company. When it was first issued, Pengyuan Credit Bureau rated it as AA; then on December 27, 2012, it was downgraded to AA,,thenonApril10,2013,toBBB+,andfinallyonMay 20, 2013, to CCC. On December 20, 2012, the stocks and bonds of Shanghai Chaori Solar Energy Science were both suspended, so we select the data of 11 super JGBS from April 20, 2012, to December 19, 2012, for empirical research. Table 1 is about the basic information of the bond. Firstly, we calculate the default probability and the price of 11 super JGBS based on the incomplete information model. Under this model, we calculate the asset value V and its volatility σ V and the company stock value E and its volatility σ e as well as its draft term μ e,respectively.thenweusematlab to estimate the parameters and get the results as shown in Table 2. C i

4 Mathematical Problems in Engineering Table 3: Moody s credit rating standard. Credit ratings Aaa Aa A Baa Ba B C Default probability (%) <0.02 0.02 0.04 0.04 0.08 0.08 0.2 0.2 1.8 1.8 8.3 >8.3 Table 3 is Moody s credit rating standard, from which we can see that, under the incomplete information model, 11 super JGBS bond has a high default probability and is already ajunkbond.butinreality,thebondratingremainsatalevel. Since 2012, the company has experienced a series of credit events such as delinquent loans, stagnation of production, and investigation by the regulators. The suspension of the bond on December 20, 2012, also verifies the truth that the default probability is very high. Investors are very concerned with the bond s credit risk, so the incomplete information model can effectively discriminate the bonds with high risk. Next, we calculate the default probability for the sample period every 5 days based on the asset s volatility under three models and get Figure 1. Combining with Table 3, we can see that, during the sample periods, the default probabilities estimated from Merton s model and the reduced-form model are significantly less than those obtained from the incomplete information model. For thetwotraditionalmodels,thebondcreditratingremains at A level. Only for some points obtained from Merton s model, the estimated credit rating falls to Baa. However, the default probability estimated from the incomplete information model has been rising, and the credit rating is falling. It indicates that the incomplete information model is most sensitive to the change of company information, so it can detect the changes in the credit quality quickly. Therefore the incomplete information model can effectively discriminate thecreditriskoflistedcompanies. We then calculate the bond s theoretical prices for the sample period every 5 days based on the asset s volatility under the three models and compare them with the real prices, which are shown Figure2. From Figure 2 we can see that all the three models can reflect the general trend of bond prices, but in comparison, the incomplete information model fits the real prices and reflects the basic movements of prices better. In addition, influenced by bad information such as reduced performance and credit downgrade, the default risk of the bond is expanding fast. For example, on December 20, the stock of Chaori Solar Energy Science Company was suspended. Based on historical information of the company, the incomplete information model reflects the surge of credit risk. Finally, we make the sensitivity analysis by letting the asset s volatility range from 0.02 to 0.06 and the risk-free interest rate from 0.0314 to 0.072. Results are shown in Figure 3. In Figure 3, the x-, y-, and z-axes represent the asset s volatility, the risk-free interest rate, and the bond price, respectively.wecanseethatthehighertheasset svolatility, the higher the bond price, and the higher the risk-free interest rate, the lower the bond price. And the bond price is more sensitive to the asset s volatility. So the asset s volatility is an important factor influencing the bond price. Default probability 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 20-04-2012 15-05-2012 09-06-2012 04-07-2012 Merton s model Reduced model 29-07-2012 23-08-2012 Time 17-09-2012 12-10-2012 06-11-2012 01-12-2012 Incomplete information model Figure 1: Default probabilities of 11 super JGBS bond under three models. Price 120 115 110 105 100 95 90 85 80 75 20-04-2012 15-05-2012 09-06-2012 04-07-2012 Real price Merton s model 29-07-2012 23-08-2012 Time 17-09-2012 12-10-2012 06-11-2012 01-12-2012 Reduced model Incomplete information model Figure 2: Comparison of the bond s prices.

Mathematical Problems in Engineering 5 Bond price 102 100 98 96 94 92 90 0.020 0.025 0.028 0.031 0.035 0.038 0.042 0.045 0.049 0.052 0.055 0.059 100 102 98 100 96 98 5. Conclusions Asset s volatility 94 96 92 94 90 92 Figure 3: Sensitivity analysis. 0.063 0.053 0.043 0.031 Risk-free interest rate In this paper, we extend the incomplete information model by combining with the credit risk premium and redefined the pricing formula for defaultable corporate bonds. By comparing the three models, we find that the incomplete information model is a better method to identify and distinguish between bond credit ratings. It can also fit bond prices better and reflect the change of company information more sensitively. Suntech s default broke the no-default phenomenon in Chinese bond market, which means that the impact of listed company s credit risk is more and more important. We think that the incomplete information model is more effective to price defaultable corporate bonds and will provide new thoughts for investors to avoid credit risk. In recent years, some researchers use the structural model and reduced-form model to study the default correlation and then price defaultable bonds and credit derivatives, such as Yu [10],P.LiandZ.Z.Li[11], and Wen and Liu [12]. In our next research, we will use the incomplete information model combined with copula function to study the default correlation and price defaultable bonds and credit derivatives. [2] R. Jarrow and S. Turnbull, Pricing derivatives on financial securities subject to credit risk, The Finance, vol.50,no. 1, pp. 53 85, 1995. [3] D. Duffie and D. Lando, Term structures of credit spreads with incomplete accounting information, Econometrica, vol. 69, no. 3,pp.622 664,2001. [4] K. Giesecke and L. Goldberg, Sequential defaults and incomplete information, Risk,vol.7,no. 1,pp.1 26, 2004. [5] X.F.Liu,The researeh on default risk of corporate bond based on incomplete information model [M.S. thesis], Jinan University, Guangzhou, China, 2008. [6] F.BlackandJ.Cox, Valuingcorporatesecurities:someeffects of bond indenture provisions, The Finance, vol. 31, no. 2, pp. 351 367, 1976. [7] F. Black and M. Scholes, The Pricing of Options and Coporate Liabilities, Joumal of Political Eeonomy, no.81,pp.637 654, 1973. [8] D. Duffie and K. J. Singleton, Modeling term structures of defaultable bonds, Review of Financial Studies, vol.12,no.4, pp. 687 720, 1999. [9] K. Giesecke, Default and information, Economic Dynamics and Control,vol.30,no.11,pp.2281 2303,2006. [10] F. Yu, Correlated defaults and the valuation of defaultable securities, in Proceedings of the 2nd International Conference on Credit Risk,pp.1 30,2004. [11] P. Li and Z. Z. Li, Change analysis for the dependence structure and dynamic pricing of basket default swaps, European Financial Management,2013. [12] F. H. Wen and Z. F. Liu, A copula-based correlation measure and its application in chinese stock market, International Journal of Information Technology and Decision Making, vol.8,no. 4, pp. 787 801, 2009. Conflict of Interests The authors declare that there is no conflict of interests regarding the publication of this paper. Acknowledgment This work was supported by the National Natural Science Foundation of China (nos. 71271015, 70971006). References [1] C. R. Merton, On the pricing of corporate debt: the risk structure of interest rates, The Finance, vol.29,no.2,pp. 449 470, 1974.

Advances in Operations Research Advances in Decision Sciences Applied Mathematics Algebra Probability and Statistics The Scientific World Journal International Differential Equations Submit your manuscripts at International Advances in Combinatorics Mathematical Physics Complex Analysis International Mathematics and Mathematical Sciences Mathematical Problems in Engineering Mathematics Discrete Mathematics Discrete Dynamics in Nature and Society Function Spaces Abstract and Applied Analysis International Stochastic Analysis Optimization