# ABSTRACT. Keywords: Dispersion; hybrid method; oscillatory problems; oscillatory solution; trigonometrically fitted ABSTRAK

Save this PDF as:

Size: px
Start display at page:

Download "ABSTRACT. Keywords: Dispersion; hybrid method; oscillatory problems; oscillatory solution; trigonometrically fitted ABSTRAK"

## Transcription

1 Sains Malaysiana 44(3)(2015): Zero-Dissipative Trigonometrically Fitted Hybrid Method for Numerical Solution of Oscillatory Problems (Kaedah Hibrid Penyuaian Trigonometri Lesapan-Sifar untuk Penyelesaian Berangka Masalah Berayun) YUSUF DAUDA JIKANTORO, FUDZIAH ISMAIL* & NORAZAK SENU ABSTRACT In this paper, an improved trigonometrically fitted zero-dissipative explicit two-step hybrid method with fifth algebraic order is derived. The method is applied to several problems where by the solutions are oscillatory in nature. Numerical results obtained are compared with existing methods in the scientific literature. The comparison shows that the new method is more effective and efficient than the existing methods of the same order. Keywords: Dispersion; hybrid method; oscillatory problems; oscillatory solution; trigonometrically fitted ABSTRAK Dalam kertas ini, suatu penyuaian trigonometri lesapan sifar kaedah hibrid dua langkah penambahbaikan peringkat kelima diterbitkan. Kaedah ini digunakan untuk beberapa masalah yang penyelesaiannya berayun. Keputusan berangka yang diperoleh dibandingkan dengan kaedah sedia ada dalam maklumat saintifik. Perbandingan tersebut menunjukkan kaedah yang yang baharu ini adalah lebih efektif dan cekap berbanding kaedah sedia ada dengan peringkat yang sama. Kata kunci: Kaedah hibrid; masalah berayun; penyelesaian berayun; penyuaian trigonmetri; serakan INTRODUCTION In the last few decades, there has been a growing interest in the research of new numerical techniques for approximating the solution of second order initial value problem of the form y" f (x, y), y(x 0 ) y 0, y'(x 0 ) y' 0, (1) which is independent on y' explicitly. This type of problem arises in different fields of science and engineering, which includes quantum mechanics, celestial mechanics, molecular dynamics, quantum chemistry, astrophysics, electronics and semi-discretizations of wave equation. Numerous numerical methods have been derived for approximating the solutions of (1), some of which are Runge-Kutta methods, Runge-Kutta Nyström methods and linear multistep methods. For the Runge-Kutta methods and other related methods specifically derived for approximating the solutions of first order initial value problems (IVPs), the second order IVPs need to be transformed to a system of first order IVPs so that the methods can be applied. In the quest for methods that best approximate the solution of (1), many authors considered different modifications on Runge-Kutta methods, multistep methods and Runge-Kutta Nyström methods, such work can be seen in Al-Khasawneh et al. (2007); Butcher (2008) Dormand et al. (1987), Franco (2006), Ming et al. (2012), Senu et al. (2010) and Simos (2002). Hybrid type methods related to multistep methods have been proposed by many authors for approximating the solutions of (1) (Simos 2012, 1999; Tsitouras 2002). Franco (2006) proposed explicit two-step hybrid methods up to algebraic order six with less computational cost by using the algebraic order conditions of twostep hybrid methods developed by Coleman (2003). In furtherance to this, Ahmad et al. (2013) proposed semiimplicit hybrid methods up to algebraic order five. In the numerical integration of oscillatory or periodic problems, consideration of dispersion and dissipation errors is very important. The higher is the order of the errors for a method, the more efficient it is for solving oscillatory problems. Highest order of dispersion achieved for constant coefficients two-step hybrid method in the literature is eight (Franco 2006). Dispersive of order infinity two-step hybrid methods which coefficients are variables depending on the frequency of the problem to be solved were derived by Ahmad et al. (2013), Dizicheh et al. (2012) and Fang and Wu (2007). In this paper we present a trigonometrically fitted two-step hybrid method based on zero-dissipative five stage fifth order hybrid method derived by Yusuf Dauda (2014). Again, the coefficient of the method depends on the frequency of the problem to be solved; hence the frequency of the problems must be known in advance. The issue of how to choose the frequencies in trigonometrically fitted techniques if the frequency is not known in advance is very difficult. This problem has been considered by Vanden Beghe et al. (2001) whereby the

2 474 frequency is tuned at every step by analyzing the local error. In their work, Ramos and Vigo-Aguiar (2010) conjectured that the frequency depends on the numerical method, the problem to be solved together with the initial value as well as the stepsize chosen. Hybrid method can be written as follows: Y i (1 + c i ) y n c i y n 1 +h 2 a ij f (t n + c j h, Y j ), i 1,, s, (2) Y n +1 2y n y n 1 +h 2 b i f (t n + c i h, Y j ), (3) where y n + 1 and are approximations to y(x n+1 ) and y(x n ) respectively. The paper is organized as follows: The dispersion and dissipation analysis of two-step hybrid methods are given in the next section. In section that follows, we look at the stability aspect of the method followed by the derivation of the method next. Numerical results and discussion is given in the next section and finally the conclusion is given in the last section. DISPERSION AND DISSIPATION ANALYSIS In this section we discuss the dispersion and dissipation of explicit hybrid methods. Consider the homogeneous test problem, y"(x) λ 2 y(x) for v > 0. By replacing f (x, y) λ 2 y into (2) and (3) we have: Y i (1 + c i )y n c i y n 1 h 2 a ij λ 2 y j, i 1,, s, (4) y n+1 2y n y n 1 h 2 b i λ 2 y i. (5) Alternatively (4) and (5) can be written in vector form as: Y (e + c)y n cy n 1 z 2 AY, y n+1 2y n y n 1 z 2 b T Y, where z hv Y (Y 1,, Y s ) T, (c 1,, c s ) T, and b e (1,, 1) T, A, z λh. From (4) we obtain: Y (I + z 2 A) -1 (e + c)y n (I + z 2 A) 1 cy n 1, (6) We substitute (6) into (5) and obtained: y n+1 (2 z 2 b T (I + z 2 A) 1 (e + c) y n (1 z 2 b T (I + z 2 A) 1 c)y n 1. (7) Rewrite (7) and the following recursion is obtained, where y n+1 T(z 2 )y n + D(z 2 )y n 1 0, (8) T(z 2 ) 2 z 2 b T (I + z 2 A) 1 (e + c) D(z 2 ) 1 z 2 b T (I + z 2 A) 1 c. The solution of the difference equation (8), (Van der Houwen & Sommeijer 1987) is given by, y n 2 c ρ m cos(ω + nϕ), (9) and the solution to the test problem is, y(x n ) 2 δ cos(ψ + nz), (10) where ρ is the amplification factor, ϕ is the phase and ω, δ, ψ are real constants. The definition formulated by Van de houwen and Sommeijer (1987) follows immediately. Definition 1: From (9) and (10), the quantity R(z) z ϕ is called dispersion error of the method. The method is said to have dispersion error of order q if R(z) O (z q+1 ). Furthermore, the quantity S(z) 1 ρ is called dissipation error of the method. And the method is said to be dissipative of order r if S(z) O (z r+1 ). STABILITY ANALYSIS In the analysis of stability of explicit hybrid methods, the test equation y"(x) λ 2 y(x) need to be solved using the method. Rewriting equation (8) gives the stability polynomial of the method as: ξ 2 T(z 2 ) ξ + D(z 2 ) 0. The magnitude of the dispersion and dissipation errors is an important feature for solving second order IVPs (1) with periodic or oscillating solutions. In such problems, it is also important that the numerical solution obtained from the difference in (8) should be periodic as it is the exact solution of the test equation y"(x) λ 2 y(x). This is true only if its coefficients satisfy the following conditions (Franco 2006),

3 475 D(z 2 ) 1, T(z 2 ) < 2, z 2 (0, ). where (0, ) is known as the interval of periodicity of the method. With this, the method is said to be zero-dissipative since S(z) 0, which implies that the order of dissipation of the method is infinity. On the other hand, if the order of dissipation is finite that is S(z) 0, then absolute stability of the method is guaranteed on the following condition; D(z 2 < 1, T(z 2 < 1 + D(z 2 ), z 2 (0, ), and the interval (0, ) is known as the interval of absolute stability of the method. The method derived in this paper is periodic because of the zero-dissipative property of the method. METHODS Trigonometrically fitted methods are generally derived to exactly approximate the solution of the IVPs whose solutions are linear combination of the functions: {x j e αx, x j e αx }, where α can be complex or a real number. Suppose G(x) e iαx, were i is an imaginary unit, is the solution of (1). Then applying (2) to G(x) generates the following recursive relations (Fang & Wu 2007), All subscripts run to m or less. To derive the fifth order five stage method, we substitute m 5 in the recursive relations to get, cos(c 3 z) 1 c 3 + c 3 cos(z) + z 2 (a 31 cos(z) + a 32 ) 0 (11) sin(c 3 z) c 3 sin(z) z 2 (a 31 sin(z)) 0 (12) cos(c 4 z) 1 c 4 +c 4 cos(z)+z 2 (a 41 cos(z)+a 42 cos(c 3 z)) 0 sin(c 4 z) c 4 sin(z) + z 2 ( a 41 sin(z) + a 43 sin(c 3 z)) 0 cos(c 5 z) 1 c 5 + c 5 cos(z) + z 2 (a 51 cos(z) + (13) (14) a 52 + a 53 cos(c 3 z) + a 54 cos(c 5 z) 0 (15) sin(c 5 z) c 5 sin(z) + z 2 ( a 51 sin(z) + a 53 sin(c 3 z) + a 54 sin(c 5 z)) 0 (16) 2cos(z) 2 + z 2 (b 1 cos(z) + b 2 + b 3 cos(c 3 z) + b 4 cos(c 4 z) + b 5 cos(c 5 z)) 0 (17) b 1 sin(z) + b 3 sin(c 3 z) + b 4 sin(c 4 z) + b 5 sin(c 5 z) 0 (18) cos(c i z) (1 + c i ) c i cos(z) z 2 sin(c i z) c i sin(z) z 2 sin(c j z), cos(c i z) 1 cos(c i z), sin(c i z) 0. cos(c j z), Equations (11)-(18) are the equations of order conditions for five stage fifth order trigonometrically fitted hybrid method that replaced equations of algebraic order conditions of the original method. To obtain the coefficients of the method we solve the system of eight equations, (11)-(18), together with some additional equations from algebraic order conditions of order five, namely; b 1 + b 2 + b 3 + b 4 + b 5 1, (19) The relations replace the equations of algebraic order conditions of two-step hybrid method, which can be solved to give the coefficients of a particular method based on existing coefficients for solving problem of the form (1). The algebraic order conditions of two-step hybrid methods derived by (Coleman 2003) is given: Order 1: Σ b i 1 Order 2: Σ b i c i 0 Order 3: Σ b i Order 4: Σ b i Σ b i a ij c j 0, b 1 + b 3 (20) b 1 + b 3 (21) These sums up to eleven equations with seventeen unknowns parameters. The equations are solved in terms of six free parameters (c 3 1, c 4 c 5 a 43 0, a 53 a 54 0), whose values are obtained from the coefficients of ), presented in Table 1 as obtained from (Dauda Yusuf 2014). Equations (19)-(21) are chosen to augment the updated (17)-(18) so that bi are not taken as free parameters. b 1 b 3 Order5: Σ b i Σ b i c i a ij c j b 5

4 476 TABLE 1. Coefficient of ): zero dissipative fifth order five stage sixth order dispersive hybrid method derived in Yusouf Dauda (2014) M cos(z)sin sin(z) sin ( z) sin z 2 sin(z) sin(z)cos(z) cos(z) sin ( z) z 2 cos(z) sin ( z) z 2 cos(z) sin z 2 cos ( z) sin(z) z 2 cos ( z) sin z 2 sin ( z) cos z 2 sin(z) cos z 2 sin ( z) z 2 sin cos ( z) sin sin ( z) sin(z) z 2 cos sin( z) sin(z) cos sin M cos(z) sin sin(z) sin sin(z)cos(z) cos(z) sin ( z) z 2 cos(z)sin(z) z 2 cos(z)sin ( z) z 2 cos(z) sin M 2 z 2 ( cos(z) sin ( z) sin(z) cos(z) cos(z) sin z 2 cos( z) sin z 2 cos sin ( z) cos ( z) sin(z) z 2 cos sin(z)

5 477 M 4 z 2 ( cos(z) sin ( z) sin(z) cos(z) cos(z) sin cos ( z) sin(z) cos ( z) sin sin ( z) sin (z) cos sin ( z) M 6 z 2 ( cos(z) sin ( z) sin(z)cos(z) cos(z) sin cos( z) sin(z) cos ( z) sin sin ( z) sin(z) cos sin ( z) cos sin(z) cos sin(z) sin sin M z 2 cos(z) sin ( z) z 2 cos(z) sin z 2 cos( z) sin(z) z 2 cos sin ( z) z 2 sin ( z) z 2 sin(z) z 2 cos sin ( z) z 2 cos sin(z) M cos(z) sin sin(z) sin z 2 sin(z) sin(z) cos(z) z 2 cos(z) sin(z) z 2 cos(z) sin z 2 sin, z 2 sin sin ( z) cos(z)sin ( z) sin(z) sin(z)cos(z) sin cos(z)sin. M 8 z 2 ( cos(z)sin( z) sin(z)cos(z) cos(z)sin ( z) sin(z) cos( z) sin sin( z) sin(z)

6 cos sin ( z) a cos sin(z) 2sin + sin(z) 3 sin(z) sin a 42 M 9 ( 8sin ( z)+sin(z))(14z z 2 cos(z)+30cos9z)). a 51 M 10 z 2 ( cos(z) sin( z) sin(z)cos(z) cos(z) sin cos ( z) sin(z) cos ( z) sin sin( z) sin(z) cos sin( z) a 52 It is important to note that the original method, that is ) needs to be recovered as z approaches zero. As, such, Taylor series expansion of the coefficients above is useful. That is b 1 b 2 b cos sin sin(z) b 4 b 5, M sin sin(z) z 2 sin(z), a 42 a 43 0, a 31 0, a 32 M cos sin(z) cos(z) sin sin(z) a 41 a z 2 cos(z)sin(z), a 52 a 31 0, a 32 It is clear that as z approaches zero the original method HM5(5,6) is recovered.

7 479 Next is to verify the algebraic order of the method. To check if the method is order five as claimed, we substitute the coefficients of the method into the algebraic order conditions up to order five and take the Taylor series expansion of each. Order 1: Order 2: Order3: Exact solution: y(x) λ π. Source: Dizicheh et al Problem 2 (Homogeneous problem) y"(x) 64y(x), y(0) 1, y'(0) 2, Exact solution is y(x) sin(8x) + cos(8x), λ 8 Source: Senu et al Problem 3 (Homogeneous problem) y"(x) 100y(x), y(0) 1, y'(0) 2 and the fitted frequency, λ 10 Exact solution is y(x) Source: Senu et al sin(10x) + cos(10x) Problem 4 (Two body problem) Order 4: Order5: Exact solution: y 1 (x) cos(x), y 2 (x) sin(x) λ 1 Source: Senu et al Problem 5 (Inhomogeneous problem) y"(x) y(x) + x, y(0) 1, y'(0) 2, Exact solution: y(x) sin(x) + cos(x) + x λ 1 Source: Al-Khasawneh et al Problem 6 (Duffing problem) It can be seen from above that as z tends to zero the algebraic order five conditions are recovered, which implies that the coefficients of trigonometrically fitted fifth order method satisfies algebraic order five conditions. RESULTS Presented in this section are numerical results of the trigonometrically fitted zero-dissipative fifth order hybrid method derived in the paper denoted by TFZDHM5. The method is applied to Problems 1-6 presented as follows: Problem 1 (Non-linear problem) y"(x) π sin(πx), y(0) 0, y'(0) 1. μ 5 is , μ 7 is , μ 9 <10 12 ; and λ 1. Exact solution: y(x) cos[(2i + 1)ux], where μ 1 is ; μ 3 is , μ 5 is Source: Van de Vyver 2005 The results are tabulated together with those of existing codes in the literature for the purpose of comparison. The interval of integration are taken to be 100, 1000 and 10000, respectively. Numerical results of trigonometrically fitted zerodissipative fifth order hybrid method applied to Problems 1-6 with different step sizes and integration intervals are

8 480 tabulated in this section. Both small and large integration intervals are considered to measure the stability of the methods when solving highly oscillatory problems. It can be seen from Tables 2-7 that the accuracy of TFZDHM5 diminishes as h grows smaller, especially for problems with small frequencies. This is because the TABLE 2. Accuracy Comparison of TFZDHM5 for problem 1 1/4 1/ ): zero dissipative fifth order five stage sixth order dispersive hybrid method derived in Yusouf Dauda (2014) TABLE 3. Accuracy Comparison of TFZDHM5 for problem ): zero dissipative fifth order five stage sixth order dispersive hybrid method derived in Yusouf Dauda (2014) TABLE 4. Accuracy Comparison of TFZDHM5 for problem ): zero dissipative fifth order five stage sixth order dispersive hybrid method derived in Yusouf Dauda (2014)

9 481 TABLE 5. Accuracy Comparison of TFZDHM5 for problem ): zero dissipative fifth order five stage sixth order dispersive hybrid method derived in Yusouf Dauda (2014) TABLE 6. Accuracy Comparison of TFZDHM5 for problem ): zero dissipative fifth order five stage sixth order dispersive hybrid method derived in Yusouf Dauda (2014) TABLE 7. Accuracy Comparison of TFZDHM5 for problem ): zero dissipative fifth order five stage sixth order dispersive hybrid method derived in Yusouf Dauda (2014)

10 482 method approaches the original method as h 0. For problems with large frequency, TFZDHM5 performs far better than the original method as well as other methods considered. CONCLUSION In this paper trigonometrically fitted hybrid method based on the existing zero-dissipative hybrid method was derived. The method was applied to highly oscillatory problems. The results obtained were compared with those of existing methods in the literature. From the numerical results obtained, we can conclude that the new trigonometrically fitted hybrid methods approximate the solution of highly oscillatory problems better than trigonometrically fitted hybrid methods that are based on higher dispersive and dissipative methods. Digits50 was used in MAPLE 16 environment for implementing all the methods on the problems. REFERENCES Ahmad, S.Z., Ismail, F., Senu, N. & Suleiman, M Zerodissipative phase-fitted hybrid method for solving oscillatory second order ordinary differential equations. Applied Mathematics and Computation 219(19): Al-Khasawneh, R.A., Ismail, F. & Suleiman, M Embedded diagonally implicit Runge Kutta Nyström 4(3) pair for solving special second-order IVPs. Applied Mathematics and Computation 190(2): Butcher, J.C Numerical Methods for Ordinary Differential Equations. New York: John Wiley & Sons. Coleman, J.P Order conditions for a class of two-step Methods for y f (x, y). IMA Journal of Numerical Analysis 23(2): Dizicheh, A.K., Ismail, F., Md. Arifin, N. & Abu Hassan, M A class of two-step hybrid trigonometrically fitted explicit Numerov-type method for second-order IVPs. Acta Comptare 1: Dormand, J.R., El-Mikkawy, M.E.A. & Prince, P.J Highorder embedded Runge-Kutta- Nystrom formulae. IMA Journal of Numerical Analysis 7(4): Fang, Y. & Wu, X A trigonometrically fitted explicit hybrid method for the numerical integration of orbital problems. Applied Mathematics and Computation 189(1): Franco, J.M A class of explicit two-step hybrid methods for second-order IVPs. Journal of Computational and Applied Mathematics 187(1): Ramos, H. & Vigo-Aguiar, J On the frequency choice in trigonometrically fitted methods. Applied Mathematics Letters 23: Ming, Q., Yang, Y. & Fang, Y An Optimized Runge-Kutta Method for the Numerical Solution of the Radial Schrödinger Equation. Mathematical Problems in Engineering 2012: Article ID Mohamad, M., Senu, N., Suleiman, M. & Ismail, F Fifth order explicit Runge-Kutta-Nyström methods for periodic initial value problems. World of Applied Sciences Journal 17: Senu, N., Suleiman, M., Ismail, F. & Othman, M Kaedah pasangan 4(3) Runge-Kutta-Nyström untuk masalah nilai awal berkala. Sains Malaysiana 39(4): Senu, N., Suleiman, M. & Ismail, F An embedded explicit Runge Kutta Nyström method for solving oscillatory problems. Physica Scripta 80(1): Simos, T.E Optimizing a hybrid two-step method for the numerical solution of the Schrödinger equation and related problems with respect to phase-lag. Journal of Applied Mathematics 2012: Article ID Simos, T.E Exponentially-fitted Runge-Kutta-Nyström method for the numerical solution of initial value problems with oscillating solutions. Applied Mathematics Letters 15(2): Simos, T.E Explicit eighth order methods for the numerical integration of initial-value problems with periodic or oscillating solutions. Computer Physics Communications 119(1): Tsitouras, Ch Explicit two-step methods for second-order linear IVPs. Computers & Mathematics with Applications 43(8): Van der Houwen, P.J. & Sommeijer, B.P Explicit Runge- Kutta (-Nyström) methods with reduced phase errors for computing oscillating solutions. SIAM Journal on Numerical Analysis 24(3): Van de Vyver, H A Runge-Kutta Nyström pair for the numerical integration of perturbed oscillators. Computer Physics Communications 167(2): Vanden Berghe, G., Ixaru L.G. & Van Daele, M Optimal implicit exponentially fitted Runge-Kutta methods. Computational Physics Comm. 140: Yusuf Dauda Jikantoro Numerical solution of special second initial value problems by hybrid type methods. MSc Thesis, Universiti Putra Malaysia (Unpublished). Department of Mathematics and Institute for Mathematical Research Universiti Putra Malaysia Serdang 43400, Selangor Darul Ehsan Malaysia *Corresponding author; Received: 7 January 2014 Accepted: 3 October 2014

### ABSTRACT. Keywords: Continuity; interpolation; monotonicity; rational Bernstein-Bézier ABSTRAK

Sains Malaysiana 40(10)(2011): 1173 1178 Improved Sufficient Conditions for Monotonic Piecewise Rational Quartic Interpolation (Syarat Cukup yang Lebih Baik untuk Interpolasi Kuartik Nisbah Cebis Demi

### ABSTRACT. Keywords: Continuity; interpolation; rational trigonometric spline; shape preserving ABSTRAK

Sains Malaysiana 45(0)(206): 557 563 A Univariate Rational Quadratic Trigonometric Interpolating Spline to Visualize Shaped Data (Sebuah Nisbah Kuadratik Trigonometri Univariat Menginterpolasikan Splin

### SOLUTION AND INTERPOLATION OF ONE-DIMENSIONAL HEAT EQUATION BY USING CRANK-NICOLSON, CUBIC SPLINE AND CUBIC B-SPLINE WAN KHADIJAH BINTI WAN SULAIMAN

SOLUTION AND INTERPOLATION OF ONE-DIMENSIONAL HEAT EQUATION BY USING CRANK-NICOLSON, CUBIC SPLINE AND CUBIC B-SPLINE WAN KHADIJAH BINTI WAN SULAIMAN A thesis submitted in fulfillment of the requirements

### Spline Solution Of Some Linear Boundary Value Problems

Applied Mathematics E-Notes, 8(2008, 171-178 c ISSN 1607-2510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Spline Solution Of Some Linear Boundary Value Problems Abdellah Lamnii,

### Review Initial Value Problems Euler s Method Summary

THE EULER METHOD P.V. Johnson School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 INITIAL VALUE PROBLEMS The Problem Posing a Problem 3 EULER S METHOD Method Errors 4 SUMMARY OUTLINE 1 REVIEW 2 INITIAL

### 4 Visualization and. Approximation

4 Visualization and Approximation b A slope field for the differential equation y tan(x + y) tan(x) tan(y). It is not always possible to write down an explicit formula for the solution to a differential

### 2 Computation with Floating-Point Numbers

2 Computation with Floating-Point Numbers 2.1 Floating-Point Representation The notion of real numbers in mathematics is convenient for hand computations and formula manipulations. However, real numbers

### Optimized Low Dispersion and Low Dissipation Runge- Kutta Algorithms in Computational Aeroacoustics

Appl. Math. Inf. Sci. 8, No. 1, 57-68 214 57 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/1.12785/amis/816 Optimized Low Dispersion Low Dissipation Runge- Kutta

### Introduction to Complex Analysis

Introduction to Complex Analysis George Voutsadakis 1 1 Mathematics and Computer Science Lake Superior State University LSSU Math 413 George Voutsadakis (LSSU) Complex Analysis October 2014 1 / 50 Outline

### Rounding Theorem the Possibility of Applying the Cryptosystems on the Decimal Numbers

Journal of Mathematics and Statistics 4 (): 5-20, 2008 ISSN 549-3644 2008 Science Publications Rounding Theorem the Possibility of Applying the Cryptosystems on the Decimal Numbers Rand Alfaris, Muhamad

### An Adaptive Stencil Linear Deviation Method for Wave Equations

211 An Adaptive Stencil Linear Deviation Method for Wave Equations Kelly Hasler Faculty Sponsor: Robert H. Hoar, Department of Mathematics ABSTRACT Wave Equations are partial differential equations (PDEs)

Quick Start with CASSY Lab Bi-05-05 About this manual This manual helps you getting started with the CASSY system. The manual does provide you the information you need to start quickly a simple CASSY experiment

### Computational Methods. Sources of Errors

Computational Methods Sources of Errors Manfred Huber 2011 1 Numerical Analysis / Scientific Computing Many problems in Science and Engineering can not be solved analytically on a computer Numeric solutions

### International Journal of Emerging Technologies in Computational and Applied Sciences (IJETCAS)

International Association of Scientific Innovation and Research (IASIR) (An Association Unifying the Sciences, Engineering, and Applied Research) International Journal of Emerging Technologies in Computational

### = w. w u. u ; u + w. x x. z z. y y. v + w. . Remark. The formula stated above is very important in the theory of. surface integral.

1 Chain rules 2 Directional derivative 3 Gradient Vector Field 4 Most Rapid Increase 5 Implicit Function Theorem, Implicit Differentiation 6 Lagrange Multiplier 7 Second Derivative Test Theorem Suppose

### Floating-point representation

Lecture 3-4: Floating-point representation and arithmetic Floating-point representation The notion of real numbers in mathematics is convenient for hand computations and formula manipulations. However,

### Edge Detections Using Box Spline Tight Frames

Edge Detections Using Box Spline Tight Frames Ming-Jun Lai 1) Abstract. We present a piece of numerical evidence that box spline tight frames are very useful for edge detection of images. Comparsion with

### A Diagonal Split-cell Model for the High-order Symplectic FDTD Scheme

PIERS ONLINE, VOL. 2, NO. 6, 2006 715 A Diagonal Split-cell Model for the High-order Symplectic FDTD Scheme Wei Sha, Xianliang Wu, and Mingsheng Chen Key Laboratory of Intelligent Computing & Signal Processing

### Truncation Errors. Applied Numerical Methods with MATLAB for Engineers and Scientists, 2nd ed., Steven C. Chapra, McGraw Hill, 2008, Ch. 4.

Chapter 4: Roundoff and Truncation Errors Applied Numerical Methods with MATLAB for Engineers and Scientists, 2nd ed., Steven C. Chapra, McGraw Hill, 2008, Ch. 4. 1 Outline Errors Accuracy and Precision

### Scientific Computing: Interpolation

Scientific Computing: Interpolation Aleksandar Donev Courant Institute, NYU donev@courant.nyu.edu Course MATH-GA.243 or CSCI-GA.22, Fall 25 October 22nd, 25 A. Donev (Courant Institute) Lecture VIII /22/25

### ENHANCEMENT OF UML-BASED WEB ENGINEERING FOR METAMODELS: HOMEPAGE DEVELOPMENT CASESTUDY KARZAN WAKIL SAID

ENHANCEMENT OF UML-BASED WEB ENGINEERING FOR METAMODELS: HOMEPAGE DEVELOPMENT CASESTUDY KARZAN WAKIL SAID A dissertation submitted in partial fulfillment of the requirements for the award of the degree

### COURSE: NUMERICAL ANALYSIS. LESSON: Methods for Solving Non-Linear Equations

COURSE: NUMERICAL ANALYSIS LESSON: Methods for Solving Non-Linear Equations Lesson Developer: RAJNI ARORA COLLEGE/DEPARTMENT: Department of Mathematics, University of Delhi Page No. 1 Contents 1. LEARNING

### ON THE TOTAL VERTEX IRREGULARITY STRENGTHS OF QUADTREES AND BANANA TREES

J. Indones. Math. Soc. Vol. 18, No. 1 (2012), pp. 31 36. ON THE TOTAL VERTEX IRREGULARITY STRENGTHS OF QUADTREES AND BANANA TREES Nurdin Department of Mathematics, Faculty of Mathematics and Natural Sciences,

### Hyperbola for Curvilinear Interpolation

Applied Mathematical Sciences, Vol. 7, 2013, no. 30, 1477-1481 HIKARI Ltd, www.m-hikari.com Hyperbola for Curvilinear Interpolation G. L. Silver 868 Kristi Lane Los Alamos, NM 87544, USA gsilver@aol.com

### Euler s Methods (a family of Runge- Ku9a methods)

Euler s Methods (a family of Runge- Ku9a methods) ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an

### MATH3016: OPTIMIZATION

MATH3016: OPTIMIZATION Lecturer: Dr Huifu Xu School of Mathematics University of Southampton Highfield SO17 1BJ Southampton Email: h.xu@soton.ac.uk 1 Introduction What is optimization? Optimization is

### Higher. The Wave Equation. The Wave Equation 146

Higher Mathematics UNIT OUTCOME 4 The Wave Equation Contents The Wave Equation 146 1 Expressing pcosx + qsinx in the form kcos(x a 146 Expressing pcosx + qsinx in other forms 147 Multiple Angles 148 4

### PS Geometric Modeling Homework Assignment Sheet I (Due 20-Oct-2017)

Homework Assignment Sheet I (Due 20-Oct-2017) Assignment 1 Let n N and A be a finite set of cardinality n = A. By definition, a permutation of A is a bijective function from A to A. Prove that there exist

### Exponential Membership Functions in Fuzzy Goal Programming: A Computational Application to a Production Problem in the Textile Industry

American Journal of Computational and Applied Mathematics 2015, 5(1): 1-6 DOI: 10.5923/j.ajcam.20150501.01 Exponential Membership Functions in Fuzzy Goal Programming: A Computational Application to a Production

### Lagrange Multipliers and Problem Formulation

Lagrange Multipliers and Problem Formulation Steven J. Miller Department of Mathematics and Statistics Williams College Williamstown, MA 01267 Abstract The method of Lagrange Multipliers (and its generalizations)

### Floating-point numbers. Phys 420/580 Lecture 6

Floating-point numbers Phys 420/580 Lecture 6 Random walk CA Activate a single cell at site i = 0 For all subsequent times steps, let the active site wander to i := i ± 1 with equal probability Random

### A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions

A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions Nira Dyn Michael S. Floater Kai Hormann Abstract. We present a new four-point subdivision scheme that generates C 2 curves.

### 2 Second Derivatives. As we have seen, a function f (x, y) of two variables has four different partial derivatives: f xx. f yx. f x y.

2 Second Derivatives As we have seen, a function f (x, y) of two variables has four different partial derivatives: (x, y), (x, y), f yx (x, y), (x, y) It is convenient to gather all four of these into

### ODE IVP. An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable:

Euler s Methods ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an initial value/condition (i.e., value

### . Tutorial Class V 3-10/10/2012 First Order Partial Derivatives;...

Tutorial Class V 3-10/10/2012 1 First Order Partial Derivatives; Tutorial Class V 3-10/10/2012 1 First Order Partial Derivatives; 2 Application of Gradient; Tutorial Class V 3-10/10/2012 1 First Order

### Math 340 Fall 2014, Victor Matveev. Binary system, round-off errors, loss of significance, and double precision accuracy.

Math 340 Fall 2014, Victor Matveev Binary system, round-off errors, loss of significance, and double precision accuracy. 1. Bits and the binary number system A bit is one digit in a binary representation

### Simulation in Computer Graphics. Particles. Matthias Teschner. Computer Science Department University of Freiburg

Simulation in Computer Graphics Particles Matthias Teschner Computer Science Department University of Freiburg Outline introduction particle motion finite differences system of first order ODEs second

### Third Order WENO Scheme on Three Dimensional Tetrahedral Meshes

COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 5, No. 2-4, pp. 86-848 Commun. Comput. Phys. February 29 Third Order WENO Scheme on Three Dimensional Tetrahedral Meshes Yong-Tao Zhang 1, and Chi-Wang Shu

### ON SWELL COLORED COMPLETE GRAPHS

Acta Math. Univ. Comenianae Vol. LXIII, (1994), pp. 303 308 303 ON SWELL COLORED COMPLETE GRAPHS C. WARD and S. SZABÓ Abstract. An edge-colored graph is said to be swell-colored if each triangle contains

### = f (a, b) + (hf x + kf y ) (a,b) +

Chapter 14 Multiple Integrals 1 Double Integrals, Iterated Integrals, Cross-sections 2 Double Integrals over more general regions, Definition, Evaluation of Double Integrals, Properties of Double Integrals

### Dynamic Computation of Runge-Kutta s Fourth-Order Algorithm for First and Second Order Ordinary Differential Equation Using Java

211 Dynamic Computation of Runge-Kutta s Fourth-Order Algorithm for First and Second Order Ordinary Differential Equation Using Java Adesola O. Anidu Samson A. Arekete Ayomide O. Adedayo and Adekunle O.

### 1. Suppose that the equation F (x, y, z) = 0 implicitly defines each of the three variables x, y, and z as functions of the other two:

Final Solutions. Suppose that the equation F (x, y, z) implicitly defines each of the three variables x, y, and z as functions of the other two: z f(x, y), y g(x, z), x h(y, z). If F is differentiable

### Locally convex topological vector spaces

Chapter 4 Locally convex topological vector spaces 4.1 Definition by neighbourhoods Let us start this section by briefly recalling some basic properties of convex subsets of a vector space over K (where

### An ε-uniform Initial Value Technique For Convection-Diffusion Singularly Perturbed Problems

An -Uniform Initial Value Technique For Convection-Diffusion Singularly Perturbed Problems RKBawa and Vinod Kumar Abstract In this paper, we have proposed an -uniform initial value technique for singularly

### used to describe all aspects of imaging process: input scene f, imaging system O, and output image g.

SIMG-716 Linear Imaging Mathematics I 03 1 Representations of Functions 1.1 Functions Function: rule that assigns one numerical value (dependent variable) to each of a set of other numerical values (independent

### CLEP Pre-Calculus. Section 1: Time 30 Minutes 50 Questions. 1. According to the tables for f(x) and g(x) below, what is the value of [f + g]( 1)?

CLEP Pre-Calculus Section : Time 0 Minutes 50 Questions For each question below, choose the best answer from the choices given. An online graphing calculator (non-cas) is allowed to be used for this section..

### used to describe all aspects of imaging process: input scene f, imaging system O, and output image g.

SIMG-716 Linear Imaging Mathematics I 03 1 Representations of Functions 1.1 Functions Function: rule that assigns one numerical value (dependent variable) to each of a set of other numerical values (independent

### A Robust Optimum Response Surface Methodology based On MM-estimator

A Robust Optimum Response Surface Methodology based On MM-estimator, 2 HABSHAH MIDI,, 2 MOHD SHAFIE MUSTAFA,, 2 ANWAR FITRIANTO Department of Mathematics, Faculty Science, University Putra Malaysia, 434,

### Hybrid Dynamic Iterations for the Solution of Initial Value Problems

Hybrid Dynamic Iterations for the Solution of Initial Value Problems Yanan Yu and Ashok Srinivasan Department of Computer Science, Florida State University, Tallahassee, FL32306 Abstract Many scientific

### Power Series and Summation

Power Series and Summation Prof. Dr. Wolfram Koepf Department of Mathematics University of Kassel koepf@mathematik.uni-kassel.de http://www.mathematik.uni-kassel.de/~koepf Tutorial ISSAC 2004 Santander,

### Higher-Order Newton-Cotes Formulas

Journal of Mathematics and Statistics 6 (): 19-0, 010 ISSN 159-6 010 Science Publications Higher-Order Newton-Cotes Formulas Pedro Americo Almeida Magalhaes Junior and Cristina Almeida Magalhaes Pontificia

### Lacunary Interpolation Using Quartic B-Spline

General Letters in Mathematic, Vol. 2, No. 3, June 2017, pp. 129-137 e-issn 2519-9277, p-issn 2519-9269 Available online at http:\\ www.refaad.com Lacunary Interpolation Using Quartic B-Spline 1 Karwan

### Midterm Exam Solutions

Midterm Exam Solutions Computer Vision (J. Košecká) October 27, 2009 HONOR SYSTEM: This examination is strictly individual. You are not allowed to talk, discuss, exchange solutions, etc., with other fellow

### NEWTON BASED HOMOTOPY OPTIMIZATION METHOD FOR SOLVING GLOBAL OPTIMIZATION PROBLEM AHMAD ZHARIF SALAMI BIN MOHD TAHA

NEWTON BASED HOMOTOPY OPTIMIZATION METHOD FOR SOLVING GLOBAL OPTIMIZATION PROBLEM AHMAD ZHARIF SALAMI BIN MOHD TAHA A dissertation submitted in partial fulfillment of the requirements for the award of

### The Eccentric Digraph of a Lintang Graph

The Eccentric Digraph of a Lintang Graph Tri Atmojo Kusmayadi * dan Fauziyyah Fathmawatie Department of Mathematics Faculty of Mathematics and Natural Sciences Sebelas Maret University, Surakarta, Indonesia

### arxiv: v1 [math.co] 17 Jan 2014

Regular matchstick graphs Sascha Kurz Fakultät für Mathematik, Physik und Informatik, Universität Bayreuth, Germany Rom Pinchasi Mathematics Dept., Technion Israel Institute of Technology, Haifa 2000,

### Module 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 6:

file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture6/6_1.htm 1 of 1 6/20/2012 12:24 PM The Lecture deals with: ADI Method file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture6/6_2.htm 1 of 2 6/20/2012

### Section 7.6 Graphs of the Sine and Cosine Functions

Section 7.6 Graphs of the Sine and Cosine Functions We are going to learn how to graph the sine and cosine functions on the xy-plane. Just like with any other function, it is easy to do by plotting points.

### A Toolbox of Level Set Methods

A Toolbox of Level Set Methods Ian Mitchell Department of Computer Science University of British Columbia http://www.cs.ubc.ca/~mitchell mitchell@cs.ubc.ca research supported by the Natural Science and

### CS321. Introduction to Numerical Methods

CS31 Introduction to Numerical Methods Lecture 1 Number Representations and Errors Professor Jun Zhang Department of Computer Science University of Kentucky Lexington, KY 40506 0633 August 5, 017 Number

### GPU Implementation of Implicit Runge-Kutta Methods

GPU Implementation of Implicit Runge-Kutta Methods Navchetan Awasthi, Abhijith J Supercomputer Education and Research Centre Indian Institute of Science, Bangalore, India navchetanawasthi@gmail.com, abhijith31792@gmail.com

### FINITE ELEMENT INVESTIGATION ON THE STRENGTH OF SEMI-RIGID EXTENDED END PLATE STEEL CONNECTION USING LUSAS SOFTWARE MOHD MAIZIZ BIN FISHOL HAMDI

FINITE ELEMENT INVESTIGATION ON THE STRENGTH OF SEMI-RIGID EXTENDED END PLATE STEEL CONNECTION USING LUSAS SOFTWARE MOHD MAIZIZ BIN FISHOL HAMDI Faculty of Civil Engineering Universiti Teknologi Malaysia

### Numerical Analysis of Shock Tube Problem by using TVD and ACM Schemes

Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics

### f xx (x, y) = 6 + 6x f xy (x, y) = 0 f yy (x, y) = y In general, the quantity that we re interested in is

1. Let f(x, y) = 5 + 3x 2 + 3y 2 + 2y 3 + x 3. (a) Final all critical points of f. (b) Use the second derivatives test to classify the critical points you found in (a) as a local maximum, local minimum,

### correlated to the Michigan High School Mathematics Content Expectations

correlated to the Michigan High School Mathematics Content Expectations McDougal Littell Algebra 1 Geometry Algebra 2 2007 correlated to the STRAND 1: QUANTITATIVE LITERACY AND LOGIC (L) STANDARD L1: REASONING

### Vectors. Section 1: Lines and planes

Vectors Section 1: Lines and planes Notes and Examples These notes contain subsections on Reminder: notation and definitions Equation of a line The intersection of two lines Finding the equation of a plane

### What is the Optimal Bin Size of a Histogram: An Informal Description

International Mathematical Forum, Vol 12, 2017, no 15, 731-736 HIKARI Ltd, wwwm-hikaricom https://doiorg/1012988/imf20177757 What is the Optimal Bin Size of a Histogram: An Informal Description Afshin

### DYNAMIC TIMETABLE GENERATOR USING PARTICLE SWARM OPTIMIZATION (PSO) METHOD TEH YUNG CHUEN UNIVERSITY MALAYSIA PAHANG

DYNAMIC TIMETABLE GENERATOR USING PARTICLE SWARM OPTIMIZATION (PSO) METHOD TEH YUNG CHUEN UNIVERSITY MALAYSIA PAHANG DYNAMIC TIMETABLE GENERATOR USING PARTICLE SWARM OPTIMIZATION (PSO) METHOD TEH YUNG

### On the order of accuracy and numerical performance of two classes of finite volume WENO schemes

On the order of accuracy and numerical performance of two classes of finite volume WENO schemes Rui Zhang, Mengping Zhang and Chi-Wang Shu November 29, 29 Abstract In this paper we consider two commonly

### Numerical Methods for (Time-Dependent) HJ PDEs

Numerical Methods for (Time-Dependent) HJ PDEs Ian Mitchell Department of Computer Science The University of British Columbia research supported by National Science and Engineering Research Council of

### r v i e w o f s o m e r e c e n t d e v e l o p m

O A D O 4 7 8 O - O O A D OA 4 7 8 / D O O 3 A 4 7 8 / S P O 3 A A S P - * A S P - S - P - A S P - - - - L S UM 5 8 - - 4 3 8 -F 69 - V - F U 98F L 69V S U L S UM58 P L- SA L 43 ˆ UéL;S;UéL;SAL; - - -

### Recursive Estimation

Recursive Estimation Raffaello D Andrea Spring 28 Problem Set : Probability Review Last updated: March 6, 28 Notes: Notation: Unless otherwise noted, x, y, and z denote random variables, p x denotes the

### Data Analysis 3. Support Vector Machines. Jan Platoš October 30, 2017

Data Analysis 3 Support Vector Machines Jan Platoš October 30, 2017 Department of Computer Science Faculty of Electrical Engineering and Computer Science VŠB - Technical University of Ostrava Table of

### A Deterministic Dynamic Programming Approach for Optimization Problem with Quadratic Objective Function and Linear Constraints

A Deterministic Dynamic Programming Approach for Optimization Problem with Quadratic Objective Function and Linear Constraints S. Kavitha, Nirmala P. Ratchagar International Science Index, Mathematical

### UNIVERSITI PUTRA MALAYSIA RELIABILITY PERFORMANCE EVALUATION AND INTEGRATION OF ROUTING ALGORITHM IN SHUFFLE EXCHANGE WITH MINUS ONE STAGE

UNIVERSITI PUTRA MALAYSIA RELIABILITY PERFORMANCE EVALUATION AND INTEGRATION OF ROUTING ALGORITHM IN SHUFFLE EXCHANGE WITH MINUS ONE STAGE NUR ARZILAWATI BINTI MD YUNUS FSKTM 2012 5 RELIABILITY PERFORMANCE

### Structure Based Classification and Kinematic Analysis of Six-Joint Industrial Robotic Manipulators

Structure Based Classification and Kinematic Analysis of Six-Joint Industrial Robotic Manipulators 5 Tuna Balkan, M. Kemal Özgören and M. A. Sahir Arıkan Open Access Database www.i-techonline.com 1. Introduction

### A Subdivision-Regularization Framework for Preventing Over Fitting of Data by a Model

Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 178-190 Applications and Applied Mathematics: An International Journal (AAM) A Subdivision-Regularization

### MATH 2400, Analytic Geometry and Calculus 3

MATH 2400, Analytic Geometry and Calculus 3 List of important Definitions and Theorems 1 Foundations Definition 1. By a function f one understands a mathematical object consisting of (i) a set X, called

### Finite element solution of multi-scale transport problems using the least squares based bubble function enrichment

Finite element solution of multi-scale transport problems using the least squares based bubble function enrichment A. Yazdani a, V. Nassehi b1 a Cranfield University, School of Applied Sciences, Cranfield,

### IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL RESERVOIR SIMULATIONS USING THE CIP SCHEME WITH THIRD-ORDER ACCURACY

PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL

### Math 253, Section 102, Fall 2006 Practice Final Solutions

Math 253, Section 102, Fall 2006 Practice Final Solutions 1 2 1. Determine whether the two lines L 1 and L 2 described below intersect. If yes, find the point of intersection. If not, say whether they

### High Order Weighted Essentially Non-Oscillatory Schemes for Convection. Dominated Problems. Chi-Wang Shu 1

High Order Weighted Essentially Non-Oscillatory Schemes for Convection Dominated Problems Chi-Wang Shu Division of Applied Mathematics, Brown University, Providence, Rhode Island 09 ABSTRACT High order

### SYSTEMS OF NONLINEAR EQUATIONS

SYSTEMS OF NONLINEAR EQUATIONS Widely used in the mathematical modeling of real world phenomena. We introduce some numerical methods for their solution. For better intuition, we examine systems of two

### BACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION

BACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION TODD F. DUPONT AND YINGJIE LIU Abstract. We propose a method that significantly

Turbulence et Génération de Bruit Equipe de recherche du Centre Acoustique LMFA, UMR CNRS 5509, Ecole Centrale de Lyon Simulation Numérique en Aéroacoustique Institut Henri Poincaré - 16 novembre 2006

### SOLVING SYSTEMS OF LINEAR INTERVAL EQUATIONS USING THE INTERVAL EXTENDED ZERO METHOD AND MULTIMEDIA EXTENSIONS

Please cite this article as: Mariusz Pilarek, Solving systems of linear interval equations using the "interval extended zero" method and multimedia extensions, Scientific Research of the Institute of Mathematics

### Primes in Classes of the Iterated Totient Function

1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 11 (2008), Article 08.1.2 Primes in Classes of the Iterated Totient Function Tony D. Noe 14025 NW Harvest Lane Portland, OR 97229 USA noe@sspectra.com

### 4 Parametrization of closed curves and surfaces

4 Parametrization of closed curves and surfaces Parametrically deformable models give rise to the question of obtaining parametrical descriptions of given pixel or voxel based object contours or surfaces,

### Discrete Optimization. Lecture Notes 2

Discrete Optimization. Lecture Notes 2 Disjunctive Constraints Defining variables and formulating linear constraints can be straightforward or more sophisticated, depending on the problem structure. The

### Development of Time-Dependent Queuing Models. in Non-Stationary Condition

The 11th Asia Pacific Industrial Engineering and Management Systems Conference Development of Time-Dependent Queuing Models in Non-Stationary Condition Nur Aini Masruroh 1, Subagyo 2, Aulia Syarifah Hidayatullah

### Partial Differential Equations

Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described

### Lecture 5: The Halting Problem. Michael Beeson

Lecture 5: The Halting Problem Michael Beeson Historical situation in 1930 The diagonal method appears to offer a way to extend just about any definition of computable. It appeared in the 1920s that it

### Computational Physics PHYS 420

Computational Physics PHYS 420 Dr Richard H. Cyburt Assistant Professor of Physics My office: 402c in the Science Building My phone: (304) 384-6006 My email: rcyburt@concord.edu My webpage: www.concord.edu/rcyburt

### Generating Functions

6.04/8.06J Mathematics for Computer Science Srini Devadas and Eric Lehman April 7, 005 Lecture Notes Generating Functions Generating functions are one of the most surprising, useful, and clever inventions

### Recent Developments in Model-based Derivative-free Optimization

Recent Developments in Model-based Derivative-free Optimization Seppo Pulkkinen April 23, 2010 Introduction Problem definition The problem we are considering is a nonlinear optimization problem with constraints:

### SLANTING EDGE METHOD FOR MODULATION TRANSFER FUNCTION COMPUTATION OF X-RAY SYSTEM FARHANK SABER BRAIM UNIVERSITI TEKNOLOGI MALAYSIA

SLANTING EDGE METHOD FOR MODULATION TRANSFER FUNCTION COMPUTATION OF X-RAY SYSTEM FARHANK SABER BRAIM UNIVERSITI TEKNOLOGI MALAYSIA SLANTING EDGE METHOD FOR MODULATION TRANSFER FUNCTION COMPUTATION OF

### 10. Network dimensioning

Partly based on slide material by Samuli Aalto and Jorma Virtamo ELEC-C7210 Modeling and analysis of communication networks 1 Contents Introduction Parameters: topology, routing and traffic Dimensioning

### Computer Vision. Coordinates. Prof. Flávio Cardeal DECOM / CEFET- MG.

Computer Vision Coordinates Prof. Flávio Cardeal DECOM / CEFET- MG cardeal@decom.cefetmg.br Abstract This lecture discusses world coordinates and homogeneous coordinates, as well as provides an overview