A COMPLEX VARIABLE BOUNDARY ELEMENT METHOD FOR AN ELLIPTIC PARTIAL DIFFERENTIAL EQUATION WITH VARIABLE COEFFICIENTS
|
|
- Deborah Harper
- 5 years ago
- Views:
Transcription
1 A COMPLEX VARIABLE BOUNDARY ELEMENT METHOD FOR AN ELLIPTIC PARTIAL DIFFERENTIAL EQUATION WITH VARIABLE COEFFICIENTS Y. S. PARK a adw.t.ang b Computatioal ad Mathematical Scieces Group Faculty of Iformatio Techology Uiversiti Malaysia Sarawak Kota Samaraha, Malaysia Abstract A boudary elemet method based o the Cauchy itegral formulae is proposed for the umerical solutio of a boudary value problem govered by a secod order elliptic partial differetial equatio with variable coefficiets. The boudary value problem has applicatios i egieerig problems ivolvig ohomogeeous media. The method reduces the boudary value problem to the task of solvig a system of liear algebraic equatios. It ca be easily implemeted o the computer. Keywords: complex variable boudary elemet method, elliptic partial differetial equatio with variable coefficiets. This is a preprit of the article i Commuicatios i Numerical Methods i Egieerig 16 (2000) The URL of the publisher is: 1
2 1 INTRODUCTION The secod order elliptic partial differetial equatio (PDE) µ K(x) φ + µ K(x) φ =0, φ = φ(x, y), (1) x x y y where K(x) is a give fuctio of the Cartesia coordiate x such that K(x) > 0 i the domai of iterest, has extesive egieerig applicatios i structural ceramics, modelig of mechaics of soils ad other particulate materials (see e.g. Brea [1] ad Gibso [2]). A boudary value problem (BVP) of iterest is to solve (1) i a fiite regio R bouded by a simple closed curve C (o the 0xy plae) subject to the coditios φ = p(x, y) for (x, y) D, (2) φ = q(x, y) for (x, y) E, (3) where p ad q are suitably prescribed fuctios, D ad E are o-itersectig curves such that C = D E, φ/ = φ ad =[ 1, 2 ] is the uit ormal outward vector to E. Clemets [3] had derived a geeral solutio of (1) i the form ( ) 1 X φ(x, y) = p Re f (x)φ (z), (4) K(x) where i = 1adz = x + iy. The fuctios f are particular solutios of the recurrece equatios 2 df +1 dx with f 0 (x) =1ad =0 = d2 f dx 2 + Λ(x)f for =0, 1, 2, 3,, (5) Λ(x) = d dx " µdk 2 /(4K)# / dx µ dk. (6) dx 2
3 The complex fuctios Φ (z) are holomorphic i the domai of iterest ad satisfy dφ dz = Φ 1(z) for =1, 2, 3,. (7) I the preset article, the solutio (4) together with the Cauchy itegral fromulae is applied to derive a umerical method for solvig approximately the BVP defied by (1), (2) ad (3). The method requires oly the boudary C of the domai R to be discretized ad is kow as the complex variable boudary elemet method (CVBEM). It reduces the BVP to solvig a system of liear algebraic equatios. The coefficiets of the liear algebraic equatios are easy to compute. Hece the method ca be easily implemeted o the computer. The proposed method should offer itself as a iterestig alterative to earlier BEM approaches (see e.g. Clemets [3] ad Ragogi [4]) for solvig certai BVPs govered by (1). The CVBEM was origially itroduced by Hromadka ad Lai [5] for BVPs govered by the two-dimesioal Laplace equatio, i.e. for the case where K(x) is a costat fuctio i (1). Recetly, Ag ad Park [6] developed adifferet versio of the method for a more geeral system of secod order elliptic PDEs with costat coefficiets. The approach of [6] which differs from that of [5] i the approximatio of the relevat complex fuctios ad i the treatmet of the boudary coditios is adopted here to solve the BVP uder cosideratio. 2 CVBEM Firstly, the boudary C is discretized as follows. Place M well-spaced out poits (x (1),y (1) ), (x (2),y (2) ),, (x (M 1),y (M 1) )ad(x (M),y (M) )iaaticlockwise order o C. Let C (k) be the straight lie segmet from (x (k),y (k) ) to (x (k+1),y (k+1) )(k =1, 2,,M), where (x (M+1),y (M+1) )=(x (1),y (1) ). We the make the approximatio C C (1) C (2) C (M 1) C (M). (8) 3
4 Repeated applicatio of the recurrece relatio (7) gives Z z 1 Φ (z) = (z t) 1 Φ 0 (t)dt for =1, 2,, (9) ( 1) a where a is a suitably chose complex umber. Let z (k) = x (k) + iy (k) (k =1, 2,,M +1), bz (p) = ξ (p) + iη (p), ξ (p) = (x (p) + x (p+1) )/2 adη (p) =(y (p) + y (p+1) )/2 (p =1, 2,,M). I (9), if we take a = z (1),z= bz (p) ad do the itegratio from a to z alog the path C (1) C (2) C (p 1) C b(p) (where C b(p) is the straight lie segmet from z (p) to bz (p) ), the Φ (bz (p) )= ( 1 Xp 1 ( 1) Z z (k+1) z (k) Z bz (p) + z (p) (bz (p) t) 1 Φ 0 (t)dt (bz (p) t) 1 Φ 0 (t)dt for p =1, 2,,M ad =1, 2,, (10) For t C (k), if we expad Φ 0 (t) as a Taylor-Maclauri series about t = bz (k), i.e. if we write Φ 0 (t) =Φ 0 (bz (k) )+(t bz (k) )Φ 0 0 (bz(k) )+ 1 2 (t bz(k) ) 2 Φ 00 0 (bz(k) )+, (11) the, after igorig terms whose magitudes are O( z (k+1) z (k) 2 )orhigher order, we fid that (10) gives (approximately) ) Φ (bz (p) )= px Γ (pk) Φ 0 (bz (k) ) for p =1, 2,,M ad =0, 1, 2,, (12) where Γ (pk) 0 = δ pk ad Ã Γ (pk) 1 = (1 δ pk ) ( 1) Z z (k+1) z (k) (bz (p) t) 1 dt + δ pk Z bz (p) z (p) (bz (p) t) 1 dt for =1, 2, 3,. (13) 4
5 Notice that Γ (pk) ca be easily evaluated. If we write Φ 0 (bz (p) )=u (p) + iv (p) where u (p) ad v (p) are real costats the use of (12) together with (4) i the coditio (2) yields ( ) where α (pk) px = u (k) X =0 f (ξ (p) )α (pk) v (k) X =0 f (ξ (p) )β (pk) q K(ξ (p) )p(ξ (p), η (p) ) if φ is specified over C (p), (14) ad β (pk) Now, from (4) ad (7), we fid that à ( ) 1 X p Re f (x)φ (z) K(x) =Re =0 ( ( 1 + i 2 ) p K(x) Φ 0 0(z)+ are real parameters give by Γ (pk) =0 = α (pk) + iβ (pk). à #) X f +1 (x) Φ (z) " 1 p + F (x) + i 2f +1 (x) p, K(x) K(x) (15) ³ where F (x) =d f (x)/ p K(x) /dx for =0, 1, 2,. Thus, to deal with the coditio i (3), we are required to evaluate Φ 0 0(z) at z o the boudary. We shall apply the Cauchy itegral formula M 2πiΦ 0 0 (z) = X Z Φ 0 (t)dt for z R. (16) (t z) 2 C (k) If we substitute (11) ito (16) ad igore terms whose magitudes are O( z (k+1) z (k) 2 ) or higher order, we obtai the approximatio M 2πiΦ 0 0 (z) = X Φ 0 (bz (k) ) q(z (k),z (k+1),z)+ir(z (k),z (k+1),z) for a fixed z R, where q ad r are real parameters defied by (17) q(a, b, c)+ir(a, b, c) = 1 b c + 1 a c. (18) 5
6 However, if we repeat the same exercise but with z bz (p) (the midpoit of the p-th straight lie segmet), we fid that πiφ 0 0 (bz(p) )= Φ 0 (bz (k) ) q(z (k),z (k+1), bz (p) )+ir(z (k),z (k+1), bz (p) ). (19) Notice the differece i the factors that is multiplied to the derivative Φ 0 0 i (17) ad (19). With (12), (15) ad (19), the coditio i (4) gives à u (k) (p) 2 q(z (k),z (k+1), bz (p) )+ (p) 1 r(z (k),z (k+1), bz (p) ) π p K(ξ (p) ) + + px where [ (p) 1, (p) σ (p) u (k) X =0 =0 α (pk) µ (p) β(pk) σ (p) à v (k) (p) 1 q(z (k),z (k+1), bz (p) ) (p) 2 r(z (k),z (k+1), bz (p) ) π p K(ξ (p) ) px X v (k) α (pk) σ (p) + β(pk) µ (p) = q(ξ (p), η (p) ) if φ is specified over C(p), (20) 2 ] is the uit ormal (outward) vector to C (p) ad µ (p) ad are real parameters defied by " à # µ (p) + iσ(p) = (p) f +1 (ξ (p) ) 1 p K(ξ (p) ) + F (ξ (p) ) + i(p) 2 f +1 (ξ (p) ) p. (21) K(ξ (p) ) The equatios (14) ad (20) give M equatios i 2M ukows u (k) ad v (k) (k =1, 2,,M). To obtai aother M equatios, we apply the Cauchy itegral formula 2πiΦ 0 (z) = Z C (k) Φ 0 (t)dt (t z) for z R. (22) 6
7 Proceedig as before, i.e. substitutig (11) ito (22) ad igorig terms whose magitudes are O( z (k+1) z (k) 2 )orhigherorder,weobtai 2πiΦ 0 (z) = (u (k) + iv (k) ) γ(z (k),z (k+1),z)+iθ(z (k),z (k+1),z) for z R, (23) where γ ad θ are real parameters defied by γ(a, b, c) =l b c l a c Θ(a, b, c) if Θ(a, b, c) [ π, π] θ(a, b, c) = Θ(a, b, c)+2π if Θ(a, b, c) [ 2π, π) Θ(a, b, c) 2π if Θ(a, b, c) (π, 2π] Θ(a, b, c) =Arg(b c) Arg(a c), (24) where Arg(z) deotes the pricipal argumet of the complex umber z. If thesolutiodomaiiscovexishape, θ(a, b, c) ca be calculated directly from à θ(a, b, c) =cos 1 b c 2 + a c 2 b a 2. (25) 2 b c a c We may take the real part of equatio (23) ad let z bz (p) to obtai 2πv (p) = u (k) γ(z (k),z (k+1), bz (p) ) v (k) θ(z (k),z (k+1), bz (p) ) ª for p =1, 2,,M. (26) The equatios (14), (20) ad (26) costitute 2M equatios from which the ukows u (k) ad v (k) (k =1, 2,,M) ca be determied. Oce the ukows are determied, we ca calculate Φ (z) ( =1, 2, )atthepoit z = bz (k) (k =1, 2,,M) usig (12) ad the at the poit z = c i the iterior of R via 2πiΦ (c) = Φ (bz (k) ) γ(z (k),z (k+1),c)+iθ(z (k),z (k+1),c). (27) 7
8 Notice that (27) is also valid for =0. Thus, all the complex fuctios eeded for the solutio of the BVP ca be determied umerically at all poits i R C. 3 A TEST PROBLEM For a test problem, we shall apply the CVBEM proposed above to solve µ (x +2) 4 φ +(x +2) 4 2 φ =0i0<x<1, 0 <y<π, (28) x x y2 subject to φ =0 for0<x<1, (29) y=0 φ =0 for0<x<1, (30) y=π φ(0,y)=0 for0<y<π, (31) φ(1,y)=cos(y) for 0 <y<π. (32) The exact solutio of the BVP ca be obtaied by the method of separatio of variables. It is give by φ(x, y) = 27 ( 3x 3+x exp( 2x)+3exp( 2x)) exp(x +2) 2e (3e 2 2) (x +2) 3 cos(y). (33) For the PDE i (28), K(x) =(x +2) 4, f 0 (x) =1, f 1 (x) = 1 (x +2) f (x) =0 for 2, (34) ad F 0 (x) = 2 (x +2) 3, F 1(x) = 3 (x +2) 4, F (x) = 0 for 2. (35) 8
9 The CVBEM described i the previous sectio is applied to solve (28) subject to (29)-(32). The umerical values of φ are the calculated at selected iterior poits ad compared with the exact values from (33). A compariso is made i Table 1 for umerical results obtaied by puttig 10 elemets o each side of the rectagular boudary ad also 40 elemets per side (i.e. for M =40adM = 120) with the exact solutio. For M =40, the largest legth of the elemets is approximately 0.31 uits, while for M = 120 the largest legth is about uits. It is obvious from Table 1 that there is sigificat improvemet i the accuracy of the umerical results whe the umber of elemets is icreased from 40 to 120. The accuracy is i geeral good except for iterior poits extremely close to the boudary. For iterior poits very close to the boudary, it is ecessary to employ a larger umber of elemets for better accuracy. For example, at the poit (0.10, 3.10) which is extremely close to the boudary y = π, the percetage error i the calculatio of φ is about 47% if the computatio is carried out usig M =40. However, usig M = 120, we are able to reduce the percetage error to 6.1%. Further calculatio usig M = 400 brigs the percetage error dow to well uder 1%. Table 1. A compariso of the CVBEM umerical results with the exact solutio at selected iterior poits. Iterior poit CVBEM CVBEM (x, y) M =40 M =120 Exact ( ) (0.40, 0.60) (0.80, 0.70) (0.50, 0.50) (0.99, 2.10) (0.10, 3.10) SUMMARY A CVBEM is proposed for the umerical solutio of BVPs govered by a secod order elliptic PDE with variable coefficiets. A geeral solutio of 9
10 the PDE, expressible i terms of a arbitrary complex fuctio that is aalytic i the domai of iterest, is available i the literature (Clemets [3]). The proposed method makes use of the Cauchy itegral formulae to costruct approximately the relevat complex fuctio that gives the required solutio of the BVP. The task of costructig the complex fuctio requires maipulatio of the ukow data o oly the boudary of the solutio domai ad ca be evetually reduced to solvig a system of liear algebraic equatios. The method ca be efficietly implemeted o the computer, as oly the boudary eeds to be discretized ad the coefficiets of the algebraic equatios are easy to compute. It is used to solve a test problem which has a explicit exact solutio. The umerical results obtaied idicates that the proposed CVBEM works. Covergece of the umerical solutio to the exact oe is also observed whe the calculatio is refied by icreasig the umber of boudary elemets. Ackowledgemet. The first author (YSP) would like to ackowledge the assistace provided by Uiversiti Malaysia Sarawak i the form of a research grat for carryig out the work reported. Refereces [1] Brea, JJ. Iterfacial studies of refractory glass-ceramic matrix/advaced SiC fiber reiforced composite. Aual Report, ONR Cotract N C-0699, 1991, Uited Techologies Research Ceter. [2] Gibso, RE. Some results cocerig displacemets ad stresses i a ohomogeeous elastic half-space. Geotechique 1967; 17: [3] Clemets, DL. A boudary itegral equatio method for the umerical solutio of a secod order elliptic partial differetial equatio with variable coefficiets. Joural of the Australia Mathematical Society (Series B) 1980; 22:
11 [4] Ragogi, R. A solutio of Darcy s flow with variable permeability by meas of BEM ad perturbatio techique. Proceedigs of the 9th World Coferece o the Boudary Elemet Method 1987 (Edited by C. A. Brebbia) Spriger-Verlag, Berli. [5] Hromadka II, TV ad Lai, C. The Complex Variable Boudary Elemet Method i Egieerig Aalysis 1987, Spriger-Verlag, Berli. [6] Ag WT ad Park, YS. CVBEM for a system of secod-order elliptic partial differetial equatios. Egieerig Aalysis with Boudary Elemets 1998 ; 21:
Boundary Element Method (BEM)
Boudary Elemet Method BEM Zora Ilievski Wedesday 8 th Jue 006 HG 6.96 TU/e Talk Overview The idea of BEM ad its advatages The D potetial problem Numerical implemetatio Idea of BEM 3 Idea of BEM 4 Advatages
More information1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y
Questio (a) A square matrix A= A is called positive defiite if the quadratic form waw > 0 for every o-zero vector w [Note: Here (.) deotes the traspose of a matrix or a vector]. Let 0 A = 0 = show that:
More informationLecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables
Lecture 4 Coformal Mappig ad Gree s Theorem Today s topics. Solvig electrostatic problems cotiued. Why separatio of variables does t always work 3. Coformal mappig 4. Gree s theorem The failure of separatio
More informationGenerating Functions for Laguerre Type Polynomials. Group Theoretic method
It. Joural of Math. Aalysis, Vol. 4, 2010, o. 48, 257-266 Geeratig Fuctios for Laguerre Type Polyomials α of Two Variables L ( xy, ) by Usig Group Theoretic method Ajay K. Shula* ad Sriata K. Meher** *Departmet
More informationMETHOD OF FUNDAMENTAL SOLUTIONS FOR HELMHOLTZ EIGENVALUE PROBLEMS IN ELLIPTICAL DOMAINS
Please cite this article as: Staisław Kula, Method of fudametal solutios for Helmholtz eigevalue problems i elliptical domais, Scietific Research of the Istitute of Mathematics ad Computer Sciece, 009,
More informationTaylor polynomial solution of difference equation with constant coefficients via time scales calculus
TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu
More informationWe are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n
Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at
More informationPRELIM PROBLEM SOLUTIONS
PRELIM PROBLEM SOLUTIONS THE GRAD STUDENTS + KEN Cotets. Complex Aalysis Practice Problems 2. 2. Real Aalysis Practice Problems 2. 4 3. Algebra Practice Problems 2. 8. Complex Aalysis Practice Problems
More informationNumerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION
alaysia Joural of athematical Scieces 3(1): 83-93 (9) umerical Coformal appig via a Fredholm Itegral Equatio usig Fourier ethod 1 Ali Hassa ohamed urid ad Teh Yua Yig 1, Departmet of athematics, Faculty
More informationn 3 ln n n ln n is convergent by p-series for p = 2 > 1. n2 Therefore we can apply Limit Comparison Test to determine lutely convergent.
06 微甲 0-04 06-0 班期中考解答和評分標準. ( poits) Determie whether the series is absolutely coverget, coditioally coverget, or diverget. Please state the tests which you use. (a) ( poits) (b) ( poits) (c) ( poits)
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
More informationSection A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics
X0/70 NATIONAL QUALIFICATIONS 005 MONDAY, MAY.00 PM 4.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.
More informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
More informationPOWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS
Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity
More informationIndian Institute of Information Technology, Allahabad. End Semester Examination - Tentative Marking Scheme
Idia Istitute of Iformatio Techology, Allahabad Ed Semester Examiatio - Tetative Markig Scheme Course Name: Mathematics-I Course Code: SMAT3C MM: 75 Program: B.Tech st year (IT+ECE) ate of Exam:..7 ( st
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationIntroduction to Optimization Techniques. How to Solve Equations
Itroductio to Optimizatio Techiques How to Solve Equatios Iterative Methods of Optimizatio Iterative methods of optimizatio Solutio of the oliear equatios resultig form a optimizatio problem is usually
More informationDECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan
Mathematical ad Computatioal Applicatios, Vol. 9, No. 3, pp. 30-40, 04 DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS Muhammad Aslam Noor, Khalida Iayat Noor ad Asif Waheed
More informationPAPER : IIT-JAM 2010
MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure
More informationFinally, we show how to determine the moments of an impulse response based on the example of the dispersion model.
5.3 Determiatio of Momets Fially, we show how to determie the momets of a impulse respose based o the example of the dispersio model. For the dispersio model we have that E θ (θ ) curve is give by eq (4).
More informationModified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations
Noliear Aalysis ad Differetial Equatios, Vol. 5, 27, o. 4, 57-7 HIKARI Ltd, www.m-hikari.com https://doi.org/.2988/ade.27.62 Modified Decompositio Method by Adomia ad Rach for Solvig Noliear Volterra Itegro-
More informationThe axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.
5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =
More informationProblem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =
Mathematics Summer Wilso Fial Exam August 8, ANSWERS Problem 1 (a) Fid the solutio to y +x y = e x x that satisfies y() = 5 : This is already i the form we used for a first order liear differetial equatio,
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More information(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is
Calculus BC Fial Review Name: Revised 7 EXAM Date: Tuesday, May 9 Remiders:. Put ew batteries i your calculator. Make sure your calculator is i RADIAN mode.. Get a good ight s sleep. Eat breakfast. Brig:
More informationMATH 1080: Calculus of One Variable II Fall 2017 Textbook: Single Variable Calculus: Early Transcendentals, 7e, by James Stewart.
MATH 1080: Calculus of Oe Variable II Fall 2017 Textbook: Sigle Variable Calculus: Early Trascedetals, 7e, by James Stewart Uit 3 Skill Set Importat: Studets should expect test questios that require a
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationDETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.
ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics
More informationStreamfunction-Vorticity Formulation
Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More informationPAijpam.eu ON TENSOR PRODUCT DECOMPOSITION
Iteratioal Joural of Pure ad Applied Mathematics Volume 103 No 3 2015, 537-545 ISSN: 1311-8080 (prited versio); ISSN: 1314-3395 (o-lie versio) url: http://wwwijpameu doi: http://dxdoiorg/1012732/ijpamv103i314
More informationTeaching Mathematics Concepts via Computer Algebra Systems
Iteratioal Joural of Mathematics ad Statistics Ivetio (IJMSI) E-ISSN: 4767 P-ISSN: - 4759 Volume 4 Issue 7 September. 6 PP-- Teachig Mathematics Cocepts via Computer Algebra Systems Osama Ajami Rashaw,
More informationNumerical Solutions of Second Order Boundary Value Problems by Galerkin Residual Method on Using Legendre Polynomials
IOSR Joural of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 11, Issue 6 Ver. IV (Nov. - Dec. 15), PP 1-11 www.iosrjourals.org Numerical Solutios of Secod Order Boudary Value Problems
More information1 Approximating Integrals using Taylor Polynomials
Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................
More information: Transforms and Partial Differential Equations
Trasforms ad Partial Differetial Equatios 018 SUBJECT NAME : Trasforms ad Partial Differetial Equatios SUBJECT CODE : MA 6351 MATERIAL NAME : Part A questios REGULATION : R013 WEBSITE : wwwharigaeshcom
More informationμ are complex parameters. Other
A New Numerical Itegrator for the Solutio of Iitial Value Problems i Ordiary Differetial Equatios. J. Suday * ad M.R. Odekule Departmet of Mathematical Scieces, Adamawa State Uiversity, Mubi, Nigeria.
More informationExact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method
Exact Solutios for a Class of Noliear Sigular Two-Poit Boudary Value Problems: The Decompositio Method Abd Elhalim Ebaid Departmet of Mathematics, Faculty of Sciece, Tabuk Uiversity, P O Box 741, Tabuki
More informationENGI Series Page 6-01
ENGI 3425 6 Series Page 6-01 6. Series Cotets: 6.01 Sequeces; geeral term, limits, covergece 6.02 Series; summatio otatio, covergece, divergece test 6.03 Stadard Series; telescopig series, geometric series,
More informationAnalysis of composites with multiple rigid-line reinforcements by the BEM
Aalysis of composites with multiple rigid-lie reiforcemets by the BEM Piotr Fedeliski* Departmet of Stregth of Materials ad Computatioal Mechaics, Silesia Uiversity of Techology ul. Koarskiego 18A, 44-100
More informationA numerical Technique Finite Volume Method for Solving Diffusion 2D Problem
The Iteratioal Joural Of Egieerig d Sciece (IJES) Volume 4 Issue 10 Pages PP -35-41 2015 ISSN (e): 2319 1813 ISSN (p): 2319 1805 umerical Techique Fiite Volume Method for Solvig Diffusio 2D Problem 1 Mohammed
More informationNumerical Methods in Fourier Series Applications
Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets
More informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
More informationConcavity Solutions of Second-Order Differential Equations
Proceedigs of the Paista Academy of Scieces 5 (3): 4 45 (4) Copyright Paista Academy of Scieces ISSN: 377-969 (prit), 36-448 (olie) Paista Academy of Scieces Research Article Cocavity Solutios of Secod-Order
More informationSOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.
SOLUTIONS TO EXAM 3 Problem Fid the sum of the followig series 2 + ( ) 5 5 2 5 3 25 2 2 This series diverges Solutio: Note that this defies two coverget geometric series with respective radii r 2/5 < ad
More informationf(w) w z =R z a 0 a n a nz n Liouville s theorem, we see that Q is constant, which implies that P is constant, which is a contradiction.
Theorem 3.6.4. [Liouville s Theorem] Every bouded etire fuctio is costat. Proof. Let f be a etire fuctio. Suppose that there is M R such that M for ay z C. The for ay z C ad R > 0 f (z) f(w) 2πi (w z)
More informationA NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS
Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet
More informationChapter 10: Power Series
Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because
More informationPoincaré Problem for Nonlinear Elliptic Equations of Second Order in Unbounded Domains
Advaces i Pure Mathematics 23 3 72-77 http://dxdoiorg/4236/apm233a24 Published Olie Jauary 23 (http://wwwscirporg/oural/apm) Poicaré Problem for Noliear Elliptic Equatios of Secod Order i Ubouded Domais
More informationMIDTERM 3 CALCULUS 2. Monday, December 3, :15 PM to 6:45 PM. Name PRACTICE EXAM SOLUTIONS
MIDTERM 3 CALCULUS MATH 300 FALL 08 Moday, December 3, 08 5:5 PM to 6:45 PM Name PRACTICE EXAM S Please aswer all of the questios, ad show your work. You must explai your aswers to get credit. You will
More informationSolution of Differential Equation from the Transform Technique
Iteratioal Joural of Computatioal Sciece ad Mathematics ISSN 0974-3189 Volume 3, Number 1 (2011), pp 121-125 Iteratioal Research Publicatio House http://wwwirphousecom Solutio of Differetial Equatio from
More informationCastiel, Supernatural, Season 6, Episode 18
13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio
More informationFinite Difference Approximation for First- Order Hyperbolic Partial Differential Equation Arising in Neuronal Variability with Shifts
Iteratioal Joural of Scietific Egieerig ad Research (IJSER) wwwiseri ISSN (Olie): 347-3878, Impact Factor (4): 35 Fiite Differece Approimatio for First- Order Hyperbolic Partial Differetial Equatio Arisig
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More information[ 11 ] z of degree 2 as both degree 2 each. The degree of a polynomial in n variables is the maximum of the degrees of its terms.
[ 11 ] 1 1.1 Polyomial Fuctios 1 Algebra Ay fuctio f ( x) ax a1x... a1x a0 is a polyomial fuctio if ai ( i 0,1,,,..., ) is a costat which belogs to the set of real umbers ad the idices,, 1,...,1 are atural
More informationThe Analysis of the Non-linear Deflection of Non-straight Ludwick type Beams Using Lie Symmetry Groups
Proceedigs of the 3 rd Iteratioal Coferece o Cotrol, Dyamic Systems, ad Robotics CDSR 16 Ottawa, Caada May 9 10, 016 Paper No. 107 DOI: 10.11159/cdsr16.107 The Aalysis of the No-liear Deflectio of No-straight
More informationEFFECTIVE WLLN, SLLN, AND CLT IN STATISTICAL MODELS
EFFECTIVE WLLN, SLLN, AND CLT IN STATISTICAL MODELS Ryszard Zieliński Ist Math Polish Acad Sc POBox 21, 00-956 Warszawa 10, Polad e-mail: rziel@impagovpl ABSTRACT Weak laws of large umbers (W LLN), strog
More informationwavelet collocation method for solving integro-differential equation.
IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 3 (arch. 5), V3 PP -7 www.iosrje.org wavelet collocatio method for solvig itegro-differetial equatio. Asmaa Abdalelah Abdalrehma
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam will cover.-.9. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for you
More informationFinite Difference Approximation for Transport Equation with Shifts Arising in Neuronal Variability
Iteratioal Joural of Sciece ad Research (IJSR) ISSN (Olie): 39-764 Ide Copericus Value (3): 64 Impact Factor (3): 4438 Fiite Differece Approimatio for Trasport Equatio with Shifts Arisig i Neuroal Variability
More informationPrinciple Of Superposition
ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give
More informationIIT JAM Mathematical Statistics (MS) 2006 SECTION A
IIT JAM Mathematical Statistics (MS) 6 SECTION A. If a > for ad lim a / L >, the which of the followig series is ot coverget? (a) (b) (c) (d) (d) = = a = a = a a + / a lim a a / + = lim a / a / + = lim
More informationb i u x i U a i j u x i u x j
M ath 5 2 7 Fall 2 0 0 9 L ecture 1 9 N ov. 1 6, 2 0 0 9 ) S ecod- Order Elliptic Equatios: Weak S olutios 1. Defiitios. I this ad the followig two lectures we will study the boudary value problem Here
More informationMonte Carlo Integration
Mote Carlo Itegratio I these otes we first review basic umerical itegratio methods (usig Riema approximatio ad the trapezoidal rule) ad their limitatios for evaluatig multidimesioal itegrals. Next we itroduce
More information(a) (b) All real numbers. (c) All real numbers. (d) None. to show the. (a) 3. (b) [ 7, 1) (c) ( 7, 1) (d) At x = 7. (a) (b)
Chapter 0 Review 597. E; a ( + )( + ) + + S S + S + + + + + + S lim + l. D; a diverges by the Itegral l k Test sice d lim [(l ) ], so k l ( ) does ot coverge absolutely. But it coverges by the Alteratig
More informationNew Version of the Rayleigh Schrödinger Perturbation Theory: Examples
New Versio of the Rayleigh Schrödiger Perturbatio Theory: Examples MILOŠ KALHOUS, 1 L. SKÁLA, 1 J. ZAMASTIL, 1 J. ČÍŽEK 2 1 Charles Uiversity, Faculty of Mathematics Physics, Ke Karlovu 3, 12116 Prague
More information5. Matrix exponentials and Von Neumann s theorem The matrix exponential. For an n n matrix X we define
5. Matrix expoetials ad Vo Neuma s theorem 5.1. The matrix expoetial. For a matrix X we defie e X = exp X = I + X + X2 2! +... = 0 X!. We assume that the etries are complex so that exp is well defied o
More informationThe Numerical Solution of Singular Fredholm Integral Equations of the Second Kind
WDS' Proceedigs of Cotributed Papers, Part I, 57 64, 2. ISBN 978-8-7378-39-2 MATFYZPRESS The Numerical Solutio of Sigular Fredholm Itegral Equatios of the Secod Kid J. Rak Charles Uiversity, Faculty of
More informationSome properties of Boubaker polynomials and applications
Some properties of Boubaker polyomials ad applicatios Gradimir V. Milovaović ad Duša Joksimović Citatio: AIP Cof. Proc. 179, 1050 (2012); doi: 10.1063/1.756326 View olie: http://dx.doi.org/10.1063/1.756326
More informationMachine Learning for Data Science (CS 4786)
Machie Learig for Data Sciece CS 4786) Lecture & 3: Pricipal Compoet Aalysis The text i black outlies high level ideas. The text i blue provides simple mathematical details to derive or get to the algorithm
More information*X203/701* X203/701. APPLIED MATHEMATICS ADVANCED HIGHER Numerical Analysis. Read carefully
X0/70 NATIONAL QUALIFICATIONS 006 MONDAY, MAY.00 PM.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.
More informationComplex Analysis Spring 2001 Homework I Solution
Complex Aalysis Sprig 2001 Homework I Solutio 1. Coway, Chapter 1, sectio 3, problem 3. Describe the set of poits satisfyig the equatio z a z + a = 2c, where c > 0 ad a R. To begi, we see from the triagle
More informationNEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE
UPB Sci Bull, Series A, Vol 79, Iss, 207 ISSN 22-7027 NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE Gabriel Bercu We itroduce two ew sequeces of Euler-Mascheroi type which have fast covergece
More informationChapter 7 Isoperimetric problem
Chapter 7 Isoperimetric problem Recall that the isoperimetric problem (see the itroductio its coectio with ido s proble) is oe of the most classical problem of a shape optimizatio. It ca be formulated
More informationNTMSCI 5, No. 1, (2017) 26
NTMSCI 5, No. 1, - (17) New Treds i Mathematical Scieces http://dx.doi.org/1.85/tmsci.17.1 The geeralized successive approximatio ad Padé approximats method for solvig a elasticity problem of based o the
More informationThe Method of Least Squares. To understand least squares fitting of data.
The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve
More informationChapter 2 Feedback Control Theory Continued
Chapter Feedback Cotrol Theor Cotiued. Itroductio I the previous chapter, the respose characteristic of simple first ad secod order trasfer fuctios were studied. It was show that first order trasfer fuctio,
More informationA PROOF OF THE TWIN PRIME CONJECTURE AND OTHER POSSIBLE APPLICATIONS
A PROOF OF THE TWI PRIME COJECTURE AD OTHER POSSIBLE APPLICATIOS by PAUL S. BRUCKMA 38 Frot Street, #3 aaimo, BC V9R B8 (Caada) e-mail : pbruckma@hotmail.com ABSTRACT : A elemetary proof of the Twi Prime
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam 4 will cover.-., 0. ad 0.. Note that eve though. was tested i exam, questios from that sectios may also be o this exam. For practice problems o., refer to the last review. This
More informationA NOTE ON BOUNDARY BLOW-UP PROBLEM OF u = u p
A NOTE ON BOUNDARY BLOW-UP PROBLEM OF u = u p SEICK KIM Abstract. Assume that Ω is a bouded domai i R with 2. We study positive solutios to the problem, u = u p i Ω, u(x) as x Ω, where p > 1. Such solutios
More informationA collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation
Iteratioal Joural of Mathematics Research. ISSN 0976-5840 Volume 9 Number 1 (017) pp. 45-51 Iteratioal Research Publicatio House http://www.irphouse.com A collocatio method for sigular itegral equatios
More informationAssignment 1 : Real Numbers, Sequences. for n 1. Show that (x n ) converges. Further, by observing that x n+2 + x n+1
Assigmet : Real Numbers, Sequeces. Let A be a o-empty subset of R ad α R. Show that α = supa if ad oly if α is ot a upper boud of A but α + is a upper boud of A for every N. 2. Let y (, ) ad x (, ). Evaluate
More informationRotationally invariant integrals of arbitrary dimensions
September 1, 14 Rotatioally ivariat itegrals of arbitrary dimesios James D. Wells Physics Departmet, Uiversity of Michiga, A Arbor Abstract: I this ote itegrals over spherical volumes with rotatioally
More informationU8L1: Sec Equations of Lines in R 2
MCVU U8L: Sec. 8.9. Equatios of Lies i R Review of Equatios of a Straight Lie (-D) Cosider the lie passig through A (-,) with slope, as show i the diagram below. I poit slope form, the equatio of the lie
More informationAn Analysis of a Certain Linear First Order. Partial Differential Equation + f ( x, by Means of Topology
Iteratioal Mathematical Forum 2 2007 o. 66 3241-3267 A Aalysis of a Certai Liear First Order Partial Differetial Equatio + f ( x y) = 0 z x by Meas of Topology z y T. Oepomo Sciece Egieerig ad Mathematics
More informationTMA4205 Numerical Linear Algebra. The Poisson problem in R 2 : diagonalization methods
TMA4205 Numerical Liear Algebra The Poisso problem i R 2 : diagoalizatio methods September 3, 2007 c Eiar M Røquist Departmet of Mathematical Scieces NTNU, N-749 Trodheim, Norway All rights reserved A
More informationDiscrete Orthogonal Moment Features Using Chebyshev Polynomials
Discrete Orthogoal Momet Features Usig Chebyshev Polyomials R. Mukuda, 1 S.H.Og ad P.A. Lee 3 1 Faculty of Iformatio Sciece ad Techology, Multimedia Uiversity 75450 Malacca, Malaysia. Istitute of Mathematical
More informationChapter 2 The Solution of Numerical Algebraic and Transcendental Equations
Chapter The Solutio of Numerical Algebraic ad Trascedetal Equatios Itroductio I this chapter we shall discuss some umerical methods for solvig algebraic ad trascedetal equatios. The equatio f( is said
More informationA widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α
Nice plottig of proteis: I A widely used display of protei shapes is based o the coordiates of the alpha carbos - - C α -s. The coordiates of the C α -s are coected by a cotiuous curve that roughly follows
More information8. Applications To Linear Differential Equations
8. Applicatios To Liear Differetial Equatios 8.. Itroductio 8.. Review Of Results Cocerig Liear Differetial Equatios Of First Ad Secod Orders 8.3. Eercises 8.4. Liear Differetial Equatios Of Order N 8.5.
More informationThe picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled
1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how
More informationSolutions to quizzes Math Spring 2007
to quizzes Math 4- Sprig 7 Name: Sectio:. Quiz a) x + x dx b) l x dx a) x + dx x x / + x / dx (/3)x 3/ + x / + c. b) Set u l x, dv dx. The du /x ad v x. By Itegratio by Parts, x(/x)dx x l x x + c. l x
More informationStopping oscillations of a simple harmonic oscillator using an impulse force
It. J. Adv. Appl. Math. ad Mech. 5() (207) 6 (ISSN: 2347-2529) IJAAMM Joural homepage: www.ijaamm.com Iteratioal Joural of Advaces i Applied Mathematics ad Mechaics Stoppig oscillatios of a simple harmoic
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationPhysics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.
Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;
More informationON SINGULAR INTEGRAL OPERATORS
ON SINGULAR INTEGRAL OPERATORS DEJENIE ALEMAYEHU LAKEW Abstract. I this paper we study sigular itegral operators which are hyper or weak over Lipscitz/Hölder spaces ad over weighted Sobolev spaces de ed
More informationUNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014
UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 6C Problem Set 4 Bejami Stahl November 6, 4 BOAS, P. 63, PROBLEM.-5 The Laguerre differetial equatio, x y + ( xy + py =, will be solved
More informationLOWER BOUNDS FOR THE BLOW-UP TIME OF NONLINEAR PARABOLIC PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
Electroic Joural of Differetial Equatios, Vol. 214 214), No. 113, pp. 1 5. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.ut.edu ftp ejde.math.txstate.edu LOWER BOUNDS FOR THE BLOW-UP
More informationFourier Series and their Applications
Fourier Series ad their Applicatios The fuctios, cos x, si x, cos x, si x, are orthogoal over (, ). m cos mx cos xdx = m = m = = cos mx si xdx = for all m, { m si mx si xdx = m = I fact the fuctios satisfy
More informationCHAPTER 1 SEQUENCES AND INFINITE SERIES
CHAPTER SEQUENCES AND INFINITE SERIES SEQUENCES AND INFINITE SERIES (0 meetigs) Sequeces ad limit of a sequece Mootoic ad bouded sequece Ifiite series of costat terms Ifiite series of positive terms Alteratig
More information