Multiquadratic methods in the numerical. treatment of the integro-differential equations of collective non-ruin; the finite time case
|
|
- Briana Chapman
- 5 years ago
- Views:
Transcription
1 Multiquadratic methods in the numerical treatment of the integro-differential equations of collective non-ruin; the finite time case Athena Makroglou Division of Mathematics and Statistics, School of Computer Science and Mathematics, University of Portsmouth, Mercantile House, Portsmouth PO7 2EG, EngWd alpha2. iso.port.ac. uk Abstract An important problem of collective non-ruin is the estimation of the probabilities R(z,t) and R(z) of the finite and ultimate non-ruin respectively, where t is time and z is the initial reserve. The governing equations are first order Volterra integro-differential equations (VIDEs). Computational results for the finite time case are presented. These are obtained by using multiquadric approximation methods applied directly to the VIDE. Comparisons are made with results obtained by indirect approaches. 1 Introduction The problem of estimating the probability of non-ruin given an initial reserve is very important in modelling a risk business such as an insurance company. In the case offinitetime non-ruin, we are concerned with the estimation of the probability R(z,t) = P[Z(a)>Q,Q<8<t\Z(Q)=z], (1.1) where, using the notation of Peters and Mangel (1996)), Z(s) denotes the risk reserve at time s and z = Z(0) is the initial risk reserve.
2 496 Boundary Element Technology The simple model of the risk business of an insurance company considered, assumes that: Risk reserve = initial reserve + total premiums - total claims. In the absence of claims, it is assumed that the reserve satisfies the deterministic equation ^ = & + TZ«, (1.2) where /% denotes the rate (constant) at which the premiums come in, and 7 is the interest rate that applies to risk reserve until a claim is made. Claim sizes at time t denoted by Xt are assumed to have distribution function B(x) P[Xt < x] and are assumed to arrive according to a Poisson process Nt with parameter A. This means that Then St = ]Ci=i ^, the accumulated claims process is a compound Poisson process. Models derived under more complicated assumptions can be found for example in Seal (1974), Bade (1983). Using the above assumptions, an integro-differential equation model is derived for the probability R(z,t), see for examples Arfwedson (1950), Knessl and Peters (1994, 1996) and Grandell (1991) for more details. This problem has been solved by a number of theoretical and numerical approaches involving use of Laplace transforms (cf. Seal (1974), Knessl and Peters ( 1994, 1996), Makroglou (1998)). Solving the VIDEs numerically is very important since it allows treatment of the problem under more complicated assumptions for B(x) and /3(Z(t)) when analytical solutions are difficult to obtain. In this paper we are concerned with the computational treatment of the problem following the approach of solving directly the VIDE in R(z, t) numerically using multiquadric (MQ) methods. Solving numerically the VIDEs directly is also very important, since it is a more efficient and simpler way of approach than the indirect one based on the use of Laplace transforms. The organization of the paper is as follows: In section 2 the form of the VIDEs obeyed by R(z, t) is given together with an overview of existing methods of solution. In section 3, the method of solving directly the VIDE by using multiquadric (MQ) methods is described. In section 4, numerical results are presented. Conclusions and some comparisons with the results of Makroglou (1998) can be found in section 5.
3 Boundary Element Technology The integro-differential equations models Starting from equation (1.1) and using probabilistic arguments (cf. Knessl and Peters (1994)) the following equation is obtained for the probability of non-ruin R(z,t): with conditions R(z, t) = 1, R(z, 0) = l,z > 0. (2.2) We also have that R(z, t) = 0, z < 0. Most of the methods applied to the solution of (2.1)-(2.2) treat special cases of the problem concerning the choice of the model for the premiums ((3(Z(t))) and the distribution function B(x). The choices for the premiums are: and 0(Z(t)) =0 = constant (2.3) 0(Z(t)) = A> + 7^(*),A),7constants. (2.4) The choices for the distribution function B(x) and its density function b(x) include the Gamma distribution with v 1 (exponential) and v = 2, i.e., where {W-'ae-" jf _>n ' "W 'f*-0 (2.5) u otherwise (i/)= Jo is the Gamma function. Existing methods for the solution of equations (2.1)-(2.2) include Seal (1974, 1978)), Knessl and Peters (1994, 1996), and Makroglou (1998). The approach of Seal (1974, 1978) in brief, is to consider the Laplace transform of R(z,t) over z values, i.e. %*(s) = f e~'*r(z,t)dz and then obtain a partial differential equation for vt(s) which is solved analytically in terms of #(0, t) and w(s) = f e~'*b(y)dy. To find R(z, t), vt(s) is inverted to obtain a formula which involves the inverse Laplace transforms of three more functions which are obtained by numerical methods. Numerical integration is also used for the evaluation of an integral occuring in the formula
4 498 Boundary Element Technology for R(z,t). The approach of Knessl and Peters (1994) can be summarized as follows: They consider the Laplace transform Q(z,s) of the function R(z,t) over time, i.e. they take e~**r(z, t)dt. (2.6) For the special case that B(x) is the exponential distribution (B'(x) = b(x) ae'"*, p, = I/a, the mean) this leads to a linear integro-differential equation in Q(z, s) which they solve analytically by conversion to a second order ordinary differential equation. Analytical inversion of Q(z,s) is then applied to obtain R(z,t). Both cases (2.3) and (2.4) for the premiums are considered. Exact and asymptotic solutions are derived for case (2.3) and exact ones for the case (2.4). Asymptotic solutions for the case (2.4) were given in Knessl and Peters (1996). Some of the formulae are complicated and involve series and integrals. Some numerical results which are obtained by numerical evaluation of the final formulae for R(z,t) using MATHEMATICA can be found in Knessl and Peters (1996). Some of the results of Seal (1974) and Knessl and Peters (1994) were verified in Harper (1995, 1996) where the numerical inversion of the analytical expression for vt(s) and Q(z,s) respectively, was done using ACM routine 619 (Piessens and Huysmans (1984)). All the above approaches depend on the ability to be able to solve analytically the equation for the Laplace transforms of the unknown function R(z,t). For more general choices of B(x) and /3(Z(t)) though this is not always possible and even for the simple choices of B(x) and 0(Z(t)), considered the form of the resulting formulae is often complicated. So in Makroglou (1998), a fully numerical method was implemented using polynomial collocation methods to solve the integro-differential equation in Q(z, s) and then numerical inversion to obtain values of R(z,t). The analytical (exact) solution provided in Knessl and Peters (1994) was used for testing purposes and computation of the true errors. In this paper a method for the direct numerical solution of equations (2.1)-(2.2) is presented which does not use Laplace transforms. In particular, multiquadric approximation methods are used (see for example Powell (1987, 1992) for an introduction to the subject) extending the line of approach of Bonzani (1997) as applied to PDEs. The description of the method is given in next section.
5 Boundary Element Technology Description of the multiquadric method for solving (2.1)-(2.2) We are going to use multiquadric approximations (cf. Kansa (1990a,b), Makroglou (1994)) for the numerical solution of (2.1)-(2.2) following the approach of Bonzani (1997) who applied sine and Lagrange type approximations to second order semilinear partial differential equations in one space dimension. For fixed t, the following approximation (see also Kansa (1990a,b), to the unknown function R(z, t) of equation (2.1) is used: n+l (t)fi,(z), (3,1) where and where 0 = z\ <... < Zn < Zn+i = 2, * > 0, and c is the multiquadric parameter. The boundary condition g(t) = R(z, t) was found using the true solution. This boundary condition was used to express dn+i(t) in (3.1) in terms of a\ ( ),..., a,n(t). In practice, a z value has to be chosen to represent the oo, so that the given boundary condition lim^oo R(z^ t) = I can be used instead. The coefficients a,j(t) are then to be estimated so that (3.1)-(3.3) satisfies the equation (2.1)-(2.2), giving afirstorder system of ordinary differential equations in aj (t) of the form where 5, F are n x n matrices with elements s-'fa + S-'M*), (3.4) Skj = Pj(zk), 771=1 kj = l,2,...,n. (3.5)
6 500 Boundary Element Technology and b(t) is given by where k + XhY ' -^ m=l & = l,2,...,n (3.6) (3.7) To obtain equations (3.5) we also used a quadrature rule with weights Wkm and stepsize h for the evaluation of the integrals J** Pj(zk - y)b(y)dy and To solve the ODE system (3.4) we also need starting values a(0). These are found using the second of the conditions (2.2) for B, that is R(z, 0) = l,z > 0 and (3.1). We thus obtain the equations (3.8) After we solve the ODE system (3.4)-(3.8) for the values of t we need, the function R(z, t) is evaluated from (3.1)-(3.3). For other methods applied to hyperbolic partial differential equations of higher order which include error estimates (but no numerical results) we refer for example to Yanik and Fairweather (1988), Pani, Thomee and Wahlbin (1992), Fairweather (1994). 4 Numerical results Numerical results are presented for the solution R(z, t),z 1.0,0 < t < 10 of (2.1)-(2.2) using with a = 1.0, A) = 2.0, A = 7 = 1.0. The true solution is given (Knessl and Peters (1994)) as, (4.1)
7 Boundary Element Technology 501 (4.2) The results of method COL-LAP in Table 1 (Makroglou (1998)) were obtained using polynomial collocation methods with h = 0.1 and collocation parameters 0o = 0,0i = 1, which makes the method identical to that resulting from the use of the trapezoidal method, in combination with the ACM routine 619 (Piessens (1984)) with tolerance equal to 1.D - 4 for the inversions of the Laplace transforms. The results of the MQ method in Table 1 were obtained using the direct multiquadric method presented in this paper. Table 4.1: Results for R(l,t) t exact R(l,t) COL - LAP true error 0.74D D D D - 4 MQ,c = true error 0.25D J D-2
8 502 Boundary Element Technology The results for the direct MQ method were obtained using Matlab 5.0 with the ODE23 rourine for the solution of the system of differential equations for obtaining the a?(f),j = 1,...,n. The MQ approximations used n = 10. The integrals were evaluated using the trapezoidal quadrature rule with mesh points Z{, i 1,..., k and stepsize h 0.1. Prom Table 1 we may see that the results obtained by the MQ method and the collocation method of Makroglou (1998) are in good agreement with the exact solution of Knessl (1994) with these obtained by the collocation method considerably smaller. The accuracy of the MQ method might improve if the integrals were evaluated to convergence. That way, one could easily change the value of n and get results for different MQ approximations easily. Also, using different values of the c parameter for different t values or certain methods of estimation (cf. Golberg, Chen and Karur (1996)) might benefit the MQ method. Tests with different ODE solvers is also worth considering. 5 Conclusions The estimation of the probability offinitetime non - ruin R(z, t) was obtained by solving directly the governing integro-differential equation by using multiquadric approximations. The direct treatment of the problem described here and also the fully numerical indirect approach of Makroglou (1998), allow the easy extension for other types of distribution function B(x) and premiums function 0(z) whenfindingthe exact solutions for R or the corresponding Laplace transforms of R is not easy or possible, and it is hoped to be useful for the actuarial society. Global multiquadric approximations were used with respect to the "space" variable z. Next line of approach to follow would be the use of multiquadric approximations with compact support (cf. Wendland (1995, 1998)) or to extend other methods which are in use for the numerical solution of first order hyperbolic partial differential equations (see for example Morton and Mayers (1994, chapter 4) for an introduction). Keywords: Partial integro-differential equations,firstorder, numerical solution, multiquadrics, actuarial risk management. References [1] Arfwedson, G., Some problems in the collective theory of risk, Skand. Aktuar. Tidskr, 33 (1950), 1-38.
9 Boundary Element Technology 503 [2] Bonzani, I., Solution of nonlinear evolution problems by parallelized collocation-interpolation methods, Computers Math. Applic., 34 (1997), [3] Bade, J.P., The ruin problem for mixed Poisson risk processes, Scand. Actuarial J., (1983), [4] Fairweather, G., Spline collocation methods for a class of hyperbolic partial integro-differential equations, SIAM J. Numer. Anal., 31 (1994), [5] Golberg, M.A., Chen, C.S., Karur, S.R., Improved multiquadric approximation for partial differential equations, Engng Anal, with Boundary Elements, 1996, [6] Grandell, J., Aspects of risk theory, Springer-Verlag, [7] Harper, W., Numerical treatment of the problem of collective ruin, Final year project, School of Computer Science and Mathematics, Div. of Mathematics and Statistics, University of Portsmouth, UK, [8] Harper, W., On the estimation of the survival probability for the problem offinitetime collective non-ruin, M.Sc. thesis, School of Computer Science and Mathematics, Div. of Mathematics and Statistics, University of Portsmouth, UK, [9] Kansa, E.J., Multiquadrics - A scattered data approximation scheme with applications to computational fluid dynamics-i, Surface approximations and partial derivatives estimates, Computers Math. Applic., 19 (1990a), [10] Kansa, E.J., Multiquadrics - A scattered data approximation scheme with applications to computational fluid dynamics-ii, Solutions to parabolic, hyperbolic and elliptic partial differential equations, Computers Math. Applic., 19 (1990b), [11] Knessl, C. and Peters, C.S., Exact and asymptotic solutions for the time dependent problem of collective ruin I, SIAM J. Appl. Math., 54 (1994), [12] Knessl, C. and Peters, C.S., Exact and asymptotic solutions for the time-dependent problem of collective ruin II, SIAM. J. Appl. Math., 56 (1996), [13] Makroglou, A., Radial basis functions in the numerical solution of nonlinear Volterra integral equations, J. Applied Science and Computations, 1 (1994),
10 504 Boundary Element Technology [14] Makroglou, A., Computer treatment of the problem of collective nonruin; the finite time case, contributed talk: 18th International Congress on Computational and Applied Mathematics (ICCAM '98), July 27 - August 1, 1998, Katholieke Universiteit Leuven, Belgium, and Conference HERCMA '98, September, 1998, AUEB, Athens, Greece. [15] Morton, K.W. and Mayers, D.F., Numerical solution of partial differential equations, Cambridge University Press, [16] Pani, A.K., Thomee, V. and Wahlbin, L.B., Numerical methods for hyperbolic and parabolic integro-differential equations, J. Integral Equations, 4 (1992), [17] Peters, C.S. and Mangel, M., New methods for the problem of collective ruin, SIAM J. Appl. Math., 50 (1990), [18] Piessens, R., Alg. 619, Automatic numerical inversion of the Laplace transform, ACM Trans. Math. Soft., 10 (1984), [19] Powell, M.J.D., Radial basis functions for multivariable interpolation: a review, pp in: Algorithms for Approximation, Mason, J.C. and Cox, M.G. (Eds), Clarendon Press, [20] Powell, M.J.D., The theory of radial basis functions in 1990, pp in: Advances in Numerical Analysis II: Wavelets, Subdivision, and Radial Basis Functions, W. Light (Ed.), Oxford University Press, [21] Seal, H.L., The numerical calculation of 7(w;,t), the probability of non-ruin in an interval (0,t), Scand. Actuarial J., (1974), [22] Seal, H.L., Survival probabilities. The goal of risk theory, John Wiley and Sons, [23] Wendland, H., Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree, Advances in Computational Mathematics, 4 (1995), [24] Wendland, H., Error estimates for interpolation by compactly supported radial basis functions of minimal degree, J. Approx. Theory, 93 (1998), [25] Yanik, E.G. and Fair weather, G., Finite element methods for parabolic and hyperbolic partial integro-differential equations, Nonlinear Analysis, Theory, Methods and Applications, 12 (1988), Acknowledgement. The author of this article would like to thank Professor C.S. Chen of the University of Nevada, for inviting her to the Conference BETECH 99 and for providing her with additional references on the theory and applications of multiquadric approximations.
Numerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method
Numerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method Y. Azari Keywords: Local RBF-based finite difference (LRBF-FD), Global RBF collocation, sine-gordon
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
More informationDiscrete Projection Methods for Integral Equations
SUB Gttttingen 7 208 427 244 98 A 5141 Discrete Projection Methods for Integral Equations M.A. Golberg & C.S. Chen TM Computational Mechanics Publications Southampton UK and Boston USA Contents Sources
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 5: Completely Monotone and Multiply Monotone Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 9: Conditionally Positive Definite Radial Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH
More informationA Point Interpolation Meshless Method for the. Numerical Solution of the Singularly Perturbed. Integral and Integro-differential Equations
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 13, 643-656 HIKARI Ltd, www.m-hikari.com A Point Interpolation Meshless Method for the Numerical Solution of the Singularly Perturbed Integral and Integro-differential
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 39: Non-Symmetric RBF Collocation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 40: Symmetric RBF Collocation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationMeshfree Approximation Methods with MATLAB
Interdisciplinary Mathematical Sc Meshfree Approximation Methods with MATLAB Gregory E. Fasshauer Illinois Institute of Technology, USA Y f? World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI
More information1 Comparison of Kansa s method versus Method of Fundamental Solution (MFS) and Dual Reciprocity Method of Fundamental Solution (MFS- DRM)
1 Comparison of Kansa s method versus Method of Fundamental Solution (MFS) and Dual Reciprocity Method of Fundamental Solution (MFS- DRM) 1.1 Introduction In this work, performances of two most widely
More informationA Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions
Applied Mathematical Sciences, Vol. 5, 2011, no. 23, 1145-1152 A Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions Z. Avazzadeh
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationIn the present work, we consider the singularly perturbed Volterra integro-differential equations (SVIDE) x K(x, t, ε, y(t)) dt, x I = [0, X], (1)
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.12(211) No.4,pp.43-441 An Approximation Algorithm for the Solution of the Singularly Perturbed Volterra Integro-differential
More informationNumerical Methods for Engineers and Scientists
Numerical Methods for Engineers and Scientists Second Edition Revised and Expanded Joe D. Hoffman Department of Mechanical Engineering Purdue University West Lafayette, Indiana m MARCEL D E К К E R MARCEL
More informationWinfried Auzinger, Ernst Karner, Othmar Koch, Ewa Weinmüller
Opuscula Mathematica Vol. 26 No. 2 26 Winfried Auzinger, Ernst Karner, Othmar Koch, Ewa Weinmüller COLLOCATION METHODS FOR THE SOLUTION OF EIGENVALUE PROBLEMS FOR SINGULAR ORDINARY DIFFERENTIAL EQUATIONS
More informationBIBLIOGRAPHY KENDALL E ATKINSON
September 28, 2017 BIBLIOGRAPHY KENDALL E ATKINSON I. BOOKS 1. K. Atkinson, D. Chien, and O. Hansen. Spectral Methods Using Multivariate Polynomials On The Unit Ball, submitted for publication. 2. K. Atkinson
More informationOn separated solutions of logistic population equation with harvesting
On separated solutions of logistic population equation with harvesting Philip Korman Department of Mathematical Sciences University of Cincinnati Cincinnati Ohio 4522-25 Abstract We provide a surprising
More informationAN ELEMENTARY PROOF OF THE OPTIMAL RECOVERY OF THE THIN PLATE SPLINE RADIAL BASIS FUNCTION
J. KSIAM Vol.19, No.4, 409 416, 2015 http://dx.doi.org/10.12941/jksiam.2015.19.409 AN ELEMENTARY PROOF OF THE OPTIMAL RECOVERY OF THE THIN PLATE SPLINE RADIAL BASIS FUNCTION MORAN KIM 1 AND CHOHONG MIN
More informationUniversity Of Calgary Department of Mathematics and Statistics
University Of Calgary Department of Mathematics and Statistics Hawkes Seminar May 23, 2018 1 / 46 Some Problems in Insurance and Reinsurance Mohamed Badaoui Department of Electrical Engineering National
More informationNUMERICAL SOLUTION FOR CLASS OF ONE DIMENSIONAL PARABOLIC PARTIAL INTEGRO DIFFERENTIAL EQUATIONS VIA LEGENDRE SPECTRAL COLLOCATION METHOD
Journal of Fractional Calculus and Applications, Vol. 5(3S) No., pp. 1-11. (6th. Symp. Frac. Calc. Appl. August, 01). ISSN: 090-5858. http://fcag-egypt.com/journals/jfca/ NUMERICAL SOLUTION FOR CLASS OF
More informationTHE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS. Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands
THE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands 1. Introduction This paper deals with initial value problems for delay
More informationFractional Spectral and Spectral Element Methods
Fractional Calculus, Probability and Non-local Operators: Applications and Recent Developments Nov. 6th - 8th 2013, BCAM, Bilbao, Spain Fractional Spectral and Spectral Element Methods (Based on PhD thesis
More informationPractical approaches to the estimation of the ruin probability in a risk model with additional funds
Modern Stochastics: Theory and Applications (204) 67 80 DOI: 05559/5-VMSTA8 Practical approaches to the estimation of the ruin probability in a risk model with additional funds Yuliya Mishura a Olena Ragulina
More informationSynchronized Queues with Deterministic Arrivals
Synchronized Queues with Deterministic Arrivals Dimitra Pinotsi and Michael A. Zazanis Department of Statistics Athens University of Economics and Business 76 Patission str., Athens 14 34, Greece Abstract
More informationA polynomial expansion to approximate ruin probabilities
A polynomial expansion to approximate ruin probabilities P.O. Goffard 1 X. Guerrault 2 S. Loisel 3 D. Pommerêt 4 1 Axa France - Institut de mathématiques de Luminy Université de Aix-Marseille 2 Axa France
More informationQuadrature Prefilters for the Discrete Wavelet Transform. Bruce R. Johnson. James L. Kinsey. Abstract
Quadrature Prefilters for the Discrete Wavelet Transform Bruce R. Johnson James L. Kinsey Abstract Discrepancies between the Discrete Wavelet Transform and the coefficients of the Wavelet Series are known
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More informationRBF Collocation Methods and Pseudospectral Methods
RBF Collocation Methods and Pseudospectral Methods G. E. Fasshauer Draft: November 19, 24 Abstract We show how the collocation framework that is prevalent in the radial basis function literature can be
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationBeyond Wiener Askey Expansions: Handling Arbitrary PDFs
Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2005) DOI: 10.1007/s10915-005-9038-8 Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Xiaoliang Wan 1 and George Em Karniadakis 1
More informationThe Application of the Poincart-Transform to the Lotka-Volterra Model
J. Math. Biology 6, 67-73 (1978) Journal of by Springer-Verlag 1978 The Application of the Poincart-Transform to the Lotka-Volterra Model S. B. Hsu Department of Mathematics, University of Utah, Salt Lake
More informationNUMERICAL METHOD FOR THE MIXED VOLTERRA-FREDHOLM INTEGRAL EQUATIONS USING HYBRID LEGENDRE FUNCTIONS
Conference Applications of Mathematics 215, in honor of the birthday anniversaries of Ivo Babuška (9), Milan Práger (85), and Emil Vitásek (85) J. Brandts, S. Korotov, M. Křížek, K. Segeth, J. Šístek,
More informationCubic B-spline Collocation Method for Fourth Order Boundary Value Problems. 1 Introduction
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.142012 No.3,pp.336-344 Cubic B-spline Collocation Method for Fourth Order Boundary Value Problems K.N.S. Kasi Viswanadham,
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationStability of the Defect Renewal Volterra Integral Equations
19th International Congress on Modelling and Simulation, Perth, Australia, 12 16 December 211 http://mssanz.org.au/modsim211 Stability of the Defect Renewal Volterra Integral Equations R. S. Anderssen,
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationFinite-time Ruin Probability of Renewal Model with Risky Investment and Subexponential Claims
Proceedings of the World Congress on Engineering 29 Vol II WCE 29, July 1-3, 29, London, U.K. Finite-time Ruin Probability of Renewal Model with Risky Investment and Subexponential Claims Tao Jiang Abstract
More informationIdentification of a memory kernel in a nonlinear parabolic integro-differential problem
FACULTY OF ENGINEERING AND ARCHITECTURE Identification of a memory kernel in a nonlinear parabolic integro-differential problem K. Van Bockstal, M. Slodička and F. Gistelinck Ghent University Department
More informationThe method of fundamental solutions for Poisson's equation M.A. Golberg Czrc/e
The method of fundamental solutions for Poisson's equation M.A. Golberg Czrc/e ABSTRACT We show how to extend the method of fundamental solutions (MFS) to solve Poisson's equation in 2D without boundary
More informationApplied Mathematics Letters. Nonlinear stability of discontinuous Galerkin methods for delay differential equations
Applied Mathematics Letters 23 21 457 461 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Nonlinear stability of discontinuous Galerkin
More informationNumerical solution of Maxwell equations using local weak form meshless techniques
Journal of mathematics and computer science 13 2014), 168-185 Numerical solution of Maxwell equations using local weak form meshless techniques S. Sarabadan 1, M. Shahrezaee 1, J.A. Rad 2, K. Parand 2,*
More informationCubic spline Numerov type approach for solution of Helmholtz equation
Journal of Linear and Topological Algebra Vol. 03, No. 01, 2014, 47-54 Cubic spline Numerov type approach for solution of Helmholtz equation J. Rashidinia a, H. S. Shekarabi a and M. Aghamohamadi a a Department
More informationSyllabus for Applied Mathematics Graduate Student Qualifying Exams, Dartmouth Mathematics Department
Syllabus for Applied Mathematics Graduate Student Qualifying Exams, Dartmouth Mathematics Department Alex Barnett, Scott Pauls, Dan Rockmore August 12, 2011 We aim to touch upon many topics that a professional
More informationThe finite-time Gerber-Shiu penalty function for two classes of risk processes
The finite-time Gerber-Shiu penalty function for two classes of risk processes July 10, 2014 49 th Actuarial Research Conference University of California, Santa Barbara, July 13 July 16, 2014 The finite
More informationOn the decay of elements of inverse triangular Toeplitz matrix
On the decay of elements of inverse triangular Toeplitz matrix Neville Ford, D. V. Savostyanov, N. L. Zamarashkin August 03, 203 arxiv:308.0724v [math.na] 3 Aug 203 Abstract We consider half-infinite triangular
More informationModeling Recurrent Events in Panel Data Using Mixed Poisson Models
Modeling Recurrent Events in Panel Data Using Mixed Poisson Models V. Savani and A. Zhigljavsky Abstract This paper reviews the applicability of the mixed Poisson process as a model for recurrent events
More informationExistence of solution and solving the integro-differential equations system by the multi-wavelet Petrov-Galerkin method
Int. J. Nonlinear Anal. Appl. 7 (6) No., 7-8 ISSN: 8-68 (electronic) http://dx.doi.org/.75/ijnaa.5.37 Existence of solution and solving the integro-differential equations system by the multi-wavelet Petrov-Galerkin
More informationA Nodal Spline Collocation Method for the Solution of Cauchy Singular Integral Equations 1
European Society of Computational Methods in Sciences and Engineering (ESCMSE) Journal of Numerical Analysis, Industrial and Applied Mathematics (JNAIAM) vol. 3, no. 3-4, 2008, pp. 211-220 ISSN 1790 8140
More informationQuarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations
MATEMATIKA, 2011, Volume 27, Number 2, 199 208 c Department of Mathematical Sciences, UTM Quarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations 1 E. Aruchunan
More informationIntroduction to Computational Stochastic Differential Equations
Introduction to Computational Stochastic Differential Equations Gabriel J. Lord Catherine E. Powell Tony Shardlow Preface Techniques for solving many of the differential equations traditionally used by
More informationIf you make reference to this version of the manuscript, use the following information:
This is the author s final, peer-reviewed manuscript as accepted for publication. The publisher-formatted version may be available through the publisher s web site or your institution s library. Inverse
More informationCAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS
CAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS Preliminaries Round-off errors and computer arithmetic, algorithms and convergence Solutions of Equations in One Variable Bisection method, fixed-point
More informationBiorthogonal Spline Type Wavelets
PERGAMON Computers and Mathematics with Applications 0 (00 1 0 www.elsevier.com/locate/camwa Biorthogonal Spline Type Wavelets Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan
More informationON NUMERICAL SOLUTIONS OF PARTIAL DIFFERENTIAL EQUATIONS BY THE DECOMPOSITION METHOD. Mustafa Inc
153 Kragujevac J. Math. 26 (2004) 153 164. ON NUMERICAL SOLUTIONS OF PARTIAL DIFFERENTIAL EQUATIONS BY THE DECOMPOSITION METHOD Mustafa Inc Department of Mathematics, Firat University, 23119 Elazig Turkiye
More informationCalculation of Cylindrical Functions using Correction of Asymptotic Expansions
Universal Journal of Applied Mathematics & Computation (4), 5-46 www.papersciences.com Calculation of Cylindrical Functions using Correction of Asymptotic Epansions G.B. Muravskii Faculty of Civil and
More informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationStudy the Numerical Methods for Solving System of Equation
Study the Numerical Methods for Solving System of Equation Ravi Kumar 1, Mr. Raj Kumar Duhan 2 1 M. Tech. (M.E), 4 th Semester, UIET MDU Rohtak 2 Assistant Professor, Dept. of Mechanical Engg., UIET MDU
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationAnalysis of the ruin probability using Laplace transforms and Karamata Tauberian theorems
Analysis of the ruin probability using Laplace transforms and Karamata Tauberian theorems Corina Constantinescu and Enrique Thomann Abstract The classical result of Cramer-Lundberg states that if the rate
More informationΨ-asymptotic stability of non-linear matrix Lyapunov systems
Available online at wwwtjnsacom J Nonlinear Sci Appl 5 (22), 5 25 Research Article Ψ-asymptotic stability of non-linear matrix Lyapunov systems MSNMurty a,, GSuresh Kumar b a Department of Applied Mathematics,
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 43: RBF-PS Methods in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 43 1 Outline
More informationSeyyed Ahmad Hosseini
Seyyed Ahmad Hosseini Curriculum Vitae Assistant Professor of Applied Mathematics-Numerical Analysis at Golestan University Gender Personal Details Male Date of birth 21 st March, 1982 Place of birth Present
More informationKatholieke Universiteit Leuven Department of Computer Science
Interpolation with quintic Powell-Sabin splines Hendrik Speleers Report TW 583, January 2011 Katholieke Universiteit Leuven Department of Computer Science Celestijnenlaan 200A B-3001 Heverlee (Belgium)
More informationNUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS
TWMS J Pure Appl Math V5, N2, 24, pp22-228 NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS S ASADI, AH BORZABADI Abstract In this paper, Haar wavelet benefits are applied to the delay
More informationCompact Local Stencils Employed With Integrated RBFs For Fourth-Order Differential Problems
Copyright 2011 Tech Science Press SL, vol.6, no.2, pp.93-107, 2011 Compact Local Stencils Employed With Integrated RBFs For Fourth-Order Differential Problems T.-T. Hoang-Trieu 1, N. Mai-Duy 1 and T. Tran-Cong
More informationBarycentric rational interpolation with no poles and high rates of approximation
Barycentric rational interpolation with no poles and high rates of approximation Michael S. Floater Kai Hormann Abstract It is well known that rational interpolation sometimes gives better approximations
More informationMultivariate Risk Processes with Interacting Intensities
Multivariate Risk Processes with Interacting Intensities Nicole Bäuerle (joint work with Rudolf Grübel) Luminy, April 2010 Outline Multivariate pure birth processes Multivariate Risk Processes Fluid Limits
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationA parallel algorithm for the heat equation with derivative boundary conditions
International Mathematical Forum, 2, 2007, no. 12, 565-574 A parallel algorithm for the heat equation with derivative boundary conditions M. Akram University College of Information Technology University
More informationUniversity of Hertfordshire Department of Mathematics. Study on the Dual Reciprocity Boundary Element Method
University of Hertfordshire Department of Mathematics Study on the Dual Reciprocity Boundary Element Method Wattana Toutip Technical Report 3 July 999 Preface The boundary Element method (BEM) is now recognised
More informationThe integrating factor method (Sect. 1.1)
The integrating factor method (Sect. 1.1) Overview of differential equations. Linear Ordinary Differential Equations. The integrating factor method. Constant coefficients. The Initial Value Problem. Overview
More informationWendland Functions A C++ code to compute them
Wendland Functions A C++ code to compute them Carlos Argáez 1, Sigurdur Hafstein 1 and Peter Giesl 2 1 Faculty of Physical Sciences, University of Iceland, 107 Reykjavík, Iceland 2 Department of Mathematics,
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationNinth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673
Ninth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673 A NEW NUMERICAL APPROACH FOR THE SOLUTION OF CONTAMINANT TRANSPORT EQUATION Mohamed El-Gamel Department of Mathematical
More informationAn Efficient Numerical Scheme for Coupled Nonlinear Burgers Equations
Appl. Math. Inf. Sci. 9, o. 1, 245-255 (2015) 245 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.12785/amis/090130 An Efficient umerical Scheme for Coupled onlinear
More informationSchemes. Philipp Keding AT15, San Antonio, TX. Quarklet Frames in Adaptive Numerical. Schemes. Philipp Keding. Philipps-University Marburg
AT15, San Antonio, TX 5-23-2016 Joint work with Stephan Dahlke and Thorsten Raasch Let H be a Hilbert space with its dual H. We are interested in the solution u H of the operator equation Lu = f, with
More informationA numerical solution of a Kawahara equation by using Multiquadric radial basis function
Mathematics Scientific Journal Vol. 9, No. 1, (013), 115-15 A numerical solution of a Kawahara equation by using Multiquadric radial basis function M. Zarebnia a, M. Takhti a a Department of Mathematics,
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More information[313 ] A USE OF COMPLEX PROBABILITIES IN THE THEORY OF STOCHASTIC PROCESSES
[313 ] A USE OF COMPLEX PROBABILITIES IN THE THEORY OF STOCHASTIC PROCESSES BY D. R. COX Received 17 September 1954 ABSTRACT. The exponential distribution is very important in the theory of stochastic
More informationIterated Defect Correction with B-Splines for a Class of Strongly Nonlinear Two-Point Boundary Value Problems
American Review of Mathematics and Statistics June 2016, Vol. 4, No. 1, pp. 31-44 ISSN: 2374-2348 (Print), 2374-2356 (Online) Copyright The Author(s). All Rights Reserved. Published by American Research
More informationProperties of BPFs for Approximating the Solution of Nonlinear Fredholm Integro Differential Equation
Applied Mathematical Sciences, Vol. 6, 212, no. 32, 1563-1569 Properties of BPFs for Approximating the Solution of Nonlinear Fredholm Integro Differential Equation Ahmad Shahsavaran 1 and Abar Shahsavaran
More informationA CLASS OF CONTINUOUS HYBRID LINEAR MULTISTEP METHODS FOR STIFF IVPs IN ODEs
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI (S.N.) MATEMATICĂ, Tomul LVIII, 0, f. A CLASS OF CONTINUOUS HYBRID LINEAR MULTISTEP METHODS FOR STIFF IVPs IN ODEs BY R.I. OKUONGHAE Abstract.
More informationA new meshless method for steady-state heat conduction problems in anisotropic and inhomogeneous media
Archive of Applied Mechanics 74 25 563--579 Springer-Verlag 25 DOI 1.17/s419-5-375-8 A new meshless method for steady-state heat conduction problems in anisotropic and inhomogeneous media H. Wang, Q.-H.
More informationD. Shepard, Shepard functions, late 1960s (application, surface modelling)
Chapter 1 Introduction 1.1 History and Outline Originally, the motivation for the basic meshfree approximation methods (radial basis functions and moving least squares methods) came from applications in
More informationAnalysis of a Bivariate Risk Model
Jingyan Chen 1 Jiandong Ren 2 July 23, 2012 1 MSc candidate, Department of Statistical and Actuarial Sciences, The University of Western Ontario, London, Ontario, Canada. 2 Associate Professor, Department
More informationMFS with RBF for Thin Plate Bending Problems on Elastic Foundation
MFS with RBF for Thin Plate Bending Problems on Elastic Foundation Qing-Hua Qin, Hui Wang and V. Kompis Abstract In this chapter a meshless method, based on the method of fundamental solutions (MFS and
More informationA NOTE ON RATIONAL OPERATOR MONOTONE FUNCTIONS. Masaru Nagisa. Received May 19, 2014 ; revised April 10, (Ax, x) 0 for all x C n.
Scientiae Mathematicae Japonicae Online, e-014, 145 15 145 A NOTE ON RATIONAL OPERATOR MONOTONE FUNCTIONS Masaru Nagisa Received May 19, 014 ; revised April 10, 014 Abstract. Let f be oeprator monotone
More informationPrecise Large Deviations for Sums of Negatively Dependent Random Variables with Common Long-Tailed Distributions
This article was downloaded by: [University of Aegean] On: 19 May 2013, At: 11:54 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
More informationSpectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions
Applied Mathematical Sciences, Vol. 1, 2007, no. 5, 211-218 Spectral Collocation Method for Parabolic Partial Differential Equations with Neumann Boundary Conditions M. Javidi a and A. Golbabai b a Department
More informationOptimal stopping of a risk process when claims are covered immediately
Optimal stopping of a risk process when claims are covered immediately Bogdan Muciek Krzysztof Szajowski Abstract The optimal stopping problem for the risk process with interests rates and when claims
More information1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,
More informationMSc Dissertation topics:
.... MSc Dissertation topics: Omar Lakkis Mathematics University of Sussex Brighton, England November 6, 2013 Office Location: Pevensey 3 5C2 Office hours: Autumn: Tue & Fri 11:30 12:30; Spring: TBA. O
More informationUpper and lower bounds for ruin probability
Upper and lower bounds for ruin probability E. Pancheva,Z.Volkovich and L.Morozensky 3 Institute of Mathematics and Informatics, the Bulgarian Academy of Sciences, 3 Sofia, Bulgaria pancheva@math.bas.bg
More informationDirect method for variational problems by using hybrid of block-pulse and Bernoulli polynomials
Direct method for variational problems by using hybrid of block-pulse and Bernoulli polynomials M Razzaghi, Y Ordokhani, N Haddadi Abstract In this paper, a numerical method for solving variational problems
More informationT u
WANG LI-LIAN Assistant Professor Division of Mathematical Sciences, SPMS Nanyang Technological University Singapore, 637616 T 65-6513-7669 u 65-6316-6984 k lilian@ntu.edu.sg http://www.ntu.edu.sg/home/lilian?
More informationarxiv: v1 [math.ap] 11 Jan 2014
THE UNIFIED TRANSFORM FOR THE MODIFIED HELMHOLTZ EQUATION IN THE EXTERIOR OF A SQUARE A. S. FOKAS AND J. LENELLS arxiv:4.252v [math.ap] Jan 24 Abstract. The Unified Transform provides a novel method for
More informationOn Finite-Time Ruin Probabilities in a Risk Model Under Quota Share Reinsurance
Applied Mathematical Sciences, Vol. 11, 217, no. 53, 269-2629 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/ams.217.7824 On Finite-Time Ruin Probabilities in a Risk Model Under Quota Share Reinsurance
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 37: RBF Hermite Interpolation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationRuin Theory. A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science at George Mason University
Ruin Theory A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science at George Mason University by Ashley Fehr Bachelor of Science West Virginia University, Spring
More informationA Dynamic Contagion Process with Applications to Finance & Insurance
A Dynamic Contagion Process with Applications to Finance & Insurance Angelos Dassios Department of Statistics London School of Economics Angelos Dassios, Hongbiao Zhao (LSE) A Dynamic Contagion Process
More information