Linear Differential Equations as a Data-Structure
|
|
- Preston Newman
- 6 years ago
- Views:
Transcription
1 Linear Differential Equations as a Data-Structure Bruno Salvy Inria & ENS de Lyon FoCM, July 14, 2017
2 Computer Algebra Effective mathematics: what can we compute exactly? And complexity: how fast? (also, how big is the result?) Systems with several million users 50+ years of algorithmic progress in computational mathematics! 1/26
3 Sources of Linear Differential Equations Generating functions in combinatorics M. Bousquet-Mélou Periods Classical elementary and special functions (small order) A. Bostan P. Lairez 2/26
4 LDEs as a Data-Structure Polynomial equations Numerical evaluation Diagonals Definite sums and integrals Linear Differential Equations Conversions Local and asymptotic expansions Proofs of identities Closed forms M. Singer Solutions called differentially finite (abbrev. D-finite) 3/26
5 A. Using Linear Differential Equations Exactly
6 A. Using Linear Differential Equations Exactly I. Numerical Values
7 Fast Computation with Linear Recurrences (70 s and 80 s) 1. Multiplication of integers is fast (Fast Fourier Transform): 2. n! millions of digits 1sec. 3. Linear recurrence: convert into 1st order recurrence on vectors and apply the same idea. Ex: in complexity O(n) by divide-and-conquer n! :=n dn/2e {z } size O(n log n) e n := nx k=0 1 k! satisfies a 2nd order rec, computed via en e n 1 = 1 n (dn/2e 1) 1 {z } size O(n log n) n n {z 0 } A(n) en 1 e n 2 Conclusion: Nth element in O(N) ops. = 1 n! A!(n) 1 0 Notation: O(n) means O(n log k n) for some k. 4/26
8 Numerical evaluation of solutions of LDEs Principle: f(x) = NX n=0 a n x n {z } fast evaluation + 1X n=n+1 a n x n {z } good bounds f solution of a LDE with coeffs in Q(x) 1. linear recurrence in N for the first sum (easy); 2. tight bounds on the tail (technical); 3. extend to C by analytic continuation. Computation on integers. No roundoff errors. arctan(1+i) Conclusion: value anywhere with N digits in Õ(N) ops. Sage code available M. Mezzarobba [Chudnovsky-Chudnovsky87;van der Hoeven99;Mezzarobba-S.10;Mezzarobba16] 5/26
9 A. Using Linear Differential Equations Exactly II. Local and Asymptotic Expansions
10 Dynamic Dictionary of Mathematical Functions User need Recent algorithmic progress Maths on the web [Benoit-Chyzak-Darrasse-Gerhold-Mezzarobba-S.2010] 6/26
11 A. Using Linear Differential Equations Exactly III. Proofs of Identities
12 Proof technique > series(sin(x)^2+cos(x)^2-1,x,4); f satisfies a LDE f,f,f, live in a finite-dim. vector space O(x 4 ) Why is this a proof? 1. sin and cos satisfy a 2nd order LDE: y +y=0; 2. their squares and their sum satisfy a 3rd order LDE; 3. the constant -1 satisfies y =0; 4. thus sin 2 +cos 2-1 satisfies a LDE of order at most 4; 5. the Cauchy-Lipschitz theorem concludes. Proofs of non-linear identities by linear algebra! 7/26
13 Mehler s identity for Hermite polynomials 1X n=0 H n (x)h n (y) un n! = exp 4u(xy u(x 2 +y 2 )) 1 4u 2 p 1 4u 2 1. Definition of Hermite polynomials: recurrence of order 2; 2. Product by linear algebra: H n+k (x)h n+k (y)/(n+k)!, k N generated over H n (x)h n (y) n! Q (x,n) by, H n+1(x)h n (y) n!, H n(x)h n+1 (y) n! recurrence of order at most 4; 3. Translate into differential equation., H n+1(x)h n+1 (y) n! 8/26
14 Guess & Prove Continued Fractions 1. Taylor expansion produces first terms (easy): arctan x = x 1 3 x x x Guess a formula (easy): a n = n2 4n Prove that the CF with these a n converges to arctan: show that H n := Q 2 n (x 2 + 1)(P n /Q n ) 0 1 = O(x n ) where P n /Q n is the nth convergent. Algo compute a LRE for H n and simplify it. gfun[contfrac] No human intervention needed. 9/26
15 It Works! This method has been applied to all explicit C-fractions in Cuyt et alii, starting from either: a Riccati equation: a q-riccati equation: y 0 = A(z)+B(z)y + C(z)y 2 y(qz) =A(z)+B(z)y(z)+C(z)y(z)y(qz) a difference Riccati equation: y(s + 1) = A(s) + B(s)y(s) + C(s)y(s)y(s + 1) It works in all cases, including Gauss s CF, Heine s q- analogue and Brouncker s CF for Gamma. In all cases, H n satisfies a recurrence of small order. In progress: 1. explain why this method works so well, [Maulat-S. 15,17] 2. classify the formulas it yields. 10/26
16 B. Conversions (LDE LDE)
17 From equations to operators D x d/dx x mult by x product composition D x x=xd x +1 S n (n n+1) n mult by n product composition S n n=(n+1)s n Ex. (erf): Taylor morphism: D x (n+1)s n ; x S n -1 produces linear recurrence from LDE D 2 x +2xD x 7! (n + 1)S n (n + 1)S n +2S 1 n (n + 1)S n =(n + 1)(n + 2)S 2 n +2n 11/26
18 Chebyshev expansions arctan Taylor Chebyshev z 1 3 z z5 + 2( p T1 (x) 2 + 1) (2 p 2 + 3) c k = 2 Z 1 T 3 (x) 3(2 p 2 + 3) + T 5 (x) 2 5(2 p 2 + 3) f(x)t k (x) p 1 x 2 dx 12/26
19 Ore fractions Generalize commutative case: R=Q -1 P with P & Q operators. B -1 A=D -1 C when ba=dc with bb=dd=lclm(b,d). Algorithms for sum and product: B -1 A+D -1 C=LCLM(B,D) -1 (ba+dc), with bb=dd=lclm(b,d) B -1 AD -1 C=(aB) -1 dc, with aa=dd=lclm(a,d). [Ore1933] 13/26
20 Application: Chebyshev expansions Taylor x n+1 =x x n x X:=S -1 (x n ) =nx n-1 d/dx D:=(n+1)S Chebyshev 2xT n (x)=t n+1 (x)+t n-1 (x) x X:=(S n +S n -1 )/2 2(1-x 2 )T n (x)=-nt n+1 (x)+nt n-1 (x) d/dx D:=(1-X 2 ) -1 n(s n -S n -1 )/2. Prop. If y is a solution of L(x,d/dx), then its Chebyshev coefficients annihilate the numerator of L(X,D). > deqarctan:=(x^2+1)*diff(y(x),x)-1: > diffeqtogfsrec(deqarctan,y(x),u(n),functions=chebyshevt(n,x)); nu(n) + 6(n + 2)u(n + 2) + (n + 4)u(n + 4) Applications to Validated Numerical Approximation [Benoit-S.09;Benoit12;BenoitJoldesMezzarobba17] M. Joldes 14/26
21 C. Computing Linear Differential Equations (Efficiently)
22 C. Computing Linear Differential Equations (Efficiently) I. Algebraic Series and Questions of Size
23 Algebraic Series can be Computed Fast P (X, Y (X)) = 0 P irreducible Wanted: the first N Taylor coefficients of Y. P x (X, Y (X)) + P y (X, Y (X)) Y 0 (X) =0 Y 0 (X) =( P x P 1 y mod P )(X, Y (X)) a polynomial Note: F sol LDE F(Y(X)) sol LDE (same argument) Y (X),Y 0 (X),Y 00 (X),...in Vect Q(X) (1,Y,Y 2,...) finite dimension (deg P) a LDE by linear algebra [Abel1827;Cockle1861;Harley1862;Tannery1875] 15/26
24 degree O(D^3) Order-Degree Curve The cost of minimality minimal LDE nice recurrence minimal recurrence degree D=deg P O(D^2) O(D^2) O(D^2) O(D) O(D) O(D) O(D) O(D^2) order O(D^2) O(D^2) O(D^3) order differential equations corresponding recurrences [Bostan-Chyzak-Lecerf-S.-Schost07;Chen-Kauers12;Chen-Jaroschek-Kauers-Singer13] 16/26
25 C. Computing Linear Differential Equations (Efficiently) II. Creative Telescoping
26 Examples X j,k nx 2 2 n n + k nx kx 3 n n + k k = k k k k j k=0 k=0 j=0 j + k r n s + n j k n + r s r ( 1) j+k =( 1) l k + l j k m j n + l m n l Z i xj 1 (ax)i 1 (ax)y 0 (x)k 0 (x) dx = ln(1 a4 ) 2 a 2 I (1 + 2xy + 4y 2 )exp 4x 2 y 2 1+4y 2 dy = H n(x) y n+1 (1 + 4y 2 ) 3 2 bn/2c! Aims: nx k=0 q k2 (q; q) k (q; q) n k = nx k= n ( 1) k q (5k2 k)/2 (q; q) n 1. Prove them automatically 2. Find the rhs given the lhs k (q; q) n+k Note: at least one free variable First: find a LDE (or LRE) 17/26
27 I(x) = Z Creative telescoping f(x, t) dt =? or U(n) = X k u(n, k) =? Input: equations (differential for f or recurrence for u). Output: equations for the sum or the integral. Ex.: U n := X k U n+1 2U n = X k n k n +1 k n 2 k = X k n +1 n +1 k k +1 {z } telescopes + n n k +1 k {z } telescopes Aim: find A(n,S n ) and B(n,k,S n,s k ) such that (A(n,S n )+Δ k B(n,k,S n,s k )) u(n,k)=0, Def: k: =S k -1. certificate then the sum telescopes, leading to A(n,S n ) U(n)=0. Integrals: differentiate under the sign and integrate by parts. 18/26
28 Telescoping Ideal T t (f) := Ann f t Q(x,t)h@ x,@ t i {z } int. by parts (certificate) \ Q(x)h@ x i. {z } di. under R First generation of algorithms relying on holonomy Restrict int. by parts to Q(x)h@ x,@ t i and use elimination. Second generation: faster using better certificates & algorithms Hypergeometric summation: dim=1 + param. Gosper. Undetermined coefficients in finite dim, Ore algebras & GB. Idem in infinite dim. X X c k (x)@ x t a i,j (x, t)@ x@ i j t 2 Ann f k i,j2s t z z F. Chyzak [Zeilberger et alii 90,91,92;Chyzak00;Chyzak-Kauers-S.09] 28 x x x y y 19/26
29 C. Computing Linear Differential Equations (Efficiently) III. 3rd Generation Creative Telescoping
30 Certificates are big n n + s n + r 2n r s C n := X r,s n ( 1) n+r+s r s s r n {z } f n,r,s (n + 2) 3 C n+2 2(2n + 3)(3n 2 + 9n + 7)C n+1 (4n + 3)(4n + 4)(4n + 5)C n = 180 kb ' 2 pages I(z) = I (1 + t 3 ) 2 dt 1 dt 2 dt 3 t 1 t 2 t 3 (1 + t 3 (1 + t 1 ))(1 + t 3 (1 + t 2 )) + z(1 + t 1 )(1 + t 2 )(1 + t 3 ) 4 z 2 (4z + 1)(16z 1)I 000 (z)+3z(128z z 1)I 00 (z)+(444z z 1)I 0 (z)+2(30z + 1)I(z) =1 080 kb ' 12 pages 3rd-generation algorithms: avoid computing the certificate 20/26
31 Periods I(t) = N := degx Q, I P(t, x) dx m Q (t, x) {z } 2Q(t,x) Q square-free Int. over a cycle where Q 0. x = (x1,..., xn ) dt := max(degt Q, degt P) degxp not too big Thm. A linear differential equation for I(t) can be computed in O(e3nN8ndt) operations in ℚ. It has order Nn and degree O(enN3ndt). tight Note: generically, the certificate has at least N n2/2 monomials. Applications to diagonals & to multiple binomial sums. [Picard1902;Dwork62;Griffiths69;Christol85;Bostan-Lairez-Salvy13;Lairez16] 21/26
32 Diagonals in this talk If F (z) = G(z) H(z) its diagonal is 1 k +1 2k k 2k k Apéry s a k : : : is a multivariate rational function with Taylor expansion F (z) = X i2n n c i z i, F (t) = X k2n c k,k,...,k t k. 1 1 x y =1+x + y +2xy + x2 + y x 2 y x (1 x y)(1 x) =1+y+1xy x2 +y x 2 y t(1 + x)(1 + y)(1 + z)(1 + y + z + yz + xyz) =1+ +5xyzt + Christol s conjecture: All differentially finite power series with integer coefficients and radius of convergence in (0, ) are diagonals. 22/26
33 Diagonals are Differentially Finite [Christol84,Lipshitz88] F (z 1,...,z d )= d 1 2 i 1 I F t dz2,z 2,...,z d z 2 z d z 2 dz d z d Thm. If F has degree d in n variables, ΔF satisfies a LDE with order d n, coeffs of degree d O(n). D-finite diag. alg. rat. Univariate power series + algo in ops. Õ(d 8n ) [Bostan-Lairez-S.13;Lairez16] 23/26
34 Multiple Binomial Sums Ex. S n = X r 0 X n ( 1) n+r+s r s 0 n s n + s n + r 2n r s s r n Thm. Diagonals binomial sums with 1 free index. > BinomSums[sumtores](S,u): ( ) defined properly 1 1 t(1 + u 1 )(1 + u 2 )(1 u 1 u 3 )(1 u 2 u 3 ) has for diagonal the generating function of S n LDE LRE [Bostan-Lairez-S.17] 24/26
35 (Non-)Commercial Algorithmes E caces en Calcul Formel z Alin Bostan Frédéric Chyzak Marc Giusti Romain Lebreton Grégoire Lecerf Bruno Salvy Éric Schost New book ( 700p.), based on our course. Freely available from our web pages, forever. Paper version before the end of y x 25/26
36 Conclusion Combinatorics Approximation Linear Differential Equations + Computer Algebra Special Functions The End Your Work? 26/26
Linear differential and recurrence equations viewed as data-structures
Linear differential and recurrence equations viewed as data-structures Bruno Salvy Inria & ENS de Lyon Leçons de Mathématiques et d'informatique d'aujourd'hui Bordeaux, Mai 2015 Computer Algebra Effective
More informationAlgorithmic Tools for the Asymptotics of Diagonals
Algorithmic Tools for the Asymptotics of Diagonals Bruno Salvy Inria & ENS de Lyon Lattice walks at the Interface of Algebra, Analysis and Combinatorics September 19, 2017 Asymptotics & Univariate Generating
More informationDiagonals: Combinatorics, Asymptotics and Computer Algebra
Diagonals: Combinatorics, Asymptotics and Computer Algebra Bruno Salvy Inria & ENS de Lyon Families of Generating Functions (a n ) 7! A(z) := X n 0 a n z n counts the number of objects of size n captures
More informationAlgorithmic Tools for the Asymptotics of Linear Recurrences
Algorithmic Tools for the Asymptotics of Linear Recurrences Bruno Salvy Inria & ENS de Lyon aaacuxicbvhlbtraejw1r2beg+aghxhrsgitvjycahkhcc7cskqsibs2vunxez3jpjyznsli+rf4bq7wbfwknz6f8torzenli5wqutvtxvklhcmo+juirl2/cfpwxu3wzt179x8mnx9+cqa2hkbcsgmpm+zacg1tfcjhsllavcbhidt+3+khn8e6yfq+litifvtouqjo0fpz4epelmacmfmvbkzgpnfwes2hw9ekwhw9cujzslxz9ofeb0li7nxz8dxjda8vaossotelowrthlkuxeibjhih/m96gwwtijzhqcixbjtxevdtmnqoksap2qjmhmqmwkxn2cpds0eeymlhrh8a6yr9d6rhyrmlynynyli6da0j/6fnaixep43qvy2geb+oqcvfq7tt0vxy4cixhjbuhbdcecks4+gpgiywnjxyoxtt+yveftkurjby1kzq19wcksgxfx1/m/3sef2o/q5rucybdln2pk1pxjyvctd9lapqmhbugrxfjtf7cdp0g52v0ccyr+d1fuxftt5m4j2f5dvs1wz5qp6r5yqm22shfcc7zeo4+uk+ke/kr/a2yiiyoopbg8h5zcnyqql3bxdu3c8=
More informationAlgorithmic Tools for the Asymptotics of Linear Recurrences
Algorithmic Tools for the Asymptotics of Linear Recurrences Bruno Salvy Inria & ENS de Lyon Computer Algebra in Combinatorics, Schrödinger Institute, Vienna, Nov. 2017 Motivation p 0 (n)a n+k + + p k (n)a
More informationThe Dynamic Dictionary of Mathematical Functions
The Dynamic Dictionary of Mathematical Functions Bruno.Salvy@inria.fr AISC 2010. July 6, 2010 Joint work with: Alexandre Benoit, Alin Bostan, Frédéric Chyzak, Alexis Darrasse, Stefan Gerhold, Marc Mezzarobba
More informationAutomatic Proofs of Identities: Beyond A=B
Bruno.Salvy@inria.fr FPSAC, Linz, July 20, 2009 Joint work with F. Chyzak and M. Kauers 1 / 28 I Introduction 2 / 28 Examples of Identities: Definite Sums, q-sums, Integrals n k=0 + 0 1 2πi n=0 ( ) H n
More informationDiagonals and Walks. Alin Bostan 1 Louis Dumont 1 Bruno Salvy 2. November 4, Introduction Diagonals Walks. 1 INRIA, SpecFun.
and Alin Bostan 1 Louis Dumont 1 Bruno Salvy 2 1 INRIA, SpecFun 2 INRIA, AriC November 4, 2014 1/16 Louis Dumont and 2/16 Louis Dumont and Motivations In combinatorics: Lattice walks w i,j : number of
More informationBon anniversaire Marc!
Bon anniversaire Marc! The full counting function of Gessel walks is algebraic Alin Bostan (INRIA) joint work with Manuel Kauers (INRIA/RISC) Journées GECKO/TERA 2008 Context: restricted lattice walks
More informationSymbolic-Numeric Tools for Multivariate Asymptotics
Symbolic-Numeric Tools for Multivariate Asymptotics Bruno Salvy Joint work with Stephen Melczer to appear in Proc. ISSAC 2016 FastRelax meeting May 25, 2016 Combinatorics, Randomness and Analysis From
More informationEfficient Experimental Mathematics for Lattice Path Combinatorics
1/45 Efficient Experimental Mathematics for Lattice Path Combinatorics Alin Bostan based on joint works with M. Bousquet-Mélou, F. Chyzak, M. van Hoeij, M. Kauers, S. Melczer, L. Pech, K. Raschel, B. Salvy
More informationNeed for transformation?
Z-TRANSFORM In today s class Z-transform Unilateral Z-transform Bilateral Z-transform Region of Convergence Inverse Z-transform Power Series method Partial Fraction method Solution of difference equations
More informationFast Computation of the Nth Term of an Algebraic Series over a Finite Prime Field
Fast Computation of the Nth Term of an Algebraic Series over a Finite Prime Field Alin Bostan 1 Gilles Christol 2 Philippe Dumas 1 1 Inria (France) 2 Institut de mathématiques de Jussieu (France) Specfun
More informationSection Taylor and Maclaurin Series
Section.0 Taylor and Maclaurin Series Ruipeng Shen Feb 5 Taylor and Maclaurin Series Main Goal: How to find a power series representation for a smooth function us assume that a smooth function has a power
More informationChapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27
Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation
More informationLucas Congruences. A(n) = Pure Mathematics Seminar University of South Alabama. Armin Straub Oct 16, 2015 University of South Alabama.
Pure Mathematics Seminar University of South Alabama Oct 6, 205 University of South Alabama A(n) = n k=0 ( ) n 2 ( n + k k k ) 2, 5, 73, 445, 3300, 89005, 2460825,... Arian Daneshvar Amita Malik Zhefan
More informationClosed Form Solutions
Closed Form Solutions Mark van Hoeij 1 Florida State University Slides of this talk: www.math.fsu.edu/ hoeij/issac2017.pdf 1 Supported by NSF 1618657 1 / 29 What is a closed form solution? Example: Solve
More informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More informationLecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain.
Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain. For example f(x) = 1 1 x = 1 + x + x2 + x 3 + = ln(1 + x) = x x2 2
More information25. Chain Rule. Now, f is a function of t only. Expand by multiplication:
25. Chain Rule The Chain Rule is present in all differentiation. If z = f(x, y) represents a two-variable function, then it is plausible to consider the cases when x and y may be functions of other variable(s).
More information1 Question related to polynomials
07-08 MATH00J Lecture 6: Taylor Series Charles Li Warning: Skip the material involving the estimation of error term Reference: APEX Calculus This lecture introduced Taylor Polynomial and Taylor Series
More informationMath Assignment 11
Math 2280 - Assignment 11 Dylan Zwick Fall 2013 Section 8.1-2, 8, 13, 21, 25 Section 8.2-1, 7, 14, 17, 32 Section 8.3-1, 8, 15, 18, 24 1 Section 8.1 - Introduction and Review of Power Series 8.1.2 - Find
More informationMathematics 324 Riemann Zeta Function August 5, 2005
Mathematics 324 Riemann Zeta Function August 5, 25 In this note we give an introduction to the Riemann zeta function, which connects the ideas of real analysis with the arithmetic of the integers. Define
More informationSome Open Problems Related to Creative Telescoping
J Syst Sci Complex (2017) 30: 154 172 Some Open Problems Related to Creative Telescoping CHEN Shaoshi KAUERS Manuel DOI: 10.1007/s11424-017-6202-9 Received: 14 September 2016 / Revised: 30 October 2016
More informationClosed Form Solutions
Closed Form Solutions Mark van Hoeij 1 Florida State University ISSAC 2017 1 Supported by NSF 1618657 1 / 26 What is a closed form solution? Example: Solve this equation for y = y(x). y = 4 x 3 (1 x) 2
More informationand modular forms Journal of Physics A: Mathematical and Theoretical PAPER
Journal of Physics A: Mathematical and Theoretical PAPER Diagonals of rational functions, pullbacked and modular forms hypergeometric functions To cite this article: 208 J. Phys. A: Math. Theor. 5 45520
More informationarxiv: v1 [cs.sc] 15 Oct 2015
ALGEBRAIC DIAGONALS AND WALKS: ALGORITHMS, BOUNDS, COMPLEXITY ALIN BOSTAN, LOUIS DUMONT, AND BRUNO SALVY arxiv:1510.04526v1 [cs.sc] 15 Oct 2015 Abstract. The diagonal of a multivariate power series F is
More informationSeries Solutions of Linear ODEs
Chapter 2 Series Solutions of Linear ODEs This Chapter is concerned with solutions of linear Ordinary Differential Equations (ODE). We will start by reviewing some basic concepts and solution methods for
More informationFree Meixner distributions and random matrices
Free Meixner distributions and random matrices Michael Anshelevich July 13, 2006 Some common distributions first... 1 Gaussian Negative binomial Gamma Pascal chi-square geometric exponential 1 2πt e x2
More informationPhilippe Flajolet and Symbolic Computation
Philippe Flajolet and Symbolic Computation Bruno Salvy 1. Combinatorics and Analysis of Algorithms 4 2. Complex Analysis 5 3. Symbolic Computation 7 While Philippe Flajolet did not publish many articles
More informationDiagonals of rational functions, pullbacked 2 F 1 hypergeometric functions and modular forms (unabrigded version)
Diagonals of rational functions, pullbacked 2 F hypergeometric functions and modular forms unabrigded version Y. Abdelaziz S. Boukraa, C. Koutschan, J-M. Maillard LPTMC, UMR 7600 CNRS, Université Pierre
More informationMath 118, Handout 4: Hermite functions and the Fourier transform. n! where D = d/dx, are a basis of eigenfunctions for the Fourier transform
The Hermite functions defined by h n (x) = ( )n e x2 /2 D n e x2 where D = d/dx, are a basis of eigenfunctions for the Fourier transform f(k) = f(x)e ikx dx on L 2 (R). Since h 0 (x) = e x2 /2 h (x) =
More informationIntroduction and Review of Power Series
Introduction and Review of Power Series Definition: A power series in powers of x a is an infinite series of the form c n (x a) n = c 0 + c 1 (x a) + c 2 (x a) 2 +...+c n (x a) n +... If a = 0, this is
More informationarxiv: v1 [cs.sc] 7 Jul 2016
Rigorous Multiple-Precision Evaluation of D-Finite Functions in SageMath Marc Mezzarobba CNRS, LIP6, Université Pierre et Marie Curie, Paris, France marc@mezzarobba.net, http://marc.mezzarobba.net/ arxiv:1607.01967v1
More informationIs Analysis Necessary?
Is Analysis Necessary? Ira M. Gessel Brandeis University Waltham, MA gessel@brandeis.edu Special Session on Algebraic and Analytic Combinatorics AMS Fall Eastern Meeting University of Connecticut, Storrs
More informationMultiple binomial sums
Multiple binomial sums Alin Bostan, Pierre Lairez, Bruno Salvy To cite this version: Alin Bostan, Pierre Lairez, Bruno Salvy. Multiple binomial sums. Journal of Symbolic Computation, Elsevier, 207, 80
More information5.4 Bessel s Equation. Bessel Functions
SEC 54 Bessel s Equation Bessel Functions J (x) 87 # with y dy>dt, etc, constant A, B, C, D, K, and t 5 HYPERGEOMETRIC ODE At B (t t )(t t ), t t, can be reduced to the hypergeometric equation with independent
More informationOn the singularities of non-linear ODEs
On the singularities of non-linear ODEs Galina Filipuk Institute of Mathematics University of Warsaw G.Filipuk@mimuw.edu.pl Collaborators: R. Halburd (London), R. Vidunas (Tokyo), R. Kycia (Kraków) 1 Plan
More informationDesingularization of Ore Operators
Desingularization of Ore Operators Shaoshi Chen 1 KLMM, AMSS Chinese Academy of Sciences 100190 Beijing, China Manuel Kauers 2 Institute for Algebra Johannes Kepler University Linz, Austria Michael F Singer
More informationq-binomial coefficient congruences
CARMA Analysis and Number Theory Seminar University of Newcastle August 23, 2011 Tulane University, New Orleans 1 / 35 Our first -analogs The natural number n has the -analog: [n] = n 1 1 = 1 + +... n
More informationExperimental mathematics and integration
Experimental mathematics and integration David H. Bailey http://www.davidhbailey.com Lawrence Berkeley National Laboratory (retired) Computer Science Department, University of California, Davis October
More informationMATHEMATICAL FORMULAS AND INTEGRALS
MATHEMATICAL FORMULAS AND INTEGRALS ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom Academic Press San Diego New York Boston London
More informationClassnotes - MA Series and Matrices
Classnotes - MA-2 Series and Matrices Department of Mathematics Indian Institute of Technology Madras This classnote is only meant for academic use. It is not to be used for commercial purposes. For suggestions
More informationFourth Week: Lectures 10-12
Fourth Week: Lectures 10-12 Lecture 10 The fact that a power series p of positive radius of convergence defines a function inside its disc of convergence via substitution is something that we cannot ignore
More informationRigorous Polynomial Approximations and Applications
Rigorous Polynomial Approximations and Applications Mioara Joldeș under the supervision of: Nicolas Brisebarre and Jean-Michel Muller École Normale Supérieure de Lyon, Arénaire Team, Laboratoire de l Informatique
More informationTrading Order for Degree in Creative Telescoping
Trading Order for Degree in Creative Telescoping Shaoshi Chen 1 Department of Mathematics North Carolina State University Raleigh, NC 27695-8205, USA Manuel Kauers 2 Research Institute for Symbolic Computation
More informationMATHEMATICAL FORMULAS AND INTEGRALS
HANDBOOK OF MATHEMATICAL FORMULAS AND INTEGRALS Second Edition ALAN JEFFREY Department of Engineering Mathematics University of Newcastle upon Tyne Newcastle upon Tyne United Kingdom ACADEMIC PRESS A Harcourt
More informationDepartment of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016
Department of Mathematics, University of California, Berkeley YOUR 1 OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016 1. Please write your 1- or 2-digit exam number on
More information1. Find the Taylor series expansion about 0 of the following functions:
MAP 4305 Section 0642 3 Intermediate Differential Equations Assignment 1 Solutions 1. Find the Taylor series expansion about 0 of the following functions: (i) f(z) = ln 1 z 1+z (ii) g(z) = 1 cos z z 2
More informationSec 2.2: Infinite Limits / Vertical Asymptotes Sec 2.6: Limits At Infinity / Horizontal Asymptotes
Sec 2.2: Infinite Limits / Vertical Asymptotes Sec 2.6: Limits At Infinity / Horizontal Asymptotes Sec 2.2: Infinite Limits / Vertical Asymptotes Sec 2.6: Limits At Infinity / Horizontal Asymptotes Infinite
More informationSec 2.2: Infinite Limits / Vertical Asymptotes Sec 2.6: Limits At Infinity / Horizontal Asymptotes
Sec 2.2: Infinite Limits / Vertical Asymptotes Sec 2.6: Limits At Infinity / Horizontal Asymptotes Sec 2.2: Infinite Limits / Vertical Asymptotes Sec 2.6: Limits At Infinity / Horizontal Asymptotes Infinite
More informationJournal of Symbolic Computation
Journal of Symbolic Computation 74 (216) 617 626 Contents lists available at ScienceDirect Journal of Symbolic Computation wwwelseviercom/locate/jsc Desingularization of Ore operators Shaoshi Chen a,1,
More informationPolynomial evaluation and interpolation on special sets of points
Polynomial evaluation and interpolation on special sets of points Alin Bostan and Éric Schost Laboratoire STIX, École polytechnique, 91128 Palaiseau, France Abstract We give complexity estimates for the
More informationHigher-order ordinary differential equations
Higher-order ordinary differential equations 1 A linear ODE of general order n has the form a n (x) dn y dx n +a n 1(x) dn 1 y dx n 1 + +a 1(x) dy dx +a 0(x)y = f(x). If f(x) = 0 then the equation is called
More informationSOLUTIONS ABOUT ORDINARY POINTS
238 CHAPTER 6 SERIES SOLUTIONS OF LINEAR EQUATIONS In Problems 23 and 24 use a substitution to shift the summation index so that the general term of given power series involves x k. 23. nc n x n2 n 24.
More informationAutomatisches Beweisen von Identitäten
Automatisches Beweisen von Identitäten Stefan Gerhold 30. Oktober 2006 The General Scheme Parametric sums summand(n, k) = answer(n) Parametric integrals k integrand(x, y)dy = answer(x) Zeilberger s Summation
More informationEXAMPLES OF PROOFS BY INDUCTION
EXAMPLES OF PROOFS BY INDUCTION KEITH CONRAD 1. Introduction In this handout we illustrate proofs by induction from several areas of mathematics: linear algebra, polynomial algebra, and calculus. Becoming
More informationConstructing Taylor Series
Constructing Taylor Series 8-8-200 The Taylor series for fx at x = c is fc + f cx c + f c 2! x c 2 + f c x c 3 + = 3! f n c x c n. By convention, f 0 = f. When c = 0, the series is called a Maclaurin series.
More informationDifferential Equations for Algebraic Functions
Differential Equations for Algebraic Functions Alin Bostan, Frédéric Chyzak, Bruno Salvy Algorithms Project, Inria Rocquencourt 78153 Le Chesnay (France) {Alin.Bostan,Frederic.Chyzak,Bruno.Salvy}@inria.fr
More informationAP Calculus Testbank (Chapter 9) (Mr. Surowski)
AP Calculus Testbank (Chapter 9) (Mr. Surowski) Part I. Multiple-Choice Questions n 1 1. The series will converge, provided that n 1+p + n + 1 (A) p > 1 (B) p > 2 (C) p >.5 (D) p 0 2. The series
More informationRational and Algebraic General Solutions of First-Order Algebraic ODEs
Contents Rational and Algebraic General Solutions of First-Order Algebraic ODEs N. Thieu Vo RISC, Johannes Kepler University, Linz, Austria KSDA Seminar, The Graduate Center, CUNY, New York April 1, 2016
More informationThe (q, t)-catalan Numbers and the Space of Diagonal Harmonics. James Haglund. University of Pennsylvania
The (q, t)-catalan Numbers and the Space of Diagonal Harmonics James Haglund University of Pennsylvania Outline Intro to q-analogues inv and maj q-catalan Numbers MacMahon s q-analogue The Carlitz-Riordan
More informationCOMPUTER ALGEBRA AND POWER SERIES WITH POSITIVE COEFFICIENTS. 1. Introduction
COMPUTER ALGEBRA AND POWER SERIES WITH POSITIVE COEFFICIENTS MANUEL KAUERS Abstract. We consider the question whether all the coefficients in the series expansions of some specific rational functions are
More informationOn integral representations of q-gamma and q beta functions
On integral representations of -gamma and beta functions arxiv:math/3232v [math.qa] 4 Feb 23 Alberto De Sole, Victor G. Kac Department of Mathematics, MIT 77 Massachusetts Avenue, Cambridge, MA 239, USA
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationEvaluating a Determinant
Evaluating a Determinant Thotsaporn Thanatipanonda Joint Work with Christoph Koutschan Research Institute for Symbolic Computation, Johannes Kepler University of Linz, Austria June 15th, 2011 Outline of
More informationSeries Solutions. 8.1 Taylor Polynomials
8 Series Solutions 8.1 Taylor Polynomials Polynomial functions, as we have seen, are well behaved. They are continuous everywhere, and have continuous derivatives of all orders everywhere. It also turns
More informationDifferential Equations and Associators for Periods
Differential Equations and Associators for Periods Stephan Stieberger, MPP München Workshop on Geometry and Physics in memoriam of Ioannis Bakas November 2-25, 26 Schloß Ringberg, Tegernsee based on: St.St.,
More informationg(t) = f(x 1 (t),..., x n (t)).
Reading: [Simon] p. 313-333, 833-836. 0.1 The Chain Rule Partial derivatives describe how a function changes in directions parallel to the coordinate axes. Now we shall demonstrate how the partial derivatives
More informationReview For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: = 0 : homogeneous equation.
Review For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: y y x y2 = 0 : homogeneous equation. x2 v = y dy, y = vx, and x v + x dv dx = v + v2. dx =
More information3.4 Introduction to power series
3.4 Introduction to power series Definition 3.4.. A polynomial in the variable x is an expression of the form n a i x i = a 0 + a x + a 2 x 2 + + a n x n + a n x n i=0 or a n x n + a n x n + + a 2 x 2
More informationA Proof of a Recursion for Bessel Moments
A Proof of a Recursion for Bessel Moments Jonathan M. Borwein Faculty of Computer Science, Dalhousie University, Halifax, NS, B3H 2W5, Canada jborwein@cs.dal.ca Bruno Salvy Algorithms Project, Inria Paris-Rocquencourt,
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationEngg. Math. I. Unit-I. Differential Calculus
Dr. Satish Shukla 1 of 50 Engg. Math. I Unit-I Differential Calculus Syllabus: Limits of functions, continuous functions, uniform continuity, monotone and inverse functions. Differentiable functions, Rolle
More informationDifferential Equations for Algebraic Functions
Differential Equations for Algebraic Functions Alin Bostan, Frédéric Chyzak, Bruno Salvy, Grégoire Lecerf, Éric Schost To cite this version: Alin Bostan, Frédéric Chyzak, Bruno Salvy, Grégoire Lecerf,
More informationPUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES. Notes
PUTNAM PROBLEMS SEQUENCES, SERIES AND RECURRENCES Notes. x n+ = ax n has the general solution x n = x a n. 2. x n+ = x n + b has the general solution x n = x + (n )b. 3. x n+ = ax n + b (with a ) can be
More informationTEST CODE: MMA (Objective type) 2015 SYLLABUS
TEST CODE: MMA (Objective type) 2015 SYLLABUS Analytical Reasoning Algebra Arithmetic, geometric and harmonic progression. Continued fractions. Elementary combinatorics: Permutations and combinations,
More informationSolutions: Problem Set 3 Math 201B, Winter 2007
Solutions: Problem Set 3 Math 201B, Winter 2007 Problem 1. Prove that an infinite-dimensional Hilbert space is a separable metric space if and only if it has a countable orthonormal basis. Solution. If
More informationCompletion Date: Monday February 11, 2008
MATH 4 (R) Winter 8 Intermediate Calculus I Solutions to Problem Set #4 Completion Date: Monday February, 8 Department of Mathematical and Statistical Sciences University of Alberta Question. [Sec..9,
More informationElementary Differential Equations, Section 2 Prof. Loftin: Practice Test Problems for Test Find the radius of convergence of the power series
Elementary Differential Equations, Section 2 Prof. Loftin: Practice Test Problems for Test 2 SOLUTIONS 1. Find the radius of convergence of the power series Show your work. x + x2 2 + x3 3 + x4 4 + + xn
More informationMATH 452. SAMPLE 3 SOLUTIONS May 3, (10 pts) Let f(x + iy) = u(x, y) + iv(x, y) be an analytic function. Show that u(x, y) is harmonic.
MATH 45 SAMPLE 3 SOLUTIONS May 3, 06. (0 pts) Let f(x + iy) = u(x, y) + iv(x, y) be an analytic function. Show that u(x, y) is harmonic. Because f is holomorphic, u and v satisfy the Cauchy-Riemann equations:
More informationEngg. Math. II (Unit-IV) Numerical Analysis
Dr. Satish Shukla of 33 Engg. Math. II (Unit-IV) Numerical Analysis Syllabus. Interpolation and Curve Fitting: Introduction to Interpolation; Calculus of Finite Differences; Finite Difference and Divided
More informationSupercongruences for Apéry-like numbers
Supercongruences for Apéry-like numbers AKLS seminar on Automorphic Forms Universität zu Köln March, 205 University of Illinois at Urbana-Champaign A(n) = n k=0 ( ) n 2 ( n + k k k ) 2, 5, 73, 445, 3300,
More informationTaylor and Maclaurin Series
Taylor and Maclaurin Series MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Background We have seen that some power series converge. When they do, we can think of them as
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 informationSingularity analysis via the iterated kernel method
FPSAC 2013 Paris, France DMTCS proc. AS, 2013, 481 492 Singularity analysis via the iterated kernel method Stephen Melczer and Marni Mishna Department of Mathematics, Simon Fraser University, Burnaby,
More informationMath Fall 05 - Lectures notes # 1 - Aug 29 (Monday)
Math 110 - Fall 05 - Lectures notes # 1 - Aug 29 (Monday) Name, class, URL (www.cs.berkeley.edu/~demmel/ma110) on board See Barbara Peavy in 967 Evans Hall for all enrollment issues. All course material
More informationOR MSc Maths Revision Course
OR MSc Maths Revision Course Tom Byrne School of Mathematics University of Edinburgh t.m.byrne@sms.ed.ac.uk 15 September 2017 General Information Today JCMB Lecture Theatre A, 09:30-12:30 Mathematics revision
More informationE[X n ]= dn dt n M X(t). ). What is the mgf? Solution. Found this the other day in the Kernel matching exercise: 1 M X (t) =
Chapter 7 Generating functions Definition 7.. Let X be a random variable. The moment generating function is given by M X (t) =E[e tx ], provided that the expectation exists for t in some neighborhood of
More informationCALCULUS JIA-MING (FRANK) LIOU
CALCULUS JIA-MING (FRANK) LIOU Abstract. Contents. Power Series.. Polynomials and Formal Power Series.2. Radius of Convergence 2.3. Derivative and Antiderivative of Power Series 4.4. Power Series Expansion
More informationChapter 1 Divide and Conquer Polynomial Multiplication Algorithm Theory WS 2015/16 Fabian Kuhn
Chapter 1 Divide and Conquer Polynomial Multiplication Algorithm Theory WS 2015/16 Fabian Kuhn Formulation of the D&C principle Divide-and-conquer method for solving a problem instance of size n: 1. Divide
More informationMATH 1372, SECTION 33, MIDTERM 3 REVIEW ANSWERS
MATH 1372, SECTION 33, MIDTERM 3 REVIEW ANSWERS 1. We have one theorem whose conclusion says an alternating series converges. We have another theorem whose conclusion says an alternating series diverges.
More informationConsidering our result for the sum and product of analytic functions, this means that for (a 0, a 1,..., a N ) C N+1, the polynomial.
Lecture 3 Usual complex functions MATH-GA 245.00 Complex Variables Polynomials. Construction f : z z is analytic on all of C since its real and imaginary parts satisfy the Cauchy-Riemann relations and
More informationDIFFERENTIATION RULES
3 DIFFERENTIATION RULES DIFFERENTIATION RULES We have: Seen how to interpret derivatives as slopes and rates of change Seen how to estimate derivatives of functions given by tables of values Learned how
More informationIRREDUCIBILITY OF CLASSICAL POLYNOMIALS AND THEIR GENERALIZATIONS
IRREDUCIBILITY OF CLASSICAL POLYNOMIALS AND THEIR GENERALIZATIONS IRREDUCIBILITY OF CLASSICAL POLYNOMIALS AND THEIR GENERALIZATIONS by Michael Filaseta IRREDUCIBILITY OF CLASSICAL POLYNOMIALS AND THEIR
More informationChapter 11 - Sequences and Series
Calculus and Analytic Geometry II Chapter - Sequences and Series. Sequences Definition. A sequence is a list of numbers written in a definite order, We call a n the general term of the sequence. {a, a
More informationSection 9.7 and 9.10: Taylor Polynomials and Approximations/Taylor and Maclaurin Series
Section 9.7 and 9.10: Taylor Polynomials and Approximations/Taylor and Maclaurin Series Power Series for Functions We can create a Power Series (or polynomial series) that can approximate a function around
More informationOverview of Computer Algebra
Overview of Computer Algebra http://cocoa.dima.unige.it/ J. Abbott Universität Kassel J. Abbott Computer Algebra Basics Manchester, July 2018 1 / 12 Intro Characteristics of Computer Algebra or Symbolic
More informationFFTs in Graphics and Vision. Homogenous Polynomials and Irreducible Representations
FFTs in Graphics and Vision Homogenous Polynomials and Irreducible Representations 1 Outline The 2π Term in Assignment 1 Homogenous Polynomials Representations of Functions on the Unit-Circle Sub-Representations
More informationLattice Green s Functions of the Higher-Dimensional Face-Centered Cubic Lattices
Lattice Green s Functions of the Higher-Dimensional Face-Centered Cubic Lattices Christoph Koutschan Research Institute for Symbolic Computation (RISC) Johannes Kepler University, A-44 Linz, Austria August
More information