Numerical Integration with the double exponential method
|
|
- Christian Chambers
- 6 years ago
- Views:
Transcription
1 Numerical Integration with the double exponential method Pascal Molin Université Paris 7 December
2 Numerical integration Number-theory context, want fast quadrature high precision ( digits) rational recognition algebraic dependencies, LLL robust scheme little knowledge of the integrand yet can assume holomorphicity around path rigorous evaluation non-vanishing, inequalities algebraic
3 Digits matter Closed form is not always the best answer maple > int(x*sin(x)/(1+cos(x)^2), x = 0.. Pi); ( ) ( i dilog dilog ) ( i dilog 2 + dilog ) i 2 1 ( ) ( 1 + i + 2 ( + arctan 1 + ) ) 1 2 π 2 1 ( ( ) ) 1 ( arctan 2 1 π + i ln 1 + ) 2 π + 1/2 π 2
4 Digits matter Closed form is not always the best answer maple > int(x*sin(x)/(1+cos(x)^2), x = 0.. Pi); gp > intnum (x=0,pi,x* sin (x) /(1+ cos (x) ^2) ) time = 13 ms. % 1 =
5 Digits matter Closed form is not always the best answer maple > int(x*sin(x)/(1+cos(x)^2), x = 0.. Pi); gp > intnum (x=0,pi,x* sin (x) /(1+ cos (x) ^2) ) time = 13 ms. % 1 = PSLQ algorithm (integer relations with precomputed values) gp > Pi ^2/4 % 1 =
6 Classical quadrature schemes Trapeze sums better than rectangles! error = O(n 2 ) for n evaluations Simpson, Newton-Cotes degree d spline interpolation on each subinterval error = O(n 2d ) for nd evaluations omberg convergence acceleration error = O(n 2k ) for 2 k/2 n evaluations (k steps) Gauss-Legendre interpolation at n chosen points, exact on degree 2n 1 polynomials error = O(2 2n )
7 DE method [Takahasi&Mori,1973] trapezoidal method over for doubly-exponential decay g(x) Me e x h 1 D log D heuristic : D digits with O(D log D) points changes of variable to obtain DE decay g(x) Me e x f (x) Me x, x f (x) M(1 + x ) α, x t = sinh(u) x = sinh(t) x = e t e t t = tanh(u) 0 f (x) Me x, x f (x) M, x [ 1, 1]
8 Empirics [Bailey 2005] best quadrature scheme for number theory (speed + accuracy) fast gp > \ p100 gp > 4* intnum (x=0,1,x* log (x +1) ) time = 14 ms. % 2 = gp > 4* intnumgauss (x=0,1,x* log (x +1) ) time = 19 ms. % 3 =
9 Empirics [Bailey 2005] best quadrature scheme for number theory (speed + accuracy) fast gp > \p 1000 realprecision = 1001 significant digits (1000 digits displayed ) gp > intnum (x=0,1,x* log (x +1) ); time = 3,246 ms. gp > intnumgauss (x=0,1,x* log (x +1) ); *** Warning : increasing stack size to *** Warning : increasing stack size to *** Warning : increasing stack size to *** Warning : increasing stack size to time = 30,423 ms.
10 Empirics [Bailey 2005] best quadrature scheme for number theory (speed + accuracy) fast accurate gp > intnum (x=0,1,x* log (x)) time = 14 ms. % 1 = gp > intnumgauss (x=0,1,x* log (x)) time = 19 ms. % 2 =
11 Empirics [Bailey 2005] best quadrature scheme for number theory (speed + accuracy) fast accurate gp > intnum (x=0,1, sqrt (1 -x ^2) ) * 4/ Pi time = 5 ms. % 1 =
12 Empirics [Bailey 2005] best quadrature scheme for number theory (speed + accuracy) fast accurate gp > intnum (x=0,1, sqrt (1 -x ^2) ) * 4/ Pi time = 5 ms. % 1 = gp > intnumgauss (x=0,1, sqrt (1 -x ^2) ) * 4/ Pi time = 18 ms. % 1 =
13 Empirics [Bailey 2005] best quadrature scheme for number theory (speed + accuracy) fast accurate robust gp > intnum (x=0,1, log (x) ^2) time = 14 ms. % 1 =
14 Empirics [Bailey 2005] best quadrature scheme for number theory (speed + accuracy) fast accurate robust gp > intnum (x=0,1, log (x) ^2) time = 14 ms. % 1 = gp > intnumgauss (x=0,1, log (x) ^2) time = 20 ms. % 1 =
15 Empirics [Bailey 2005] best quadrature scheme for number theory (speed + accuracy) fast accurate robust gp > intnum (x = 0, [oo, -I], x ^2* sin (x)) time = 29 ms. % 1 =
16 Bugs? singularities near path gp > \ p100 gp > intnum (x=-oo,oo,1/(1+ x ^2) ) / Pi time = 16 ms. % 2 =
17 Bugs? singularities near path gp > \ p100 gp > intnum (x=-oo,oo,1/(1+ x ^2) ) / Pi time = 16 ms. % 2 = gp > intnum (t=-oo,oo,1/(1+( t +10) ^2) ) / Pi time = 16 ms. % 1 =
18 Bugs? singularities near path + = lim X X X gp > intnum (x=[ -oo,2],[ oo,2], exp (-x ^2) ) / sqrt (Pi) time = 14 ms. % 1 =
19 Bugs? singularities near path + = lim X X X gp > intnum (x=[ -oo,2],[ oo,2], exp (-x ^2) ) / sqrt (Pi) time = 14 ms. % 1 = gp > exp ( -20^2) % 1 = E -174
20 Bugs? singularities near path + = lim X X X gp > intnum (x=[ -oo,2],[ oo,2], exp (-x ^2) ) / sqrt (Pi) time = 14 ms. % 1 = gp > intnum (x= -20,20, exp (-x ^2) ) / sqrt (Pi) time = 12 ms. % 1 =
21 Bugs? singularities near path + = lim X X X gp > intnum (x=[ -oo,2],[ oo,2], exp (-x ^2) ) / sqrt (Pi) time = 14 ms. % 1 = gp > intnum (x= -50,50, exp (-x ^2) ) / sqrt (Pi) time = 12 ms. % 1 =
22 Goal : rigorous convergence Theorem Let f : [ 1, 1] C such that f has holomorphic continuation to B(0, 2) = { z < 2}. Then 1 1 f (x) dx for x k, w k given by where h = n k= n w k f (x k ) x k = tanh(sinh(kh)) and w k = log(n+2) n+2+2/7 log(n). 8 sup f 2 10n log(n+2) (1) B(0,2) h cosh(kh) cosh(sinh(kh)) 2 (2)
23 Comparison with Gauss Theorem Let f : [ 1, 1] C such that f has holomorphic continuation to B(0, 2) = { z < 2}. Then 1 1 f (x) dx for x k, w k given by n k= n ( w k f (x k ) π sup f D(0,2) ) 2 2n (3) P n (x k ) = 0 and w k = where P n are Legendre polynomials { P 0 (x) = 0, P 1 (x) = 1 2 (n + 1)P n(x k )P n+1 (x k ) (k + 1)P k+1 (x) = (2k + 1)xP k (x) kp k 1 (x) (4) (5)
24 Comparison with Gauss Integral of usual functions to D digits. method Gauss DE # points D D log D cost points D 3 log D D 2 log D cost integral D 2 log 2 D D 2 log 3 D many nice integrals single nice integral unknown behaviour
25 Theory
26 Poisson formula g(x) = O(x 2 ), g C 2 (), ĝ(x ) = e 2iπXt g(t) dt h k >n g(kh) + h } {{ } 1 g(k) = ĝ(k) k Z k Z n k= n g(kh) = g + k Z ĝ( k h ) }{{} 2 nh h nh
27 Poisson formula g(x) = O(x 2 ), g C 2 (), ĝ(x ) = e 2iπXt g(t) dt h k >n g(kh) + h } {{ } 1 g(k) = ĝ(k) k Z k Z n k= n g(kh) = g + k Z ĝ( k h ) }{{} 2 Errors 1(nh) and 2( 1 h ). Summation interesting when both g and ĝ vanish quickly (exponentially).
28 g(x) = sin x x, ĝ(x) = 1 [ 1 2π, 1 2π ] k sin(k) k g(x) = e x2, ĝ(x) = e πx2 = Poisson summation sin(x) x dt 1(nh) + 2(h) e D n e D 1(nh) + 2(h) e D n D π g(x) = e a cosh(x), ĝ(x) = 2K i2πx (a) e π2 x x 1(nh) + 2(h) e D n D log(d) π 2
29 g(x) = sin x x, ĝ(x) = 1 [ 1 2π, 1 2π ] k sin(mk) k = g(x) = e x2, ĝ(x) = e πx2 sin(mx) x Poisson summation dt for m = 1, (nh) + 2(h) e D n e D 1(nh) + 2(h) e D n D π g(x) = e a cosh(x), ĝ(x) = 2K i2πx (a) e π2 x x 1(nh) + 2(h) e D n D log(d) π 2
30 Bad news I g and ĝ cannot have both arbitrary decay. Theorem (Uncertainty principle) Let g : with g = 0. { g(x) = O(e α 1 x β1 ) If ĝ(x) = O(e α 2 x β2 ) then 1 β β 2 1. In particular e D n D 1 β β 2. Optimal case : the gaussian g(x) = exp( σx 2 ).
31 Paley-Wiener theory The decay of ĝ corresponds to the regularity of g : ĝ has finite support g entire of order 1 ĝ has more than exponential decay g entire g holomorphic on a strip + i[ t, t] ĝ(x) = O(e 2πtx )
32 Control of ĝ Theorem If g : C has holomorphic continuation to a strip τ = + i[ τ, τ] such that g( iτ) 1 + g( + iτ) 1 = M 2 (τ) < ; g(x ± it) x ± 0 uniformly on t < τ ; then ĝ(x) M 2 (τ)e 2πτ x.
33 Control of ĝ Theorem If g : C has holomorphic continuation to a strip τ = + i[ τ, τ] such that g( iτ) 1 + g( + iτ) 1 = M 2 (τ) < ; g(x ± it) x ± 0 uniformly on t < τ ; then ĝ(x) M 2 (τ)e 2πτ x. e 2πxτ e 2iπxt g(t + iτ) e 2iπxt g(t) iτ
34 Control of ĝ Theorem If g : C has holomorphic continuation to a strip τ = + i[ τ, τ] such that g( iτ) 1 + g( + iτ) 1 = M 2 (τ) < ; g(x ± it) x ± 0 uniformly on t < τ ; then ĝ(x) M 2 (τ)e 2πτ x. Goal : find best decay for g such that it remains holomorphic on a strip.
35 Bad news II Théorème If g : C satisfies g(x) M 1 e αe β x on (DE) ; g has bounded holomorphic continuation on a strip τ with βτ > π 2 then g = 0. In particular cannot expect more than DE decay for g for DE functions : e D n D log(d) π 2
36 Good news DE decay extends near the axis. Theorem (Phrägmen-Lindelöf) If g is holomorphic on a strip τ and satisfies g = O(e αe β x ) sur ; g M sur τ ; then for all t < τ, g(x ± it) Me α te β x, with α t = α(cos(βt) sin(βt) tan(βτ) ). easy to control ĝ from DE hypothesis
37 Main theorem Let g : such that 1 g has analytic continuation to a strip τ = + i] τ, τ[ ; 2 g(x) M 1 e αe β x on with α, β > 0 ; 3 g(z) M 2 on τ ; then for any D > 0 there are explicit values n, h with n (D + log M 2) log(d + log M 1 ) 2πτβ and h 2πτ D + log M 2 such that g h n k= n g(kh) e D.
38 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 τ τ iτ I ϕ( τ ) τ f holomorphic on τ.
39 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 f (x) M 1 e α x β τ iτ ϕ(t) = sinh(t) = et e t ϕ (t) = cosh(t) 2
40 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 τ sin(τ) τ τ iτ ϕ(t) = sinh(t) = et e t ϕ (t) = cosh(t) 2
41 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 τ f (x) M 1 e α x β sin(τ) τ f (z) M 2 e A z γ τ iτ ϕ(t) = sinh(sinh(t)) ϕ (t) = cosh(t) cosh(sinh(t))
42 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 τ iτ f (x) M1 1+ x α ϕ(t) = sinh(sinh(t)) ϕ (t) = cosh(t) cosh(sinh(t))
43 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 τ iτ ϕ(t) = tanh(λ sinh(t)) ϕ (t) = λ cosh(t) cosh 2 (λ sinh(t))
44 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 τ Ym(τ) M 2 f (z) (1+ z ) 1+γ τ iτ X m(τ) f (x) M1 (1+ x ) α ϕ(t) = tanh(λ sinh(t)) ϕ (t) = λ cosh(t) cosh 2 (λ sinh(t))
45 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 a f M 1 b τ iτ ϕ(t) = e t αe βt ϕ (t) = (1 + αβe βt )ϕ(t)
46 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 τ iτ ϕ(t) = e t αe βt ϕ (t) = (1 + αβe βt )ϕ(t)
47 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 Y m(τ) f M 2 f M 1 a b X m(τ) τ iτ
48 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 τ iτ 0 f M 1 e αxβ
49 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 τ iτ
50 DE method : regularity hypothesis Problem : DE setting : f, 1, 2 g = f (ϕ).ϕ, 1, 2 τ τ f (z) M 2 e λ z 0 τ iτ f M 1 e αxβ Xm
51 Summary τ Ym(τ) M 2 f (z) (1+ z ) 1+γ f (x) M 1 e α x β sin(τ) τ f (z) M 2 e A z γ τ X m(τ) f (x) M1 (1+ x ) α τ f (z) M 2 e λ z f M 2 Ym(τ) τ 0 f M 1 a f M 1 e αxβ b Xm(τ) Xm Under these hypothesis, DE method rigorously evaluates f to D digits using O( D log D 2πτ ) points.
52 Use cases
53 A list of integrals z dz sh(sh(t)) poles / 1+z 2 τ τ
54 Well placed poles
55 A list of integrals z dz sh(sh(t)) poles / 1+z 2 τ dz sh(sh(t)) pole 1+(z 10) τ, τ > τ
56 Small tau value
57 Small tau value? \ p1000 realprecision = 1001 significant digits (1000 digits displayed )? f(z) =1/(1+ z ^2) ;? D =1000* log (10) ;? show ( asinh ( asinh (10+ I))) % 6 = * I? tau =0.03; h= pas_h_shsh (tau,1,1,d); show (h) % 7 = E -5? n= ceil ( asinh ( asinh (2*10^ ) )/h); show (n) %8 = ? Pi - integration_shsh (z->f(z -10),h,n,1) % 9 = E -1001
58 A list of integrals z dz sh(sh(t)) poles / 1+z 2 τ dz sh(sh(t)) pole 1+(z 10) τ, τ > dz 100 th( π 1+z 2 2 sh(t)) pole τ, τ > 0.06 τ
59 A list of integrals z dz sh(sh(t)) poles / 1+z 2 τ dz sh(sh(t)) pole 1+(z 10) τ, τ > dz 100 th( π 1+z 2 2 sh(t)) pole τ, τ > 0.06 e z2 dz sh(t) pole 1+e τ, τ iz 1 τ
60 Line of poles
61 A list of integrals z dz sh(sh(t)) poles / 1+z 2 τ dz sh(sh(t)) pole 1+(z 10) τ, τ > dz 100 th( π 1+z 2 2 sh(t)) pole τ, τ > 0.06 e z2 dz sh(t) pole 1+e τ, τ iz 1 Γ(1 + iz) dz sh(t) Γ for τ < 0 τ
62 A list of integrals z dz sh(sh(t)) poles / 1+z 2 τ dz sh(sh(t)) pole 1+(z 10) τ, τ > dz 100 th( π 1+z 2 2 sh(t)) pole τ, τ > 0.06 e z2 dz sh(t) pole 1+e τ, τ iz 1 Γ(1 + iz) dz sh(t) Γ for τ < 0 e ch(z)+i ch(2z) dz t huge L 1 ( ± iτ) sin(z) dz z sh(sh(t)) sin(ie t ) explodes τ
63 Extensions
64 Meromorphic functions If has a pole ρ τ, with residue r ρ : ĝ(x ) = ±iτ For f ĝ( k h ) 2iπ k Z ρ Z τ where ε z is the sign of Im z. e 2iπτXt g(t) dt + 2iπe 2iπX ρ r ρ ϕ 1 (ρ) ε z r ρ e ε z 2iπz/h 1 2M 1 e 2πτ/h 1
65 Taking residues into account? \ p1000 realprecision = 1001 significant digits (1000 digits displayed )? f(z) =1/(1+ z ^2) ;? tau =Pi /2.2;? h =2* Pi*tau /(1000* log (10) + log (2/ cos ( tau ))); show (h) % 6 = ? n= ceil ( asinh ( asinh (2*10^1000) )/h) %7 = 2168? Pi - integration_shsh (f,h,n,1) %8 = 0.E -1001? diff =Pi - integration_shsh (z->f(z -15),h,n,1) ; show ( diff ) % 9 = E -14? ho =[15+I,15 -I]; es =[ -I/2,I /2];? diff + poles_shsh (ho,es,h) % 11 = E E -1014* I
66 A list of integrals z dz sh(sh(t)) poles / 1+z 2 τ dz sh(sh(t)) pole 1+(z 10) τ, τ > dz 100 th( π 1+z 2 2 sh(t)) pole τ, τ > 0.06 e z2 dz sh(t) pole 1+e τ, τ iz 1 Γ(1 + iz) dz sh(t) Γ for τ < 0 e ch(z)+i ch(2z) dz t huge L 1 ( ± iτ) sin(z) dz z sh(sh(t)) sin(ie t ) explodes τ
67 Shifting path 0 0 θ δ gamma function on vertical lines
68 Shifting path θ gamma function on vertical lines δ integrate on the middle of a nice strip
69 A list of integrals z dz sh(sh(t)) poles / 1+z 2 τ dz sh(sh(t)) pole 1+(z 10) τ, τ > dz 100 th( π 1+z 2 2 sh(t)) pole τ, τ > 0.06 e z2 dz sh(t) pole 1+e τ, τ iz 1 Γ(1 + iz) dz sh(t) Γ for τ < 0 e ch(z)+i ch(2z) dz t huge L 1 ( ± iτ) sin(z) dz z sh(sh(t)) sin(ie t ) explodes τ
70 Examples
71 erfc(x) = i e z2 π e x2 z ix dz Error function
72 ( erfc(x) = hx 1 π e x2 erfc(x) = i π e x2 x k 1 comparison with MPF e k2 h 2 k 2 h 2 + x 2 Error function e z2 z ix dz ) e 2πx/h + O(e π2 /h 2 ) gp/mpfr 4 MPF 50 times faster precision integration 400 times faster
73 Incomplete gamma function Γ inc (s, x) = x ts e t dt t for s C and x /. X 0
74 Incomplete gamma function Γ inc (s, x) = x ts e t dt t for s C and x /. X τ 0
75 Incomplete gamma function Γ inc (s, x) = x ts e t dt t for s C and x /. X τ 0
76 Incomplete gamma function Γ inc (s, x) = x ts e t dt t for s C and x /. X 0
77 Incomplete gamma function Γ inc (s, x) = x ts e t dt t for s C and x /. t s 1 e t M(X, τ, s) explicite τ X θ 0 Γ inc (s, X ) to prec D using n 2D log D π 2 evaluations.
78 Conclusion DE integration is not always the best method (better convergence for Gauss, provided big precomputations and well known integrand) still very easy to implement and nice convergence rigorous integration possible with basic bounds on the integrand otherwise efficient heuristic procedures (guess the error term from first iterations)
79 eferences Hidetosi Takahasi and Masatake Mori. Double exponential formulas for numerical integration. Kyoto University. esearch Institute for Mathematical Sciences. Publications, 9 : , Masatake Mori. Discovery of the double exponential transformation and its developments. Publ. es. Inst. Math. Sci., 41(4) : , David H. Bailey, Karthik Jeyabalan, and Xiaoye S. Li. A comparison of three high-precision quadrature schemes. Experiment. Math., 14(3) : , Pascal Molin. Numerical integration and L functions computations. Theses, Université Sciences et Technologies - Bordeaux I, October 2010.
Experimental 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 informationComplex Analysis, Stein and Shakarchi The Fourier Transform
Complex Analysis, Stein and Shakarchi Chapter 4 The Fourier Transform Yung-Hsiang Huang 2017.11.05 1 Exercises 1. Suppose f L 1 (), and f 0. Show that f 0. emark 1. This proof is observed by Newmann (published
More informationSummation of Series and Gaussian Quadratures
Summation of Series Gaussian Quadratures GRADIMIR V. MILOVANOVIĆ Dedicated to Walter Gautschi on the occasion of his 65th birthday Abstract. In 985, Gautschi the author constructed Gaussian quadrature
More informationApril 15 Math 2335 sec 001 Spring 2014
April 15 Math 2335 sec 001 Spring 2014 Trapezoid and Simpson s Rules I(f ) = b a f (x) dx Trapezoid Rule: If [a, b] is divided into n equally spaced subintervals of length h = (b a)/n, then I(f ) T n (f
More informationNUMERICAL QUADRATURE: THEORY AND COMPUTATION
NUMERICAL QUADRATURE: THEORY AND COMPUTATION by Lingyun Ye Submitted in partial fulfillment of the requirements for the degree of Master of Computer Science at Dalhousie University Halifax, Nova Scotia
More informationSection 6.6 Gaussian Quadrature
Section 6.6 Gaussian Quadrature Key Terms: Method of undetermined coefficients Nonlinear systems Gaussian quadrature Error Legendre polynomials Inner product Adapted from http://pathfinder.scar.utoronto.ca/~dyer/csca57/book_p/node44.html
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationA Numerical Integration Formula Based on the Bessel Functions
Publ. RIMS, Kyoto Univ. 41 (25), 949 97 A Numerical Integration Formula Based on the Bessel Functions By Hidenori Ogata Abstract In this paper, we discuss the properties of a quadrature formula with the
More informationNumerical integration in arbitrary-precision ball arithmetic
/ 24 Numerical integration in arbitrary-precision ball arithmetic Fredrik Johansson (LFANT, Bordeaux) Journées FastRelax, Inria Sophia Antipolis 7 June 208 Numerical integration in Arb 2 / 24 New code
More informationPOISSON SUMMATION AND PERIODIZATION
POISSON SUMMATION AND PERIODIZATION PO-LAM YUNG We give some heuristics for the Poisson summation formula via periodization, and provide an alternative proof that is slightly more motivated.. Some heuristics
More informationz b k P k p k (z), (z a) f (n 1) (a) 2 (n 1)! (z a)n 1 +f n (z)(z a) n, where f n (z) = 1 C
. Representations of Meromorphic Functions There are two natural ways to represent a rational function. One is to express it as a quotient of two polynomials, the other is to use partial fractions. The
More informationMath 715 Homework 1 Solutions
. [arrier, Krook and Pearson Section 2- Exercise ] Show that no purely real function can be analytic, unless it is a constant. onsider a function f(z) = u(x, y) + iv(x, y) where z = x + iy and where u
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 informationComputing Fresnel Integrals via Modified Trapezium Rules
School of Mathematical and Physical Sciences Department of Mathematics and Statistics Preprint MPS-01-0 17 September 01 Computing Fresnel Integrals via Modified Trapezium Rules by Mohammad Alazah, Simon
More informationMath 341: Probability Seventeenth Lecture (11/10/09)
Math 341: Probability Seventeenth Lecture (11/10/09) Steven J Miller Williams College Steven.J.Miller@williams.edu http://www.williams.edu/go/math/sjmiller/ public html/341/ Bronfman Science Center Williams
More informationA Comparison of Three High-Precision Quadrature Schemes
A Comparison of Three High-Precision Quadrature Schemes David H. Bailey, Karthik Jeyabalan, and Xiaoye S. Li CONTENTS 1. Introduction 2. The ARPREC Software 3. The Three quadrature Schemes 4. The Euler-Maclaurin
More information1 Discussion on multi-valued functions
Week 3 notes, Math 7651 1 Discussion on multi-valued functions Log function : Note that if z is written in its polar representation: z = r e iθ, where r = z and θ = arg z, then log z log r + i θ + 2inπ
More information17 The functional equation
18.785 Number theory I Fall 16 Lecture #17 11/8/16 17 The functional equation In the previous lecture we proved that the iemann zeta function ζ(s) has an Euler product and an analytic continuation to the
More informationDavid A. Stephens Department of Mathematics and Statistics McGill University. October 28, 2006
556: MATHEMATICAL STATISTICS I COMPUTING THE HYPEBOLIC SECANT DISTIBUTION CHAACTEISTIC FUNCTION David A. Stephens Department of Mathematics and Statistics McGill University October 8, 6 Abstract We give
More informationFast Evaluation of Special Functions by the Modified Trapezium Rule
Fast Evaluation of Special Functions by the Modified Trapezium Rule MOHAMMAD AL AZAH Department of Mathematics and Statistics University of Reading This dissertation is submitted for the degree of Doctor
More informationMATH 311: COMPLEX ANALYSIS CONTOUR INTEGRALS LECTURE
MATH 3: COMPLEX ANALYSIS CONTOUR INTEGRALS LECTURE Recall the Residue Theorem: Let be a simple closed loop, traversed counterclockwise. Let f be a function that is analytic on and meromorphic inside. Then
More informationExercises involving elementary functions
017:11:0:16:4:09 c M K Warby MA3614 Complex variable methods and applications 1 Exercises involving elementary functions 1 This question was in the class test in 016/7 and was worth 5 marks a) Let z +
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 informationIntroduction to Numerical Analysis
Introduction to Numerical Analysis S. Baskar and S. Sivaji Ganesh Department of Mathematics Indian Institute of Technology Bombay Powai, Mumbai 400 076. Introduction to Numerical Analysis Lecture Notes
More informationComplex Analysis, Stein and Shakarchi Meromorphic Functions and the Logarithm
Complex Analysis, Stein and Shakarchi Chapter 3 Meromorphic Functions and the Logarithm Yung-Hsiang Huang 217.11.5 Exercises 1. From the identity sin πz = eiπz e iπz 2i, it s easy to show its zeros are
More informationMath 3313: Differential Equations Second-order ordinary differential equations
Math 3313: Differential Equations Second-order ordinary differential equations Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Mass-spring & Newton s 2nd law Properties
More informationPrinciple of Mathematical Induction
Advanced Calculus I. Math 451, Fall 2016, Prof. Vershynin Principle of Mathematical Induction 1. Prove that 1 + 2 + + n = 1 n(n + 1) for all n N. 2 2. Prove that 1 2 + 2 2 + + n 2 = 1 n(n + 1)(2n + 1)
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationb n x n + b n 1 x n b 1 x + b 0
Math Partial Fractions Stewart 7.4 Integrating basic rational functions. For a function f(x), we have examined several algebraic methods for finding its indefinite integral (antiderivative) F (x) = f(x)
More informationOn the Use of Conformal Maps for the Acceleration of Convergence of the Trapezoidal Rule and Sinc Numerical Methods
On the Use of Conformal Maps for the Acceleration of Convergence of the Trapezoidal Rule and Sinc Numerical Methods Richard Mikael Slevinsky and Sheehan Olver Mathematical Institute, University of Oxford
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 informationarxiv: v2 [cs.na] 2 Jan 2016
On Approximating Univariate NP-ard Integrals Ohad Asor and Avishy Carmi 1 Jan 2016 arxiv:1512.08716v2 cs.na] 2 Jan 2016 Abstract 1 Approximating I #PART = n 0 k=1 cosx kπt)dt to within an accuracy of 2
More informationSolutions to Complex Analysis Prelims Ben Strasser
Solutions to Complex Analysis Prelims Ben Strasser In preparation for the complex analysis prelim, I typed up solutions to some old exams. This document includes complete solutions to both exams in 23,
More informationSuggested Homework Solutions
Suggested Homework Solutions Chapter Fourteen Section #9: Real and Imaginary parts of /z: z = x + iy = x + iy x iy ( ) x iy = x #9: Real and Imaginary parts of ln z: + i ( y ) ln z = ln(re iθ ) = ln r
More informationProblems for MATH-6300 Complex Analysis
Problems for MATH-63 Complex Analysis Gregor Kovačič December, 7 This list will change as the semester goes on. Please make sure you always have the newest version of it.. Prove the following Theorem For
More informationFinal Year M.Sc., Degree Examinations
QP CODE 569 Page No Final Year MSc, Degree Examinations September / October 5 (Directorate of Distance Education) MATHEMATICS Paper PM 5: DPB 5: COMPLEX ANALYSIS Time: 3hrs] [Max Marks: 7/8 Instructions
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationMS 3011 Exercises. December 11, 2013
MS 3011 Exercises December 11, 2013 The exercises are divided into (A) easy (B) medium and (C) hard. If you are particularly interested I also have some projects at the end which will deepen your understanding
More informationMath Final Exam.
Math 106 - Final Exam. This is a closed book exam. No calculators are allowed. The exam consists of 8 questions worth 100 points. Good luck! Name: Acknowledgment and acceptance of honor code: Signature:
More informationPart IB. Further Analysis. Year
Year 2004 2003 2002 2001 10 2004 2/I/4E Let τ be the topology on N consisting of the empty set and all sets X N such that N \ X is finite. Let σ be the usual topology on R, and let ρ be the topology on
More informationMath 259: Introduction to Analytic Number Theory Functions of finite order: product formula and logarithmic derivative
Math 259: Introduction to Analytic Number Theory Functions of finite order: product formula and logarithmic derivative This chapter is another review of standard material in complex analysis. See for instance
More informationMath 259: Introduction to Analytic Number Theory Functions of finite order: product formula and logarithmic derivative
Math 259: Introduction to Analytic Number Theory Functions of finite order: product formula and logarithmic derivative This chapter is another review of standard material in complex analysis. See for instance
More informationOn the number of ways of writing t as a product of factorials
On the number of ways of writing t as a product of factorials Daniel M. Kane December 3, 005 Abstract Let N 0 denote the set of non-negative integers. In this paper we prove that lim sup n, m N 0 : n!m!
More informationHigh-speed high-accuracy computation of an infinite integral with unbounded and oscillated integrand. Takuya OOURA. February 2012
RIMS-1741 High-speed high-accuracy computation of an infinite integral with unbounded and oscillated integrand By Takuya OOURA February 2012 RESEARCH INSTITUTE FOR MATHEMATICAL SCIENCES KYOTO UNIVERSITY,
More information2 Write down the range of values of α (real) or β (complex) for which the following integrals converge. (i) e z2 dz where {γ : z = se iα, < s < }
Mathematical Tripos Part II Michaelmas term 2007 Further Complex Methods, Examples sheet Dr S.T.C. Siklos Comments and corrections: e-mail to stcs@cam. Sheet with commentary available for supervisors.
More informationExamination paper for TMA4215 Numerical Mathematics
Department of Mathematical Sciences Examination paper for TMA425 Numerical Mathematics Academic contact during examination: Trond Kvamsdal Phone: 93058702 Examination date: 6th of December 207 Examination
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory
More informationThe Calculus of Residues
hapter 7 The alculus of Residues If fz) has a pole of order m at z = z, it can be written as Eq. 6.7), or fz) = φz) = a z z ) + a z z ) +... + a m z z ) m, 7.) where φz) is analytic in the neighborhood
More informationHigher Order Linear Equations
C H A P T E R 4 Higher Order Linear Equations 4.1 1. The differential equation is in standard form. Its coefficients, as well as the function g(t) = t, are continuous everywhere. Hence solutions are valid
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationPhysics 307. Mathematical Physics. Luis Anchordoqui. Wednesday, August 31, 16
Physics 307 Mathematical Physics Luis Anchordoqui 1 Bibliography L. A. Anchordoqui and T. C. Paul, ``Mathematical Models of Physics Problems (Nova Publishers, 2013) G. F. D. Duff and D. Naylor, ``Differential
More informationConformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.
Lent 29 COMPLEX METHODS G. Taylor A star means optional and not necessarily harder. Conformal maps. (i) Let f(z) = az + b, with ad bc. Where in C is f conformal? cz + d (ii) Let f(z) = z +. What are the
More informationTitle Fast computation of Goursat s infin accuracy Author(s) Ooura, Takuya Citation Journal of Computational and Applie 249: 1-8 Issue Date 2013-09 URL http://hdl.handle.net/2433/173114 Right 2013 Elsevier
More informationMath 185 Fall 2015, Sample Final Exam Solutions
Math 185 Fall 2015, Sample Final Exam Solutions Nikhil Srivastava December 12, 2015 1. True or false: (a) If f is analytic in the annulus A = {z : 1 < z < 2} then there exist functions g and h such that
More informationGENG2140, S2, 2012 Week 7: Curve fitting
GENG2140, S2, 2012 Week 7: Curve fitting Curve fitting is the process of constructing a curve, or mathematical function, f(x) that has the best fit to a series of data points Involves fitting lines and
More informationIII. Consequences of Cauchy s Theorem
MTH6 Complex Analysis 2009-0 Lecture Notes c Shaun Bullett 2009 III. Consequences of Cauchy s Theorem. Cauchy s formulae. Cauchy s Integral Formula Let f be holomorphic on and everywhere inside a simple
More information8.3 Partial Fraction Decomposition
8.3 partial fraction decomposition 575 8.3 Partial Fraction Decomposition Rational functions (polynomials divided by polynomials) and their integrals play important roles in mathematics and applications,
More informationAnalysis Comprehensive Exam, January 2011 Instructions: Do as many problems as you can. You should attempt to answer completely some questions in both
Analysis Comprehensive Exam, January 2011 Instructions: Do as many problems as you can. You should attempt to answer completely some questions in both real and complex analysis. You have 3 hours. Real
More informationSimultaneous Gaussian quadrature for Angelesco systems
for Angelesco systems 1 KU Leuven, Belgium SANUM March 22, 2016 1 Joint work with Doron Lubinsky Introduced by C.F. Borges in 1994 Introduced by C.F. Borges in 1994 (goes back to Angelesco 1918). Introduced
More informationIntegration. Topic: Trapezoidal Rule. Major: General Engineering. Author: Autar Kaw, Charlie Barker.
Integration Topic: Trapezoidal Rule Major: General Engineering Author: Autar Kaw, Charlie Barker 1 What is Integration Integration: The process of measuring the area under a function plotted on a graph.
More informationSecond Midterm Exam Name: Practice Problems March 10, 2015
Math 160 1. Treibergs Second Midterm Exam Name: Practice Problems March 10, 015 1. Determine the singular points of the function and state why the function is analytic everywhere else: z 1 fz) = z + 1)z
More informationAccurate Multiple-Precision Gauss-Legendre Quadrature
Accurate Multiple-Precision Gauss-Legendre Quadrature Laurent Fousse Université Henri-Poincaré Nancy 1 laurent@komite.net Abstract Numerical integration is an operation that is frequently available in
More informationNumerical integration and differentiation. Unit IV. Numerical Integration and Differentiation. Plan of attack. Numerical integration.
Unit IV Numerical Integration and Differentiation Numerical integration and differentiation quadrature classical formulas for equally spaced nodes improper integrals Gaussian quadrature and orthogonal
More informationarxiv: v2 [math.na] 23 Jul 2013
Error control of a numerical formula for the Fourier transform by Ooura s continuous Euler transform and fractional FFT Ken ichiro Tanaka a a School of Systems Information Science, Future University Hakodate,
More informationFolland: Real Analysis, Chapter 8 Sébastien Picard
Folland: Real Analysis, Chapter 8 Sébastien Picard Problem 8.3 Let η(t) = e /t for t >, η(t) = for t. a. For k N and t >, η (k) (t) = P k (/t)e /t where P k is a polynomial of degree 2k. b. η (k) () exists
More information1 + z 1 x (2x y)e x2 xy. xe x2 xy. x x3 e x, lim x log(x). (3 + i) 2 17i + 1. = 1 2e + e 2 = cosh(1) 1 + i, 2 + 3i, 13 exp i arctan
Complex Analysis I MT333P Problems/Homework Recommended Reading: Bak Newman: Complex Analysis Springer Conway: Functions of One Complex Variable Springer Ahlfors: Complex Analysis McGraw-Hill Jaenich:
More informationA REVIEW OF RESIDUES AND INTEGRATION A PROCEDURAL APPROACH
A REVIEW OF RESIDUES AND INTEGRATION A PROEDURAL APPROAH ANDREW ARHIBALD 1. Introduction When working with complex functions, it is best to understand exactly how they work. Of course, complex functions
More informationQualifying Exam Complex Analysis (Math 530) January 2019
Qualifying Exam Complex Analysis (Math 53) January 219 1. Let D be a domain. A function f : D C is antiholomorphic if for every z D the limit f(z + h) f(z) lim h h exists. Write f(z) = f(x + iy) = u(x,
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More information4 Uniform convergence
4 Uniform convergence In the last few sections we have seen several functions which have been defined via series or integrals. We now want to develop tools that will allow us to show that these functions
More informationComputation of the omega function
Computation of the omega function David H. Bailey http://www.davidhbailey.com Lawrence Berkeley National Laboratory (retired) University of California, Davis, Department of Computer Science 1 / 17 Mordell-Tornheim-Witten
More informationA NOTE ON RATIONAL OPERATOR MONOTONE FUNCTIONS. Masaru Nagisa. Received May 19, 2014 ; revised April 10, (Ax, x) 0 for all x C n.
Scientiae Mathematicae Japonicae Online, e-014, 145 15 145 A NOTE ON RATIONAL OPERATOR MONOTONE FUNCTIONS Masaru Nagisa Received May 19, 014 ; revised April 10, 014 Abstract. Let f be oeprator monotone
More informationUsing approximate functional equations to build L functions
Using approximate functional equations to build L functions Pascal Molin Université Paris 7 Clermont-Ferrand 20 juin 2017 Example : elliptic curves Consider an elliptic curve E /Q of conductor N and root
More informationNumerical Methods I: Numerical Integration/Quadrature
1/20 Numerical Methods I: Numerical Integration/Quadrature Georg Stadler Courant Institute, NYU stadler@cims.nyu.edu November 30, 2017 umerical integration 2/20 We want to approximate the definite integral
More informationMATH 185: COMPLEX ANALYSIS FALL 2009/10 PROBLEM SET 9 SOLUTIONS. and g b (z) = eπz/2 1
MATH 85: COMPLEX ANALYSIS FALL 2009/0 PROBLEM SET 9 SOLUTIONS. Consider the functions defined y g a (z) = eiπz/2 e iπz/2 + Show that g a maps the set to D(0, ) while g maps the set and g (z) = eπz/2 e
More informationOn the Bilateral Laplace Transform of the positive even functions and proof of the Riemann Hypothesis. Seong Won Cha Ph.D.
On the Bilateral Laplace Transform of the positive even functions and proof of the Riemann Hypothesis Seong Won Cha Ph.D. Seongwon.cha@gmail.com Remark This is not an official paper, rather a brief report.
More informationPHYS-4007/5007: Computational Physics Course Lecture Notes Appendix G
PHYS-4007/5007: Computational Physics Course Lecture Notes Appendix G Dr. Donald G. Luttermoser East Tennessee State University Version 7.0 Abstract These class notes are designed for use of the instructor
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 informationPart IB. Complex Analysis. Year
Part IB Complex Analysis Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2018 Paper 1, Section I 2A Complex Analysis or Complex Methods 7 (a) Show that w = log(z) is a conformal
More informationAnalyse 3 NA, FINAL EXAM. * Monday, January 8, 2018, *
Analyse 3 NA, FINAL EXAM * Monday, January 8, 08, 4.00 7.00 * Motivate each answer with a computation or explanation. The maximum amount of points for this exam is 00. No calculators!. (Holomorphic functions)
More informationEvaluation of integrals
Evaluation of certain contour integrals: Type I Type I: Integrals of the form 2π F (cos θ, sin θ) dθ If we take z = e iθ, then cos θ = 1 (z + 1 ), sin θ = 1 (z 1 dz ) and dθ = 2 z 2i z iz. Substituting
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 informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 17, 2017 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More informationTopic 4 Notes Jeremy Orloff
Topic 4 Notes Jeremy Orloff 4 auchy s integral formula 4. Introduction auchy s theorem is a big theorem which we will use almost daily from here on out. Right away it will reveal a number of interesting
More informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
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 informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationDetermine for which real numbers s the series n>1 (log n)s /n converges, giving reasons for your answer.
Problem A. Determine for which real numbers s the series n> (log n)s /n converges, giving reasons for your answer. Solution: It converges for s < and diverges otherwise. To see this use the integral test,
More informationThe arithmetic geometric mean (Agm)
The arithmetic geometric mean () Pictures by David Lehavi This is part of exercise 5 question 4 in myc> calculus class Fall 200: a = c a n+ =5a n 2 lim n an = a = a an + bn a n+ = 2 b = b b n+ = a nb n
More informationEXTENDED RECIPROCAL ZETA FUNCTION AND AN ALTERNATE FORMULATION OF THE RIEMANN HYPOTHESIS. M. Aslam Chaudhry. Received May 18, 2007
Scientiae Mathematicae Japonicae Online, e-2009, 05 05 EXTENDED RECIPROCAL ZETA FUNCTION AND AN ALTERNATE FORMULATION OF THE RIEMANN HYPOTHESIS M. Aslam Chaudhry Received May 8, 2007 Abstract. We define
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationMTH 3102 Complex Variables Final Exam May 1, :30pm-5:30pm, Skurla Hall, Room 106
Name (Last name, First name): MTH 02 omplex Variables Final Exam May, 207 :0pm-5:0pm, Skurla Hall, Room 06 Exam Instructions: You have hour & 50 minutes to complete the exam. There are a total of problems.
More informationIn this chapter we study several functions that are useful in calculus and other areas of mathematics.
Calculus 5 7 Special functions In this chapter we study several functions that are useful in calculus and other areas of mathematics. 7. Hyperbolic trigonometric functions The functions we study in this
More informationSolutions to practice problems for the final
Solutions to practice problems for the final Holomorphicity, Cauchy-Riemann equations, and Cauchy-Goursat theorem 1. (a) Show that there is a holomorphic function on Ω = {z z > 2} whose derivative is z
More information1 Euler s idea: revisiting the infinitude of primes
8.785: Analytic Number Theory, MIT, spring 27 (K.S. Kedlaya) The prime number theorem Most of my handouts will come with exercises attached; see the web site for the due dates. (For example, these are
More informationMathematics of Physics and Engineering II: Homework answers You are encouraged to disagree with everything that follows. P R and i j k 2 1 1
Mathematics of Physics and Engineering II: Homework answers You are encouraged to disagree with everything that follows Homework () P Q = OQ OP =,,,, =,,, P R =,,, P S = a,, a () The vertex of the angle
More informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
More informationPSI Lectures on Complex Analysis
PSI Lectures on Complex Analysis Tibra Ali August 14, 14 Lecture 4 1 Evaluating integrals using the residue theorem ecall the residue theorem. If f (z) has singularities at z 1, z,..., z k which are enclosed
More information