arxiv: v1 [math.na] 21 Nov 2012
|
|
- Nora Owen
- 6 years ago
- Views:
Transcription
1 The Fourier Transforms of the Chebyshev and Legendre Polynomials arxiv: v1 [math.na] 1 Nov 1 A S Fokas School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 138, USA Permanent address: Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge CB3 WA S A Smitheman School of Biosciences, University of Nottingham, Sutton Bonington Campus, Loughborough, Leicestershire LE1 5RD, UK September 1 Abstract Analytic expressions for the Fourier transforms of the Chebyshev and Legendre polynomials are derived, and the latter is used to find a new representation for the half-order Bessel functions. The numerical implementation of the so-called unified method in the interior of a convex polygon provides an example of the applicability of these analytic expressions. 1 Introduction The importance of Chebyshev and Legendre polynomials in numerical analysis is well known. The finite Fourier transform plays a crucial role in mathematics and applications. In spite of these facts it appears that the finite Fourier transform of Chebyshev and Legendre polynomials has not been constructed. Here, we present explicit formulae for these transforms. The finite Fourier transform of the Legendre polynomials can be expressed in terms of the half-order Bessel functions, thus a direct application of our results is the construction of an explicit representation for the half-order Bessel functions. Our results are motivated from the recent work of B Fornberg and collaborators [1], [6] on the numerical implementation of the unified method of []-[4] to linear elliptic partial differential equations PDEs formulated in the interior of a convex polygon. In this case, the unified method yields a simple algebraic equation, the so-called global relation, which couples the finite Fourier transforms of the given boundary data with the finite Fourier transforms of the unknown boundary values. For the determination of these boundary values one has to choose a appropriate basis functions and b suitable collocation 1
2 points in the Fourier space. Several such choices have appeared in the literature [1], [5]-[1]; it appears that the best choice is a the unknown boundary values are expanded in terms of Legendre functions and b the collocation points are chosen on on rays introduced in [5], [1]. We note that Chebyshev and Legendre functions give rise to equivalent finite bases of L, and hence either choice will result in the same numerical method; however, in general the conditioning of the resulting linear systems will differ as conditioning is not invariant under matrix column operations. As an example of the applicability of these analytic formulae we apply this approach to the simplest possible polygon, namely to a square, and we use the explicit construction of the finite Fourier transform of the Legendre functions instead of the usual representation in terms of the Bessel functions. Fourier transforms of the Chebyshev and Legendre polynomials Theorem 1 Let ˆT m denote the finite Fourier transform of the Chebyshev polynomial T m x, i.e. 1 ˆT m = e x T m xdx, C, m=,1,,...,.1 1 where T m x denotes the Chebyshev polynomial Then, T m x=cosmcos 1 x, 1<x<1, m=,1,,...., m=1, ˆT m = 1 m+1 1 m, m=,,3, Furthermore, [ ] m+1 ˆT m = αn m e e + 1n+m n=1 n n, C {}, m=,1,,...,.4 where the coefficients αn m are defined as follows: α m 1 = 1 m, α m = 1 m+1 m,.5 α m n = 1m+n 1 n m m n+ n+k 3 n+k 3 m j, n=3,4,...,m+1..6
3 Proof Making in.1 the change of variables x = cos w, we find π ˆT m = e cosw sinwcosmwdw, C, m=,1,...,.7 from which the results for = follow. Using the identity sinwe cosw = 1 d dw and employing integration by parts in.7 we find e cosw,.8 ˆT m = 1 [ ] e 1 m e + mk m, C {}, m=,1,...,.9 where K m z= π e zcosw sinmwdw, z C {}, m=,1,....1 We now derive an analytic expression for K m z. It is clear that K z=, and on replacing by z in.8 we find K 1 z= 1 z ezcosw π w= The function K m z satisfies the relation = 1 z ez e z, z C {}..11 K m+1 z= m z K mz+k z+ z ez + 1 e z, For the derivation of this equation we will use the identity Thus, z C {}, m=1,,....1 sinm+1w=sinw+sinwcosmw..13 π K m+1 z= e zcosw sinm+1wdw = K z+ π sinwe zcosw cosmwdw..14 Using integration by parts we find that the above integral can be rewritten in the following form: π sinwe zcosw cosmwdw = 1 π z ezcosw cosmw m z π = 1 z ez + 1 e z m z K mz. e zcosw sinmwdw.15 3
4 We now show that for m=,1,..., K m z = m z ez + 1 m e z + 1 z n=1 z n 1 n e z + 1 m e z m n n+ n+ m j,.16 from which.4 follows using.9. Clearly, this relation is valid for m= and m=1. Furthermore, the function K m z in.16 satisfies the relation.1 if, and only if, 1 z m n=1 z n 1 n e z + 1 m+1 e z m n+1 n+ n+ = m z ez + 1 m e z m n z 1 n e z + 1 m e z m n n+ n=1 z + 1 m n 1 n e z + 1 e z n+ m n 1 z z n=1 n+ n+ m+1 j m j j..17 To prove.17, we compare coefficients of each power of z on each side of this equation. The coefficient of z on the left-hand side of.17 is given by e z + 1 m+1 e z m = e z + 1 m e z km+1 k [mm+1 13 mm+1m+1 ],.18 the coefficient of z on the right-hand side of.17 is given by [ ] e z + 1 m e z m + m k k = e z + 1 m e z [m + m 13 ] m m 3,.19 and the expressions in.18,.19 are equivalent. The coefficient of z m+1 on the left-hand side of.17 is given by m 1 m e z + 1 m+1 e z m m+1 j= m 1 m e z + 1 m+1 e z m!,. 4
5 the coefficient of z m+1 on the right-hand side of.17 is given by m e z + 1 m e z m j = m m e z + 1 m+1 e z!,.1 and the expressions in.,.1 are equivalent. The coefficient of z m on the left-hand side of.17 is given by 1 e z + 1 m+1 e z = 1 e z + e z = m m! 1 e z + e z, [ m+k m+1 j+m m+k m j= m+1 j m+1 j the coefficient of z m on the right-hand side of.17 is given by ]. m e z + 1 m e z = e z + e z [ m = m e z + e z!, m+k 3 m j+ m+k 3 j= m j m j ].3 and the expressions in.,.3 are equivalent. Finally, we consider z n+1, n=, 3,..., m. Omitting the common factor of n 1 n e z + 1 e z for convenience, the difference of the coefficients of the two sides of equation.17 is given by m n+1 [ m Furthermore, n+k n+k n+k m n+ = n+k m j [ n+ m j m n 1 n+ + n+ n+ mn n+ m k+1 n+ n+ m+1 j ] n+k m+1 j ] j..4 m j,.5 5
6 and m n 1 = = n+ m n+1 n+k 3 k 3 k=3 m n+1 n+ n+ j n+k 3 j k m k+1 n+k n+n+k m j..6 Hence the above difference is given by n+ m n+1 = n = n where a k = [ mn n+ m k+1 k 1 n+n+k n+ 1 n+ m n+1 m n [ m n+k 3m k+1 n+k [ n+k k ] n+k kn+k n+ n+kn+ +1 m j m j ] = ] ] m j ] n+k m n 1 n 1! a k,.7 m!n+k!n+k, k= 1,,..., m n..8 m n k!! Hence a m n+1 k = a k. If m n is odd, then m n+1 n+ =,.9 and hence a m n+1 =. Thus the coefficients of z n+1 of the two sides of.17 are equal for n=, 3,..., m. This establishes the validity of.17, which completes the proof of.4. 6
7 Theorem Let ˆP m denote the finite Fourier transform of the Legendre polynomial P m x, i.e. 1 ˆP m = e x P m xdx, C, m=,1,,...,.3 1 where P m x denotes the Legendre polynomial Then, d m P m x= 1 m m! dx mx 1 m, 1<x<1, m=,1,, {, m=, ˆP m =, m=1,,....3 Furthermore, [ ] m+1 ˆP m = βn m e e + 1m+n n=1 n n, C {}, m=,1,,...,.33 where the coefficients βn m are defined as follows: β m 1 = 1 if m even,.34a m+n 3 βn m = m+n n 1 m+n βn m = 1 n 1 m+n 3 + R m n+3 n, m+n 3 m, { 3,5,...,m+1 if m even, n=,4,...,m+1 if m odd,.34b {,4,...,m if m even, R m n +1, m+n 1m, n= 1,3,...,m if m odd,.34c where t s! R s,t r= r k+1, s t, s, s t, and = k=s s t!t! t 1, s> t,, s< t..35 Proof We first consider.3 for m=: Using the identity 1 ˆP = P xdx= m+1p m x= d dx P m+1x P x, 1<x<1, m=1,,...,.37 as well as the equations P m 1=1 and P m 1= 1 m for m =,1,... we find the second line of.3: 1 ˆP m = P m xdx= 1 1 m+1 [P m+1x P x] 1 1 =..38 7
8 For m=, 1, P m =T m and thus ˆP m = ˆT m ; hence.4 establishes the validity of.33. Furthermore, using.37 and integrating by parts we find ˆP m+1 = i m+1 ˆP m + ˆP, C {}, m=1,, For convenience, we consider the following equivalent forms of.33 for m even and odd respectively: ˆP m = 1 m m 1 1 j m j+1 ψ m j 1 m j + χ R m j 1 j+1,m j 1 m + ψ+ 1 m m χr 1,m,.4 and ˆP m = i where 1+ 1 i 1 m 3 R 1, m χ im+ ψ 1 j m j+1 ψ m j 1 m j + χ R m j 1 j+1,m j 1 m,.41 ψ= i e e and χ=ψ e + e..4 8
9 For odd m,.39 is equivalent to the following expression: 1 m+1 +ψ+ = 1 + ψ+ 1 j m j R j+1,m jm+1 m j+3 m+1 1 m+1 m 3 R 1,m m+1χ m j m j+ 1 m j ψ+ χ m j 1 j m j R j+1,m j m j m j m j 1 ψ+ χ m j 1 m j 1 1 m+1 R 1,m χ 1 1+ R 1, m χ 1 m+1m+ψ 1 m 3 Furthermore, m 3 = m+1 1 j m j 1 R j+1,m j 1m χ+m j+1 m j 1 1 j m j 1 j 1 j= m j ψ. m j m j m j 1 ψ m j 1 m j + χ R m j j+1,m j m j+1 ψ m j+ 1 + χ R m j j,m j
10 It is straightforward to show that R j+1,m j m+1=m j+1r j+1,m j 1 m, R j,m j 1 = j+ 1R j+1,m j 1 m..45 Hence for odd m,.39 is equivalent to the following equation: 1 m+1 + = m+1 1 m+1 m j m+1 m jm j+1r j+1,m j 1m m j m j+3 ψ+ m j+ 1 m+1r 1, mχ 1 j j= m j j+ 1R j+1,m j 1m m j+1 ψ+ m j+ 1 R,m mχ 1 1+ R 1, m χ 1 m+1m+ψ + m+1 1 m+1 m 3 1 j m j R j+1,m j 1m m j 1 χ+m j+1 m j χ m j χ m j ψ. m j.46 1
11 Hence, 1 m+1 m 3 j= + 1 j m j R j+1,m j 1m m j+1 1 m j m j+3 j+ 1 m j+ 1 m j+ 1 m+1 ψ m j m j m+1 χ m j 1 m j + m j+1 j+ 1 m j R,m m[m+1ψ+χ m+1m χ+ψ] 1 1 m+ 8 m+4m+1ψ+ 1 m+1χ 1 m 8 m+m 3ψ+ 1 1 χ m+1 + m+1r 1, m+ 3 1 m+1 R,m m m m m+1m+ψ =. Equation.47 follows from the identities 1+ 1 m j m j+3 j+ 1 m j+ 1 m j+ 1 m j m j+1 j+ 1 m j m+1 m j m+1 R 1, m m j m j 1 χ.47 =,.48 =,.49 11
12 m+1 1 and the coefficients of m+1 R 1, mχ, 1 ψr,m m, 1 χr,m m, 1 1 ψ and χ on the left hand side 8 of.47 vanishing: m+1 m+1=, m+1 m+1=, m+1m 3=, m+m+1m+4 mm 3 4m+1=, m+1m+ m m+1=. The validity of.39 for even m can be demonstrated similarly, completing the proof of equation An explicit representation for the half-order Bessel functions Lemma 3 Let J m+ 1 denote the half-order Bessel function, i.e. J m+ 1 = 1 π π + 1m+1 π cos m+ 1 τ sinτ dτ exp sinhτ m+ 1 τ dτ, C, m=,1,, Then J m+ 1 =, m=,1,,..., 3. and J m+ 1 admits the following explicit representation: J m+ 1= 1 m+1 π βn m n=1 [ e e + 1 n+m i n m n 1 i n m n 1 C {}, m=,1,,..., 3.3 ], where the coefficients β m n are defined by equations.34. Proof Recall that the finite Fourier transforms of the Legendre polynomials ˆP m can be expressed in the form ˆP m = 1 π i m J m This equation, together with the results of Theorem, immediately gives equation
13 y S 4 S 1 S 3 x S Figure 1: The square with sides of length. 4 The modified Helmholtz equation in the interior of a square Let ux, y satisfy the modified Helmholtz equation u x + u 4u=, 4.1 y in the interior of a convex polygon. Associated with each side {S j } n 1 of this n-gon, define the following functions {[ û j = e z z u y ] [ u dx+ u x ] } u dy, S j j = 1,,...,n, C {}, 4. and {[ ũ j = e z z u y ] [ u dx+ u x + 1 ] } u dy, S j j = 1,,...,n, C {}. 4.3 Then, the functions{û j } n 1 and{ũ j} n 1 satisfy the global relations [1] n û j =, n ũ j =, C {}. 4.4 For brevity of presentation, we consider the simplest possible polygon, namely the square with corners at Figure 1 1,1, 1,1, 1, 1, 1,1. For the sides S 1, S, S 3, S 4, we have respectively z= 1+iy, z=x i, z=1+iy, z=x+i
14 Hence, taking into account the orientations of the sides, we find the following expressions: û 1 =e + [ 1 1 e + y 1 u 1 x ] u 1 dy, 4.6a +1 û =e + [ 1 +1 e x 1 u y ] u dx, 4.6b 1 û 3 =e + [ 1 +1 e + y 1 u 3 x ] u 3 dy, 4.6c 1 û 4 =e + [ 1 1 e x 1 u 4 y ] u 4 dx. 4.6d +1 Let ˆD j and ˆN j denote the parts of û j corresponding to the Dirichlet and Neumann boundary values. Then, û 1 = e + 1 ˆN e + 1 ˆD 1, 4.7a û = e + 1 ˆN e + 1 ˆD, 4.7b û 3 =e + 1 ˆN e + 1 ˆD 3, 4.7c û 4 =e + 1 ˆN e + 1 ˆD d For simplicity, we consider the symmetric Dirichlet boundary value problem: Then, Thus, u 1 = u 1,y= cosh1cosh 3y+cosh 3coshy, 1<y<1; 4.8a u 3 = u1,y = cosh1cosh 3y+cosh 3coshy, 1<y<1; 4.8b u = ux, 1= cosh1cosh 3x+cosh 3coshx, 1<x<1; 4.8c u 4 = ux,1 = cosh1cosh 3x+cosh 3coshx, 1<x<1. 4.8d ux,y=u x,y, ux,y=ux, y, ux,y=uy,x. 4.9 u 3 x = u 1 x, u 4 y = u y, u y = u 1 x. 4.1 x,y y,x Hence, the first of the global relations 4.4 becomes cos 1 ˆN 1 +cos 1 ˆN 1 = 1 sin 1 ˆD sin 1 ˆD 1, C {}. 4.11a 14
15 The simplest way to obtain the second of the global relations 4.4 is to take the Schwartz conjugate of equation 4.11a i.e. to take the complex conjugate of 4.11a and then to replace with. This yields the equation cos = 1 1 ˆN 1 +cos 1 sin 1 ˆD 1 + ˆN 1 1 sin 1 ˆD 1, C {}. 4.11b Using N basis functions to approximate u 1 x, equations 4.11 yield equations for N unknowns. Regarding the numerical solution of the global relations 4.11, Fourier basis functions [5], [7] [1], as well as Chebyshev and Legendre polynomials [1], [6], [9], have been used to approximate u 1 x. In most of the earlier papers the collocation points C {} were chosen to lie on the rays in the complex -plane which are parallel to the edges of the polygon and its reflection in the imaginary axis. Recently B Fornberg and collaborators introduced the use of the so-called Halton nodes, [1], [6]. It appears that the most efficient numerical method involves a approximating the data u 1 x in terms of Legendre polynomials following Fornberg and b using the collocation points employed in our earlier work [1] where ideas of A G Sifalakis and collaborators [9] for the Laplace equation were extended to the modified Helmholtz equation. Furthermore, in order to ensure that the collocation matrix remains well conditioned as the number N of basis functions increases, it is important following Fornberg to i divide each row, as well as each column, of the collocation matrix by its l 1 -norm [1] and ii to over-determine the linear system by choosing the number of collocation points to be about the same as the number of unknowns. Numerical experiments suggest that Legendre polynomials yield spectral accuracy rather than the algebraic accuracy found with a Fourier basis. Furthermore, as a result of choosing the collocation points to be on the above rays, the semi-block circulant structure of the collocation matrix for regular polygons, demonstrated for the Laplace equation in [8], is preserved in modified Helmholtz equation as well. Plots of the relative error E defined in [1], as well as of the matrix condition number as a function of N, for N/, N, 3N/ and N collocation points, are presented in Figure. The rectangular collocation matrix was inverted by using the backslash command in Matlab. It is clear that over-determining the linear system by a factor of is sufficient to achieve very good matrix conditioning. 5 Acknowledgements ASF is grateful to EPSRC, UK and to Onassis foundation, USA for partial support. SAS wishes to thank Prof B Fornberg for many useful discussions. 15
16 relative error, E # coll pts = N/ # coll pts = N # coll pts = 3N/ # coll pts = N condition number of collocation matrix # coll pts = N/ # coll pts = N # coll pts = 3N/ # coll pts = N no of basis functions, N no of basis functions, N Figure : Numerical solution of the global relations 4.11 for the symmetric Dirichlet boundary value problem 4.8. References [1] Davis, C.-I.: Numerical tests of the Fokas method for Helmholtz-type partial differential equations: Dirichlet to Neumann Maps, Masters Thesis, University of Colorado at Boulder 8 [] Fokas, A.S.: A unified transform method for solving linear and certain nonlinear PDEs. Proc. R. Soc. A 453, [3] Fokas, A.S.: On the integrability of linear and nonlinear partial differential equations. J. Math. Phys. 41, [4] Fokas, A.S.: Two-dimensional linear PDEs in a convex polygon. Proc. R. Soc. Lond. Ser. A 457, [5] Fokas, A.S., Flyer, N., Smitheman S.A., Spence E.A.: A semi-analytical numerical method for solving evolution and elliptic partial differential equations. J. Comput. Appl. Math. 7, [6] Fornberg, B., Flyer, N.: A numerical implementation of Fokas boundary integral approach: Laplace s equation on a polygonal domain. Proc. R. Soc. A 467, [7] Fulton, S.R., Fokas, A.S., Xenophontos, C.A.: An analytical method for linear elliptic PDEs and its numerical implementation. J. Comput. Appl. Math. 167, [8] Saridakis, Y.G., Sifalakis, A.G., Papadopoulou, E.P.: Efficient numerical solution of the generalized Dirichlet-Neumann map for linear elliptic PDEs in regular polygon domains. J. Comput. Appl. Math. 36, [9] Sifalakis, A.G., Fokas, A.S., Fulton, S.R., Saridakis, Y.G.: The generalized Dirichlet-Neumann map for linear elliptic PDEs and its numerical implementation. J. Comput. Appl. Math. 19,
17 [1] Sifalakis, A.G., Fulton, S.R., Papadopoulou, E.P., Saridakis, Y.G.: Direct and iterative solution of the generalized Dirichlet-Neumann map for elliptic PDEs on square domains. J. Comput. Appl. Math. 7, [11] Sifalakis, A.G., Papadopoulou, E.P., Saridakis, Y.G.: Numerical study of iterative methods for the solution of the Dirichlet-Neumann map for linear elliptic PDEs on regular polygon domains. Int. J. Appl. Math. Comput. Sci. 4, [1] Smitheman, S.A., Spence, E.A., Fokas, A.S.: A spectral collocation method for the Laplace and modified Helmholtz equations in a convex polygon. IMA J. Numer. Anal. 3,
arxiv: 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 informationBoundary value problems and medical imaging
Journal of Physics: Conference Series OPEN ACCESS Boundary value problems and medical imaging To cite this article: Athanasios S Fokas and George A Kastis 214 J. Phys.: Conf. Ser. 49 1217 View the article
More informationDirect and Iterative Solution of the Generalized Dirichlet-Neumann Map for Elliptic PDEs on Square Domains
Direct and Iterative Solution of the Generalized Dirichlet-Neumann Map for Elliptic PDEs on Square Domains A.G.Sifalakis 1, S.R.Fulton, E.P. Papadopoulou 1 and Y.G.Saridakis 1, 1 Applied Mathematics and
More informationThe elliptic sinh-gordon equation in the half plane
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 8 25), 63 73 Research Article The elliptic sinh-gordon equation in the half plane Guenbo Hwang Department of Mathematics, Daegu University, Gyeongsan
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationWell-posedness of boundary-value problems for the linear Benjamin-Bona-Mahony equation
Well-posedness of boundary-value problems for the linear Benjamin-Bona-Mahony equation Vishal Vasan and Bernard Deconinck Department of Applied Mathematics University of Washington Seattle WA 98195 July
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationNonlinear Integral Equation Formulation of Orthogonal Polynomials
Nonlinear Integral Equation Formulation of Orthogonal Polynomials Eli Ben-Naim Theory Division, Los Alamos National Laboratory with: Carl Bender (Washington University, St. Louis) C.M. Bender and E. Ben-Naim,
More informationDouble Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids
Double Layer Potentials on Polygons and Pseudodifferential Operators on Lie Groupoids joint work with Hengguang Li Yu Qiao School of Mathematics and Information Science Shaanxi Normal University Xi an,
More informationA FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS
Proceedings of ALGORITMY 2005 pp. 222 229 A FAST SOLVER FOR ELLIPTIC EQUATIONS WITH HARMONIC COEFFICIENT APPROXIMATIONS ELENA BRAVERMAN, MOSHE ISRAELI, AND ALEXANDER SHERMAN Abstract. Based on a fast subtractional
More informationA Generalization of both the Method of Images and of the Classical Integral Transforms
A Generalization of both the Method of Images and of the Classical Integral Transforms Athanassios S. Fokas 1 and Daniel ben-avraham 2 1 Institute for Nonlinear Studies, Clarkson University, Potsdam, NY
More informationA direct solver for elliptic PDEs in three dimensions based on hierarchical merging of Poincaré-Steklov operators
(1) A direct solver for elliptic PDEs in three dimensions based on hierarchical merging of Poincaré-Steklov operators S. Hao 1, P.G. Martinsson 2 Abstract: A numerical method for variable coefficient elliptic
More informationSurface Tension Effect on a Two Dimensional. Channel Flow against an Inclined Wall
Applied Mathematical Sciences, Vol. 1, 007, no. 7, 313-36 Surface Tension Effect on a Two Dimensional Channel Flow against an Inclined Wall A. Merzougui *, H. Mekias ** and F. Guechi ** * Département de
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
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 informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationThe method of successive approximations for exact solutions of Laplace equation and of heat-like and wave-like equations with variable coefficients
The method of successive approximations for exact solutions of Laplace equation and of heat-like and wave-like equations with variable coefficients T. Zhanlav and D. Khongorzul National University of Mongolia,
More informationA fast and well-conditioned spectral method: The US method
A fast and well-conditioned spectral method: The US method Alex Townsend University of Oxford Leslie Fox Prize, 24th of June 23 Partially based on: S. Olver & T., A fast and well-conditioned spectral method,
More informationPDEs, Homework #3 Solutions. 1. Use Hölder s inequality to show that the solution of the heat equation
PDEs, Homework #3 Solutions. Use Hölder s inequality to show that the solution of the heat equation u t = ku xx, u(x, = φ(x (HE goes to zero as t, if φ is continuous and bounded with φ L p for some p.
More informationWRT in 2D: Poisson Example
WRT in 2D: Poisson Example Consider 2 u f on [, L x [, L y with u. WRT: For all v X N, find u X N a(v, u) such that v u dv v f dv. Follows from strong form plus integration by parts: ( ) 2 u v + 2 u dx
More informationA collocation method for solving some integral equations in distributions
A collocation method for solving some integral equations in distributions Sapto W. Indratno Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA sapto@math.ksu.edu A G Ramm
More informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationThe Method of Fundamental Solutions applied to the numerical calculation of eigenfrequencies and eigenmodes for 3D simply connected domains
ECCOMAS Thematic Conference on Meshless Methods 2005 C34.1 The Method of Fundamental Solutions applied to the numerical calculation of eigenfrequencies and eigenmodes for 3D simply connected domains Carlos
More informationConnection to Laplacian in spherical coordinates (Chapter 13)
Connection to Laplacian in spherical coordinates (Chapter 13) We might often encounter the Laplace equation and spherical coordinates might be the most convenient 2 u(r, θ, φ) = 0 We already saw in Chapter
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationINVARIANT SUBSPACES FOR CERTAIN FINITE-RANK PERTURBATIONS OF DIAGONAL OPERATORS. Quanlei Fang and Jingbo Xia
INVARIANT SUBSPACES FOR CERTAIN FINITE-RANK PERTURBATIONS OF DIAGONAL OPERATORS Quanlei Fang and Jingbo Xia Abstract. Suppose that {e k } is an orthonormal basis for a separable, infinite-dimensional Hilbert
More informationKdV equation obtained by Lie groups and Sturm-Liouville problems
KdV equation obtained by Lie groups and Sturm-Liouville problems M. Bektas Abstract In this study, we solve the Sturm-Liouville differential equation, which is obtained by using solutions of the KdV equation.
More informationIBVPs for linear and integrable nonlinear evolution PDEs
IBVPs for linear and integrable nonlinear evolution PDEs Dionyssis Mantzavinos Department of Applied Mathematics and Theoretical Physics, University of Cambridge. Edinburgh, May 31, 212. Dionyssis Mantzavinos
More informationHigher-Order Compact Finite Element Method
Higher-Order Compact Finite Element Method Major Qualifying Project Advisor: Professor Marcus Sarkis Amorn Chokchaisiripakdee Worcester Polytechnic Institute Abstract The Finite Element Method (FEM) is
More informationNumerische Mathematik
Numer. Math. (997) 76: 479 488 Numerische Mathematik c Springer-Verlag 997 Electronic Edition Exponential decay of C cubic splines vanishing at two symmetric points in each knot interval Sang Dong Kim,,
More informationNumerical Solution of Two-Dimensional Volterra Integral Equations by Spectral Galerkin Method
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 159-174 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 Numerical Solution of Two-Dimensional Volterra
More informationf (t) K(t, u) d t. f (t) K 1 (t, u) d u. Integral Transform Inverse Fourier Transform
Integral Transforms Massoud Malek An integral transform maps an equation from its original domain into another domain, where it might be manipulated and solved much more easily than in the original domain.
More informationNumerical solution of the 2-D Poisson equation on an irregular domain with Robin boundary conditions
Numerical solution of the 2-D Poisson equation on an irregular domain with Robin boundary conditions Z. Jomaa C. Macaskill August 8, 28 Abstract We describe a 2-D finite difference algorithm for inverting
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationAnalytical solution for determination the control parameter in the inverse parabolic equation using HAM
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 2 (December 2017, pp. 1072 1087 Applications and Applied Mathematics: An International Journal (AAM Analytical solution
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationSloshing problem in a half-plane covered by a dock with two equal gaps
Sloshing prolem in a half-plane covered y a dock with two equal gaps O. V. Motygin N. G. Kuznetsov Institute of Prolems in Mech Engineering Russian Academy of Sciences St.Petersurg, Russia STATEMENT OF
More informationAlgebraic Spectral Relations forelliptic Quantum Calogero-MoserProblems
Journal of Nonlinear Mathematical Physics 1999, V.6, N 3, 263{268. Letter Algebraic Spectral Relations forelliptic Quantum Calogero-MoserProblems L.A. KHODARINOVA yz and I.A. PRIKHODSKY z y Wessex Institute
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationOn an uniqueness theorem for characteristic functions
ISSN 392-53 Nonlinear Analysis: Modelling and Control, 207, Vol. 22, No. 3, 42 420 https://doi.org/0.5388/na.207.3.9 On an uniqueness theorem for characteristic functions Saulius Norvidas Institute of
More informationequations in the semiclassical regime
Adaptive time splitting for nonlinear Schrödinger equations in the semiclassical regime Local error representation and a posteriori error estimator Workshop, Numerical Analysis of Evolution Equations Vill,
More informationQuasistatic Nonlinear Viscoelasticity and Gradient Flows
Quasistatic Nonlinear Viscoelasticity and Gradient Flows Yasemin Şengül University of Coimbra PIRE - OxMOS Workshop on Pattern Formation and Multiscale Phenomena in Materials University of Oxford 26-28
More informationAND NONLINEAR SCIENCE SERIES. Partial Differential. Equations with MATLAB. Matthew P. Coleman. CRC Press J Taylor & Francis Croup
CHAPMAN & HALL/CRC APPLIED MATHEMATICS AND NONLINEAR SCIENCE SERIES An Introduction to Partial Differential Equations with MATLAB Second Edition Matthew P Coleman Fairfield University Connecticut, USA»C)
More informationarxiv:quant-ph/ v1 11 Nov 1995
Quantization of electromagnetic fields in a circular cylindrical cavity arxiv:quant-ph/950v Nov 995 K. Kakazu Department of Physics, University of the Ryukyus, Okinawa 903-0, Japan Y. S. Kim Department
More informationHarmonic Analysis for Star Graphs and the Spherical Coordinate Trapezoidal Rule
Harmonic Analysis for Star Graphs and the Spherical Coordinate Trapezoidal Rule Robert Carlson University of Colorado at Colorado Springs September 24, 29 Abstract Aspects of continuous and discrete harmonic
More informationFast Structured Spectral Methods
Spectral methods HSS structures Fast algorithms Conclusion Fast Structured Spectral Methods Yingwei Wang Department of Mathematics, Purdue University Joint work with Prof Jie Shen and Prof Jianlin Xia
More informationNON-EXTINCTION OF SOLUTIONS TO A FAST DIFFUSION SYSTEM WITH NONLOCAL SOURCES
Electronic Journal of Differential Equations, Vol. 2016 (2016, No. 45, pp. 1 5. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu NON-EXTINCTION OF
More informationFourier series on triangles and tetrahedra
Fourier series on triangles and tetrahedra Daan Huybrechs K.U.Leuven Cambridge, Sep 13, 2010 Outline Introduction The unit interval Triangles, tetrahedra and beyond Outline Introduction The topic of this
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationNumerical techniques for linear and nonlinear eigenvalue problems in the theory of elasticity
ANZIAM J. 46 (E) ppc46 C438, 005 C46 Numerical techniques for linear and nonlinear eigenvalue problems in the theory of elasticity Aliki D. Muradova (Received 9 November 004, revised April 005) Abstract
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 informationExact fundamental solutions
Journées Équations aux dérivées partielles Saint-Jean-de-Monts, -5 juin 998 GDR 5 (CNRS) Exact fundamental solutions Richard Beals Abstract Exact fundamental solutions are known for operators of various
More information8 A pseudo-spectral solution to the Stokes Problem
8 A pseudo-spectral solution to the Stokes Problem 8.1 The Method 8.1.1 Generalities We are interested in setting up a pseudo-spectral method for the following Stokes Problem u σu p = f in Ω u = 0 in Ω,
More information100 CHAPTER 4. SYSTEMS AND ADAPTIVE STEP SIZE METHODS APPENDIX
100 CHAPTER 4. SYSTEMS AND ADAPTIVE STEP SIZE METHODS APPENDIX.1 Norms If we have an approximate solution at a given point and we want to calculate the absolute error, then we simply take the magnitude
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationFractional order operators on bounded domains
on bounded domains Gerd Grubb Copenhagen University Geometry seminar, Stanford University September 23, 2015 1. Fractional-order pseudodifferential operators A prominent example of a fractional-order pseudodifferential
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More information(IV.C) UNITS IN QUADRATIC NUMBER RINGS
(IV.C) UNITS IN QUADRATIC NUMBER RINGS Let d Z be non-square, K = Q( d). (That is, d = e f with f squarefree, and K = Q( f ).) For α = a + b d K, N K (α) = α α = a b d Q. If d, take S := Z[ [ ] d] or Z
More informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
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 informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationSPECTRAL METHODS: ORTHOGONAL POLYNOMIALS
SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS 31 October, 2007 1 INTRODUCTION 2 ORTHOGONAL POLYNOMIALS Properties of Orthogonal Polynomials 3 GAUSS INTEGRATION Gauss- Radau Integration Gauss -Lobatto Integration
More informationLinear and Multilinear Algebra. Linear maps preserving rank of tensor products of matrices
Linear maps preserving rank of tensor products of matrices Journal: Manuscript ID: GLMA-0-0 Manuscript Type: Original Article Date Submitted by the Author: -Aug-0 Complete List of Authors: Zheng, Baodong;
More informationPROBLEMS IN UNBOUNDED CYLINDRICAL DOMAINS
PROBLEMS IN UNBOUNDED CYLINDRICAL DOMAINS PATRICK GUIDOTTI Mathematics Department, University of California, Patrick Guidotti, 103 Multipurpose Science and Technology Bldg, Irvine, CA 92697, USA 1. Introduction
More informationQualifying Examination
Summer 24 Day. Monday, September 5, 24 You have three hours to complete this exam. Work all problems. Start each problem on a All problems are 2 points. Please email any electronic files in support of
More informationCompatible Hamiltonian Operators for the Krichever-Novikov Equation
arxiv:705.04834v [math.ap] 3 May 207 Compatible Hamiltonian Operators for the Krichever-Novikov Equation Sylvain Carpentier* Abstract It has been proved by Sokolov that Krichever-Novikov equation s hierarchy
More informationChebyshev polynomials for solving two dimensional linear and nonlinear integral equations of the second kind
Volume 31, N. 1, pp. 127 142, 2012 Copyright 2012 SBMAC ISSN 0101-8205 / ISSN 1807-0302 (Online) www.scielo.br/cam Chebyshev polynomials for solving two dimensional linear and nonlinear integral equations
More informationON THE EXISTENCE OF THREE SOLUTIONS FOR QUASILINEAR ELLIPTIC PROBLEM. Paweł Goncerz
Opuscula Mathematica Vol. 32 No. 3 2012 http://dx.doi.org/10.7494/opmath.2012.32.3.473 ON THE EXISTENCE OF THREE SOLUTIONS FOR QUASILINEAR ELLIPTIC PROBLEM Paweł Goncerz Abstract. We consider a quasilinear
More informationCONVEXITY OF INTEGRAL MEANS OF SUBHARMONIC FUNCTIONS
PROCEEDINGS OF THE AMERICAIN MATHEMATICAL SOCIETY Volume 60, October 1976 CONVEXITY OF INTEGRAL MEANS OF SUBHARMONIC FUNCTIONS JANG-MEI G. WU' Abstract. We study the convexity of integral means of subharmonic
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationStarting from Heat Equation
Department of Applied Mathematics National Chiao Tung University Hsin-Chu 30010, TAIWAN 20th August 2009 Analytical Theory of Heat The differential equations of the propagation of heat express the most
More informationClassification of partial differential equations and their solution characteristics
9 TH INDO GERMAN WINTER ACADEMY 2010 Classification of partial differential equations and their solution characteristics By Ankita Bhutani IIT Roorkee Tutors: Prof. V. Buwa Prof. S. V. R. Rao Prof. U.
More informationThe Helmholtz Equation
The Helmholtz Equation Seminar BEM on Wave Scattering Rene Rühr ETH Zürich October 28, 2010 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness
More informationBessel Functions Michael Taylor. Lecture Notes for Math 524
Bessel Functions Michael Taylor Lecture Notes for Math 54 Contents 1. Introduction. Conversion to first order systems 3. The Bessel functions J ν 4. The Bessel functions Y ν 5. Relations between J ν and
More informationDiffusion on the half-line. The Dirichlet problem
Diffusion on the half-line The Dirichlet problem Consider the initial boundary value problem (IBVP) on the half line (, ): v t kv xx = v(x, ) = φ(x) v(, t) =. The solution will be obtained by the reflection
More informationHere we used the multiindex notation:
Mathematics Department Stanford University Math 51H Distributions Distributions first arose in solving partial differential equations by duality arguments; a later related benefit was that one could always
More informationDunkl operators and Clifford algebras II
Translation operator for the Clifford Research Group Department of Mathematical Analysis Ghent University Hong Kong, March, 2011 Translation operator for the Hermite polynomials Translation operator for
More informationAlgebraic Properties of Dual Toeplitz Operators on Harmonic Hardy Space over Polydisc
Journal of Mathematical Research with Applications Nov., 2016, Vol. 36, No. 6, pp. 703 710 DOI:10.3770/j.issn:2095-2651.2016.06.009 Http://jmre.dlut.edu.cn Algebraic Properties of Dual Toeplitz Operators
More informationOn the convergence rate of a difference solution of the Poisson equation with fully nonlocal constraints
Nonlinear Analysis: Modelling and Control, 04, Vol. 9, No. 3, 367 38 367 http://dx.doi.org/0.5388/na.04.3.4 On the convergence rate of a difference solution of the Poisson equation with fully nonlocal
More informationNumerical Implementation of Fourier Integral-Transform Method for the Wave Equations
Numerical Implementation of Fourier Integral-Transform Method for the Wave Equations Michail D. Todorov Faculty of Applied Mathematics and Computer Science Technical University of Sofia, Bulgaria (Work
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georgia Tech PHYS 612 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #5 due October 2, 2012 == show all your work for maximum credit, == put labels, title,
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationSome aspects of the multivariable Mahler measure
Some aspects of the multivariable Mahler measure Séminaire de théorie des nombres de Chevaleret, Institut de mathématiques de Jussieu, Paris, France June 19th, 2006 Matilde N. Lalín Institut des Hautes
More informationDIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION
Meshless Methods in Science and Engineering - An International Conference Porto, 22 DIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION Robert Schaback Institut für Numerische und Angewandte Mathematik (NAM)
More informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More informationAn Efficient Numerical Method for Solving. the Fractional Diffusion Equation
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 1-12 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 An Efficient Numerical Method for Solving the Fractional
More informationLinearization of Two Dimensional Complex-Linearizable Systems of Second Order Ordinary Differential Equations
Applied Mathematical Sciences, Vol. 9, 2015, no. 58, 2889-2900 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.4121002 Linearization of Two Dimensional Complex-Linearizable Systems of
More informationProperties of the Scattering Transform on the Real Line
Journal of Mathematical Analysis and Applications 58, 3 43 (001 doi:10.1006/jmaa.000.7375, available online at http://www.idealibrary.com on Properties of the Scattering Transform on the Real Line Michael
More informationE.G. KALNINS AND WILLARD MILLER, JR. The notation used for -series and -integrals in this paper follows that of Gasper and Rahman [3].. A generalizati
A NOTE ON TENSOR PRODUCTS OF -ALGEBRA REPRESENTATIONS AND ORTHOGONAL POLYNOMIALS E.G. KALNINSy AND WILLARD MILLER, Jr.z Abstract. We work out examples of tensor products of distinct generalized s`) algebras
More informationThe Helically Reduced Wave Equation as a Symmetric Positive System
Utah State University DigitalCommons@USU All Physics Faculty Publications Physics 2003 The Helically Reduced Wave Equation as a Symmetric Positive System Charles G. Torre Utah State University Follow this
More informationSpectral Methods and Inverse Problems
Spectral Methods and Inverse Problems Omid Khanmohamadi Department of Mathematics Florida State University Outline Outline 1 Fourier Spectral Methods Fourier Transforms Trigonometric Polynomial Interpolants
More informationThe Journal of Nonlinear Sciences and Applications
J. Nonlinear Sci. Appl. 2 (2009), no. 3, 195 203 The Journal of Nonlinear Sciences Applications http://www.tjnsa.com ON MULTIPOINT ITERATIVE PROCESSES OF EFFICIENCY INDEX HIGHER THAN NEWTON S METHOD IOANNIS
More information