Closed form approximations for spread options

This article expresses the price of a spread option as the sum of the prices of two compound options. One compound option is to exchange vanilla call options on the two underlying assets and the other is to exchange the corresponding put options. This way we derive a new closed form approximation for the price of a European spread option and a corresponding approximation for each of its price, volatility and correlation hedge ratios. Our approach has many advantages over existing analytical approximations, which have limited validity and an indeterminacy that renders them of little practical use. The compound exchange option approximation for European spread options is then extended to American spread options on assets that pay dividends or incur costs. Simulations quantify the accuracy of our approach; we also present an empirical application to the American crack spread options that are traded on NYMEX. For illustration, we compare our results with those obtained using the approximation attributed to Kirk (1996, Correlation in energy markets. In: V. Kaminski (Ed.), Managing Energy Price Risk, pp. 71–78 (London: Risk Publications)), which is commonly used by traders.

Exact closed-form solutions of the Dirac equation for a class of effective Morse potentials

Exact bounded solutions for a fermion subject to exponential scalar potential in 1 + 1 dimensions are found in closed form. We discuss the existence of zero modes which are related to the ultrarelativistic limit of the Dirac equation and are responsible for the induction of a fractional fermion number on the vacuum.

Exact closed-form solutions of the Dirac equation with a scalar exponential potential

The problem of a fermion subject to a general scalar potential in a two-dimensional world for nonzero eigenenergies is mapped into a Sturm-Liouville problem for the upper component of the Dirac spinor. In the specific circumstance of an exponential potential, we have an effective Morse potential which reveals itself as an essentially relativistic problem. Exact bound solutions are found in closed form for this problem. The behaviour of the upper and lower components of the Dirac spinor is discussed in detail, particularly the existence of zero modes. (c) 2005 Elsevier B.v. All rights reserved.

Exact Solution of the Elastica with Transverse Shear Effects

A new derivation of Euler's Elastica with transverse shear effects included is presented. The elastic potential energy of bending and transverse shear is set up. The work of the axial compression force is determined. The equation of equilibrium is derived using the variation of the total potential. Using substitution of variables an exact solution is derived. The equation is transcendental and does not have a closed form solution. It is evaluated in a dimensionless form by using a numerical procedure. Finally, numerical examples of laminates made of composite material (fiber reinforced) and sandwich panels are provided.

DROMO formulation for planar motions: Solution to the Tsien Problem

The two-body problem subject to a constant radial thrust is analyzed as a planar motion. The description of the problem is performed in terms of three perturbation methods: DROMO and two others due to Deprit. All of them rely on Hansen?s ideal frame concept. An explicit, analytic, closed-form solution is obtained for this problem when the initial orbit is circular (Tsien problem), based on the DROMO special perturbation method, and expressed in terms of elliptic integral functions. The analytical solution to the Tsien problem is later used as a reference to test the numerical performance of various orbit propagation methods, including DROMO and Deprit methods, as well as Cowell and Kustaanheimo?Stiefel methods.

Some examples of the simultaneous use of closed-form and numerical solutions in shell structural analysis

The computational advantages of the use of different approaches -numerical and analytical ones- to the analysis of different parts of the same shell structure are discussed. Examples of large size problems that can be reduced to those more suitable to be handled by a personal related axisyrometric finite elements, local unaxisymmetric shells, geometric quasi-regular shells, infinite elements and homogenization techniques are described

Closed Form Solutions of Integrable Nonlinear Evolution Equations

2010 Mathematics Subject Classification: 35Q55.

On Queueinginventory Models – product form solution; reservation, Cancellation and common life time

Queueing theory is the mathematical study of ‘queue’ or ‘waiting lines’ where an item from inventory is provided to the customer on completion of service. A typical queueing system consists of a queue and a server. Customers arrive in the system from outside and join the queue in a certain way. The server picks up customers and serves them according to certain service discipline. Customers leave the system immediately after their service is completed. For queueing systems, queue length, waiting time and busy period are of primary interest to applications. The theory permits the derivation and calculation of several performance measures including the average waiting time in the queue or the system, mean queue length, traffic intensity, the expected number waiting or receiving service, mean busy period, distribution of queue length, and the probability of encountering the system in certain states, such as empty, full, having an available server or having to wait a certain time to be served.

Flow of a micropolar fluid bounded by a stretching sheet

We consider boundary layer flow of a micropolar fluid driven by a porous stretching sheet. A similarity solution is defined, and numerical solutions using Runge-Kutta and quasilinearisation schemes are obtained. A perturbation analysis is also used to derive analytic solutions to first order in the perturbing parameter. The resulting closed form solutions involve relatively complex expressions, and the analysis is made more tractable by a combination of offline and online work using a computational algebra system (CAS). For this combined numerical and analytic approach, the perturbation analysis yields a number of benefits with regard to the numerical work. The existence of a closed form solution helps to discriminate between acceptable and spurious numerical solutions. Also, the expressions obtained from the perturbation work can provide an accurate description of the solution for ranges of parameters where the numerical approaches considered here prove computationally more difficult.

An analytical and experimental study of the vibration response of a clamped ribbed plate

An analytical solution is presented in this paper for the vibration response of a ribbed plate clamped on all its boundary edges by employing a travelling wave solution. A clamped ribbed plate test rig is also assembled in this study for the experimental investigation of the ribbed plate response and to provide verification results to the analytical solution. The dynamic characteristics and mode shapes of the ribbed plate are measured and compared to those obtained from the analytical solution and from finite element analysis (FEA). General good agreements are found between the results. Discrepancies between the computational and experimental results at low and high frequencies are also discussed. Explanations are offered in the study to disclose the mechanism causing the discrepancies. The dependency of the dynamic response of the ribbed plate on the distance between the excitation force and the rib is also investigated experimentally. It confirms the findings disclosed in a previous analytical study [T. R. Lin and J. Pan, A closed form solution for the dynamic response of finite ribbed plates. Journal of the Acoustical Society of America 119 (2006) 917-925] that the vibration response of a clamped ribbed plate due to a point force excitation is controlled by the plate stiffness when the source is more than a quarter plate bending wavelength away from the rib and from the plate boundary. The response is largely affected by the rib stiffness when the source location is less than a quarter bending wavelength away from the rib.

An improved computational method in structural dynamics

Purpose – In structural, earthquake and aeronautical engineering and mechanical vibration, the solution of dynamic equations for a structure subjected to dynamic loading leads to a high order system of differential equations. The numerical methods are usually used for integration when either there is dealing with discrete data or there is no analytical solution for the equations. Since the numerical methods with more accuracy and stability give more accurate results in structural responses, there is a need to improve the existing methods or develop new ones. The paper aims to discuss these issues. Design/methodology/approach – In this paper, a new time integration method is proposed mathematically and numerically, which is accordingly applied to single-degree-of-freedom (SDOF) and multi-degree-of-freedom (MDOF) systems. Finally, the results are compared to the existing methods such as Newmark’s method and closed form solution. Findings – It is concluded that, in the proposed method, the data variance of each set of structural responses such as displacement, velocity, or acceleration in different time steps is less than those in Newmark’s method, and the proposed method is more accurate and stable than Newmark’s method and is capable of analyzing the structure at fewer numbers of iteration or computation cycles, hence less time-consuming. Originality/value – A new mathematical and numerical time integration method is proposed for the computation of structural responses with higher accuracy and stability, lower data variance, and fewer numbers of iterations for computational cycles.

Scalable Bayesian Computation for intractable likelihoods in image analysis

The inverse temperature hyperparameter of the hidden Potts model governs the strength of spatial cohesion and therefore has a substantial influence over the resulting model fit. The difficulty arises from the dependence of an intractable normalising constant on the value of the inverse temperature, thus there is no closed form solution for sampling from the distribution directly. We review three computational approaches for addressing this issue, namely pseudolikelihood, path sampling, and the approximate exchange algorithm. We compare the accuracy and scalability of these methods using a simulation study.

On the efficacy of Fourier series approximations for pricing European options

This paper investigates several competing procedures for computing the prices of vanilla European options, such as puts, calls and binaries, in which the underlying model has a characteristic function that is known in semi-closed form. The algorithms investigated here are the half-range Fourier cosine series, the half-range Fourier sine series and the full-range Fourier series. Their performance is assessed in simulation experiments in which an analytical solution is available and also for a simple affine model of stochastic volatility in which there is no closed-form solution. The results suggest that the half-range sine series approximation is the least effective of the three proposed algorithms. It is rather more difficult to distinguish between the performance of the halfrange cosine series and the full-range Fourier series. However there are two clear differences. First, when the interval over which the density is approximated is relatively large, the full-range Fourier series is at least as good as the half-range Fourier cosine series, and outperforms the latter in pricing out-of-the-money call options, in particular with maturities of three months or less. Second, the computational time required by the half-range Fourier cosine series is uniformly longer than that required by the full-range Fourier series for an interval of fixed length. Taken together,these two conclusions make a case for pricing options using a full-range range Fourier series as opposed to a half-range Fourier cosine series if a large number of options are to be priced in as short a time as possible.

Biaxial-load effects on isopachics for a crack at the interface of dissimilar media

The plane problem of two dissimilar materials, bonded together and containing a crack along their common interface, which were subjected to a biaxial load at infinity, is examined by giving a closed-form expression for the first stress invariant of the normal stresses, which is equally valid everywhere, near to, and far from, the crack-tip region. This exact expression for the first-stress invariant is compared by constructing the respective isopachic-fringe patterns, to the approximate expression with non-singular terms, due to the biaxiality factor, for the same quantity. Significant differences between respective isopachic-patterns were found and their dependence on the elastic properties of both materials and the applied loads was demonstrated. The relative errors between the computedK I - andK II -components by using the approximate expression for the first stress-invariant and the accurate one, derived from closed-form solution along either isopachic-fringes or along circles and radii from the crack-tip have been given, indicating in some cases large discrepancies between exact and approximate solutions.