Exact simulation of final, minimal and maximal values of Brownian motion and jump-diffusions with applications to option pricing

We introduce a method for generating (W_x,T^(m,s),m_x,T^(m,s),M_x,T^(m,s))(W_{x,T}^{(\mu,\sigma)},m_{x,T}^{(\mu,\sigma)},M_{x,T}^{(\mu,\sigma)}) , where W_x,T^(m,s)W_{x,T}^{(\mu,\sigma)} denotes the final value of a Brownian motion starting in x with drift μ and volatility σ at some final time T, m_x,T^(m,s) = inf_{0 ￡ t ￡ T}W_x,t^(m,s)m_{x,T}^{(\mu,\sigma)} = {\rm inf}_{0\leq t \leq T}W_{x,t}^{(\mu,\sigma)} and M_x,T^(m,s) = sup_{0 ￡ t ￡ T} W_x,t^(m,s)M_{x,T}^{(\mu,\sigma)} = {\rm sup}_{0\leq t \leq T} W_{x,t}^{(\mu,\sigma)} . By using the trivariate distribution of (W_x,T^(m,s),m_x,T^(m,s),M_x,T^(m,s))(W_{x,T}^{(\mu,\sigma)},m_{x,T}^{(\mu,\sigma)},M_{x,T}^{(\mu,\sigma)}) , we obtain a fast method which is unaffected by the well-known random walk approximation errors. The method is extended to jump-diffusion models. As sample applications we include Monte Carlo pricing methods for European double barrier knock-out calls with continuous reset conditions under both models. The proposed methods feature simple importance sampling techniques for variance reduction.     

Existence and uniqueness of periodic solutions for a class of nonlinear n ‐th order differential equations with delays

In this paper, we use the coincidence degree theory to establish new results on the existence and uniqueness of T ‐periodic solutions for a class of nonlinear n ‐th order differential equations with delays of the form x⁽ⁿ⁾(t) + f (x^{(n‐ 1)}(t)) + g (t, x (t ‐ τ (t))) = p (t).     

Global existence and time decay of small solutions to the landau-ginzburg type equations

We study the Cauchy problem for the nonlinear dissipative equations (0.1) uo∂u-αδu + Β|u|^2/n u = 0,x ∃ Rⁿ,t } 0,u(0,x) = u0(x),x ∃ Rⁿ, where α,Β ∃ C, ℜα 0. We are interested in the dissipative case ℜα 0, and ℜδ(α,Β)≥ 0, θ = |∫ u0(x)dx| ⊋ 0, where δ(α, Β) = ##|α|^n-1n^n/2 / ((n + 1)|α|² + α^{2 n/2}. Furthermore, we assume that the initial data u₀ ∃ L^p are such that (1 + |x|)^αu0 ∃ L¹, with sufficiently small norm ∃ = (1 + |x|)^α u₀ 1 + u₀ p, wherep 1, α ∃ (0,1). Then there exists a unique solution of the Cauchy problem (0.1)u(t, x) ∃ C ((0, ∞); L^∞) ∩ C ([0, ∞); L¹ ∩ L^p) satisfying the time decay estimates for allt0 u(t)||_∞ Cɛt^-n/2(1 + η log 〈t〉)^-n/2, if hg = θ^2/n 2_π ℜδ(α, Β) 0; u(t)||_∞ Cɛt^-n/2(1 + Μ log 〈t〉)^-n/4, if η = 0 and Μ = θ^4/n 4π)² (ℑδ(α, Β))² ℜ((1 + 1/n) υ1-1 υ2) 0; and u(t)||_∞ Cɛt^-n/2(1 + κ log 〈t〉)^-n/6, if η = 0, Μ = 0, κ 0, where υl,l = 1,2 are defined in (1.2), κ is a positive constant defined in (2.31).     

Topological entropy and blocking cost for geodesics in Riemannian manifolds

For a pair of points x, y in a compact, Riemannian manifold M let n _t (x, y) (resp. s _t (x, y)) be the number of geodesic segments with length ≤ t joining these points (resp. the minimal number of point obstacles needed to block these geodesic segments). We study relationships between the growth rates of n _t (x, y) and s _t (x, y) as t → ∞. We obtain lower bounds on s _t (x, y) in terms of the topological entropy h(M) and the fundamental group π ₁(M). For instance, we show that if h(M) > 0 then s _t grows exponentially, with the rate at least h(M)/2. This strengthens earlier results on blocking of geodesics (Burns and Gutkin Discrete Contin Dyn Syst 21:403–413, 2008; Lafont and Schmidt Geom Topol 11:867–887, 2007), and puts them in a new perspective.      

Identification of the unknown diffusion coefficient in a quasi‐linear parabolic equation by semigroup approach with mixed boundary conditions

In this article, a semigroup approach is presented for the mathematical analysis of the inverse coefficient problems of identifying the unknown diffusion coefficient k(u(x, t)) in the quasi‐linear parabolic equation u_t(x, t)=(k(u(x, t))u_x(x, t))_x, with Dirichlet boundary conditions u_x(0, t)=ψ₀, u(1, t)=ψ₁. The main purpose of this work is to analyze the distinguishability of the input–output mappings Φ[·] : ??→C¹[0, T], Ψ[·] : ??→C¹[0, T] using semigroup theory. In this article, it is shown that if the null space of semigroups T(t) and S(t) consists of only a zero function, then the input–output mappings Φ[·] and Ψ[·] have the distinguishability property. Copyright © 2008 John Wiley & Sons, Ltd.     

On the approximate behavior of the posterior distribution for an extreme multivariate observation

The behavior of the posterior for a large observation is considered. Two basic situations are discussed; location vectors and natural parameters.Let X = (X₁, X₂, …, X_n) be an observation from a multivariate exponential distribution with that natural parameter Θ = (Θ₁, Θ₂, …, Θ_n). Let θ_x^* be the posterior mode. Sufficient conditions are presented for the distribution of Θ − θ_x^* given X = x to converge to a multivariate normal with mean vector 0 as |x| tends to infinity. These same conditions imply that E(Θ | X = x) − θ_x^* converges to the zero vector as |x| tends to infinity.The posterior for an observation X = (X₁, X₂, …, X_n is considered for a location vector Θ = (Θ₁, Θ₂, …, Θ_n) as x gets large along a path, γ, in Rⁿ. Sufficient conditions are given for the distribution of γ(t) − Θ given X = γ(t) to converge in law as t → ∞. Slightly stronger conditions ensure that γ(t) − E(Θ | X = γ(t)) converges to the mean of the limiting distribution.These basic results about the posterior mean are extended to cover other estimators. Loss functions which are convex functions of absolute error are considered. Let δ be a Bayes estimator for a loss function of this type. Generally, if the distribution of Θ − E(Θ | X = γ(t)) given X = γ(t) converges in law to a symmetric distribution as t → ∞, it is shown that δ(γ(t)) − E(Θ | X = γ(t)) → 0 as t → ∞.     

Exact Rates of Convergence of Functional Limit Theorems for Csörgő;-Révész Increments of a Wiener Process

Let {W(t); t≥ 0} be a standard Wiener process and S be the Strassen set of functions. We investigate the exact rates of convergence to zero (as T→∞) of the variables $ \sup _{{0 \leqslant t \leqslant T - \alpha _{T} }} \inf _{{f \in S}} \sup _{{0 \leqslant x \leqslant 1}} {\left| {Y_{{t,T}} {\left( x \right)} - f{\left( x \right)}} \right|} Let {W(t); t≥ 0} be a standard Wiener process and S be the Strassen set of functions. We investigate the exact rates of convergence to zero (as T→∞) of the variables sup_{0≤ t ≤ T − aT} inf _f∈S sup_{0≤ x ≤1}|Y _t,T (x) −f(x)| and inf_{0≤ t ≤ T−aT} sup_{0≤ x ≤1}|Y _t,T (x−f(x)| for any given f∈S, where Y _t,T (x) = (W(t+xa _T ) −W(t)) (2a _T (log Ta _T ⁻¹ + log log T))^−1/2. We establish a relation between how small the increments are and the functional limit results of Cs?rg{\H o}-Révész increments for a Wiener process. Similar results for partial sums of i.i.d. random variables are also given. Received September 10, 1999, Accepted June 1, 2000     

A best approximation for the solution of one‐dimensional variable‐coefficient Burgers' equation

In this article, an iterative method for the approximate solution to one‐dimensional variable‐coefficient Burgers' equation is proposed in the reproducing kernel space W_(3,2). It is proved that the approximation w_n(x,t) converges to the exact solution u(x,t) for any initial function w₀(x,t) ε W_(3,2), and the approximate solution is the best approximation under a complete normal orthogonal system . Moreover the derivatives of w_n(x,t) are also uniformly convergent to the derivatives of u(x,t).© 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Semigroup approach for identification of the unknown diffusion coefficient in a quasi‐linear parabolic equation

This article presents a semigroup approach for the mathematical analysis of the inverse coefficient problems of identifying the unknown coefficient k(u(x,t)) in the quasi‐linear parabolic equation u_t(x,t)=(k(u(x,t))u_x(x,t))_x, with Dirichlet boundary conditions u(0,t)=ψ₀, u(1,t)=ψ₁. The main purpose of this paper is to investigate the distinguishability of the input–output mappings Φ[?]:?? →C¹[0,T], Ψ[?]:??→C¹[0,T] via semigroup theory. In this paper, it is shown that if the null space of the semigroup T(t) consists of only zero function, then the input–output mappings Φ[?] and Ψ[?] have the distinguishability property. It is also shown that the types of the boundary conditions and the region on which the problem is defined play an important role in the distinguishability property of these mappings. Moreover, under the light of measured output data (boundary observations) f(t):=k(u(0,t))u_x(0,t) or/and h(t):=k(u(1,t))u_x(1,t), the values k(ψ₀) and k(ψ₁) of the unknown diffusion coefficient k(u(x,t)) at (x,t)=(0,0) and (x,t)=(1,0), respectively, can be determined explicitly. In addition to these, the values k_u(ψ₀) and k_u(ψ₁) of the unknown coefficient k(u(x,t)) at (x,t)=(0,0) and (x,t)=(1,0), respectively, are also determined via the input data. Furthermore, it is shown that measured output data f(t) and h(t) can be determined analytically by an integral representation. Hence the input–output mappings Φ[?]:??→ C¹[0,T], Ψ[?]:??→C¹[0,T] are given explicitly in terms of the semigroup. Copyright © 2007 John Wiley & Sons, Ltd.     

On a Problem of Erd?s Involving the Largest Prime Factor of n

Let P(n) denote the largest prime factor of an integer n (N(x) = x (2+O(?{log₂ x/logx} ) )ò₂^xr(logx/logt) [(logt)/(t²)] d t,N(x) = x \left(2+O\left(\sqrt{\log_{2}\,x/\!\log x}\,\right) \right)\int_2^x\rho(\log x/\!\log t) {\log t\over t^2} {\rm d} t,     

On the existence and stability of periodic orbits in non ideal problems: General results

In this work, motivated by non-ideal mechanical systems, we investigate the following O.D.E. [(x)\dot] = f (x) + eg (x, t) + e²[^(g)] (x, t, e){\dot{x} = f (x) + \varepsilon g (x, t) + \varepsilon^{2}\widehat{g} (x, t, \varepsilon)} , where x ? W ì \mathbbRⁿ{x \in \Omega \subset \mathbb{R}^n} , g,[^(g)]{g,\widehat{g}} are T periodic functions of t and there is a ₀ ∈ Ω such that f ( a ₀) = 0 and f ′( a ₀) is a nilpotent matrix. When n = 3 and f (x) = (0, q (x ₃) , 0) we get results on existence and stability of periodic orbits. We apply these results in a non ideal mechanical system: the Centrifugal Vibrator. We make a stability analysis of this dynamical system and get a characterization of the Sommerfeld Effect as a bifurcation of periodic orbits.     

Une Remarque sur une Formule Sommatoire liée à la Somme des Chiffres

Let g ≥ 2 be an integer, and let s(n) be the sum of the digits of n in basis g. Let f(n) be a complex valued function defined on positive integers, such that ?_{n ￡ x} f(n)=o(x)\sum_{n\le x} f(n)=o(x) . We propose sufficient conditions on the function f to deduce the equality ?_{n ￡ x} f(s(n))=o(x)\sum_{n\le x} f(s(n))=o(x) . Applications are given, for instance, on the equidistribution mod 1 of the sequence (s(n))^α, where α is a positive real number.     

The scaling window of the 2‐SAT transition

We consider the random 2‐satisfiability (2‐SAT) problem, in which each instance is a formula that is the conjunction of m clauses of the form x∨y, chosen uniformly at random from among all 2‐clauses on n Boolean variables and their negations. As m and n tend to infinity in the ratio m/n→α, the problem is known to have a phase transition at α_c=1, below which the probability that the formula is satisfiable tends to one and above which it tends to zero. We determine the finite‐size scaling about this transition, namely the scaling of the maximal window W(n, δ)=(α_?(n,δ), α₊(n,δ)) such that the probability of satisfiability is greater than 1?δ for α<α_? and is less than δ for α>α₊. We show that W(n,δ)=(1?Θ(n^?1/3), 1+Θ(n^?1/3)), where the constants implicit in Θ depend on δ. We also determine the rates at which the probability of satisfiability approaches one and zero at the boundaries of the window. Namely, for m=(1+ε)n, where ε may depend on n as long as |ε| is sufficiently small and |ε|n^1/3 is sufficiently large, we show that the probability of satisfiability decays like exp(?Θ(nε³)) above the window, and goes to one like 1?Θ(n^?1|ε|^?3 below the window. We prove these results by defining an order parameter for the transition and establishing its scaling behavior in n both inside and outside the window. Using this order parameter, we prove that the 2‐SAT phase transition is continuous with an order parameter critical exponent of 1. We also determine the values of two other critical exponents, showing that the exponents of 2‐SAT are identical to those of the random graph. © 2001 John Wiley & Sons, Inc. Random Struct. Alg., 18: 201–256 2001     

Sur l'approximation de fonctions intégrables sur [0, 1] par des polynômes de Bernstein modifies

We study here a new kind of modified Bernstein polynomial operators on L¹(0, 1) introduced by J. L. Durrmeyer in [4]. We define for f integrable on [0, 1] the modified Bernstein polynomial M_n f: M_nf(x) = (n + 1) ∑ⁿ_{k = o}P_nk(x)∝¹₀ P_nk(t) f(t) dt. If the derivative d^r f/dx^r with r 0 is continuous on [0, 1], d^r/dx^rM_n f converge uniformly on [0,1] and sup_xε[0,1] ¦M_n f(x) − f(x)¦ 2ω_f(1/trn) if ω_f is the modulus of continuity of f. If f is in Sobolev space W_l,p(0, 1) with l 0, p 1, M_n f converge to f in w^l,p(0, 1).     

Multiparameter acoustic imaging in the born approximation

Samples of biological tissue are modelled as inhomogeneous fluids with density ?(X) and sound speed c(x) at point x. The samples are contained in the sphere |x| ? δ and it is assumed that ?(x) ? ?₀ = 1 and c(x) ? c₀ = 1 for |x| ? δ, and |γ_n(x)| ? 1 and |?γ?(x)| ? 1 where γ?(x) = ?(x) ? 1 and γ_n(x) = c^?2(x) ? 1. The samples are insonified by plane pulses s(x · θ₀ – t) where x = |θ₀| = 1 and the scattered pulse is shown to have the form |x|^?1 e_s(|x| – t, θ, θ₀) in the far field, where x = |x| θ. The response e_s(τ, θ, θ₀) is measurable. The goal of the work is to construct the sample parameters γn and γ? from e_s(τ, θ, θ₀) for suitable choiches of s, θ and θ₀. In the limiting case of constant density: γ?(x)_? 0 it is shown that Where δ represents the Dirac δ and S² is the unit sphere |θ| = 1. Analogous formulas, based on two sets of measurements, are derived for the case of variable c(x) and ?(x).     

An Application of a Mountain Pass Theorem

We are concerned with the following Dirichlet problem: −Δu(x) = f(x, u), x∈Ω, u∈H ¹ ₀(Ω), (P) where f(x, t) ∈C (×ℝ), f(x, t)/t is nondecreasing in t∈ℝ and tends to an L ^∞-function q(x) uniformly in x∈Ω as t→ + ∞ (i.e., f(x, t) is asymptotically linear in t at infinity). In this case, an Ambrosetti-Rabinowitz-type condition, that is, for some θ > 2, M > 0, 0 > θF(x, s) ≤f(x, s)s, for all |s|≥M and x∈Ω, (AR) is no longer true, where F(x, s) = ∫ ^s ₀ f(x, t)dt. As is well known, (AR) is an important technical condition in applying Mountain Pass Theorem. In this paper, without assuming (AR) we prove, by using a variant version of Mountain Pass Theorem, that problem (P) has a positive solution under suitable conditions on f(x, t) and q(x). Our methods also work for the case where f(x, t) is superlinear in t at infinity, i.e., q(x) ≡ +∞. Received June 24, 1998, Accepted January 14, 2000.     

On a nonlinear stochastic integral equation with application to control systems

The aim of the present paper is to study a nonlinear stochastic integral equation of the form

Representation and approximation of solutions to semilinear Volterra equations with delay

In this paper we use a theorem of Crandall and Pazy to provide the product integral representation of the nonlinear evolution operator associated with solutions to the semilinear Volterra equation: x(?)(t) = W(t, τ) ?(0) + ∝_τ^tW(t, s)F(s, x_s(?)) ds.Here the kernel W(t, s) is a linear evolution operator on a Banach space X; I is an interval of the form [?r, 0] or (?∞, 0] and F is a nonlinear mapping of R × C(I, X) into X. The abstract theory is applied to examples of partial functional differential equations.     

Spectral Criteria of Abstract Sequences; Application to Difference Equations

We suppose that M is a closed subspace of l ^∞(J, X), the space of all bounded sequences {x(n)} _n?J ? X, where J ? {Z⁺,Z} and X is a complex Banach space. We define the M-spectrum σ_M (u) of a sequence u ? l ^∞(J,X). Certain conditions will be supposed on both M and σ_M (u) to insure the existence of u ? M. We prove that if u is ergodic, such that σ_M (u,) is at most countable and, for every λ ? σ_M (u), the sequence e^?iλnu(n) is ergodic, then u ? M. We apply this result to the operator difference equationu(n + 1) = Au(n) + ψ(n), n ? J,and to the infinite order difference equation Σ ^r _k=1 a_k (u(n + k) ? u(n)) + Σ _{s ? Z}?(n ? s)u(s) = h(n), n?J, where ψ?l ^∞(Z,X) such that ψ| _J ? M, A is the generator of a C ₀-semigroup of linear bounded operators {T(t)} _t>0 on X, h ? M, ? ? l ¹(Z) and a_k ?C. Certain conditions will be imposed to guarantee the existence of solutions in the class M.