Integral equations, eigenvalue, function interpolation
|
|
- Deborah Hicks
- 5 years ago
- Views:
Transcription
1 Integrl equtions, eigenvlue, function interpoltion Mrcin Chrząszcz Monte Crlo methods, 26 My, / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 1/19 19
2 Integrl equtions, introduction Fredholm integrl eqution of the second order: ϕ(x) = f(x) + K(x, y)ϕ(y)dy The f nd K re known functions K is clled kernel The CHALLENGE: find the ϕ tht obeys the bove equtions There re NO numericl tht cn solve this type of equtions! Different methods hve to be used depending on the f nd K functions The MC lgorithm: construct probbilistic lgorithm which hs n expected vlue the solution of the bove equtions There re mny wys to build this! We ssume tht the Neumnn series converges! 2 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 2/19 19
3 Integrl equtions, pproximtions The following steps pproximte the Fredholm eqution: ϕ 0(x) = 0, ϕ 1(x) = f(x) + K(x, y)ϕ 0(y)dy = f(x) ϕ 2(x) = f(x) + K(x, y)ϕ 1(y)dy = f(x) + K(x, y)f(y)dy ϕ 3(x) = f(x) + K(x, y)dy = f(x) + K(x, y)f(y)dy + K(x, y)k(y, z)f(z)dydz Now we put the following nottions: One gets: K (1) (x, y) = K(x, y) K (2) (x, y) = ϕ 3 (x) = f(x) + K (1) (x, y)f(y)dy + K(x, t)k(t, y)dt K (2) (x, y)f(y)dy 3 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 3/19 19
4 Integrl equtions, pproximtions Continuing this process: nd the n-th pproximtion: K (n) (x, y) = K(x, t)k (n 1) (t, y)dt ϕ n(x) = f(x) + K (1) (x, y)f(y)dy + K (2) (x, y)f(y)dy+ Now going with the Neumnn series: n : ϕ(x) = lim n ϕ n(x) = f(x) + + K (n) (x, y)f(y)dy n i=1 K (n) (x, y)f(y)dy The bove series converges only inside the squre: x, y b for: K(x, y) 2 dxdy < 1 4 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 4/19 19
5 Integrl equtions, lgorithm The rndom wlk of prticle hppens on the intervl (, b): In the t = 0 the prticle is in the position x 0 = x If the prticle t time t = n 1 is in the x n 1 position then in time t = n the position is: x n = x n 1 + ξ n The numbers ξ 1, ξ 2, re independent rndom numbers generted from ρ pdf The prticle stops the wlk once it reches the position or b The prticle life time is n when x n nd x n b The expected life time is given by the eqution: where: τ(x) = ρ 1 (x) + ρ q(x) = x [1 + τ(y)] ρ(y x)dy ρ(y)dy + ρ(y)dy b x is the probbility of prticle nnihiltion in the time t = 1 5 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 5/19 19
6 Integrl equtions, lgorithm 2 The bove cn be trnsformed: τ(x) = 1 + τ(y)ρ(x y)dy (1) Now if p(x) is the probbility tht the prticle in time t = 0 ws in position x gets nnihilted becuse it crosses the border The probbility obeys the nlogous eqution: where p(x) = ρ(x) + ρ(x) = x p(y)ρ(y x)dy (2) ρ(y)dy is the probbility of nnihilting the prticle in the first wlk For the functions τ nd ρ we got the integrl Fredholm eqution So the bove rndom wlk cn be be used to solve the Equtions 1 nd 2 6 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 6/19 19
7 Integrl equtions, lgorithm 3 The ρ(x) is the pdf of rndom vribles ξ n We observe the rndom wlk of the prticle The trjectory: γ = (x 0, x 1, x 2,, x n) This mens for t = 0, 1, 2, n 1 nd x n or x n b Additionlly we mrk: γ r = (x 0, x 1,, x r ), r n We defined rndom vrible: n S(x) = V (γ r)f(x r 1) where r=1 V (γ 0 ) = 1, V (γ r) = K(xr 1, xr) ρ(x r x r 1) V (γr 1) One cn prove tht E [S(x)] treted s function of x vrible is the solution to the integrl eqution 7 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 7/19 19
8 Integrl equtions, lgorithm 4 We define new rndom vrible: where ρ r (x) is defined s: ρ 1 (x) = ρ r (x) = x ρ(y)dy + c r(x) = + b x { V (γn r )f(x n r ) ρ(x n r, ) r n, 0, r > n ρ(y)dy, ρ(x 1 x)ρ(x 2 x 1 )ρ(x r 1 x r 2 )ρ 1 (x r 1 )dx 1 dx r 1 is the probbility tht the prticle tht is t given time in the x coordinte will survive r moments One cn prove tht E [c r (x)] treted s function of x vrible is the solution to the integrl eqution 8 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 8/19 19
9 Integrl equtions, generl remrk There is generl trick: Any integrl eqution cn be trnsformed to liner eqution using qudrtic form If done so one cn use the lgorithms form lecture 8 to solve it Bullet prove solution! 9 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 9/19 19
10 Eigenvlue problem The Eigenvlue problem is to find λ tht obeys the eqution: H x = λ x For simplicity we ssume there the biggest Eigenvlue is singulr nd it s rel The numericl method is bsiclly n itertive procedure to find the biggest Eigenvlue: We choose rndomly vector x 0 The m vector we choose ccordingly to formul: x m = H x m 1 /λ m where λ m is choose such tht n ( x m) j = 1 j=1 the ( x ) j is the j coordinte of the x vector, j = 1, 2, 3,, n The set λ m is converging to the lrgest Eigenvlue of the H mtrix 10 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 10/19 19
11 Eigenvlue problem From the bove we get: λ 1λ 2λ m( x j) = (H m n x 0) j; λ 1λ 2λ m = (H m x 0) j j=1 For big k nd m > k one gets: n j=1 (Hm x 0) j n j=1 (Hk x 0 ) j = λ k+1 λ k+2 λ m λ m k from which: λ [ n (Hm ] x j=1 0 ) m k 1 j n j=1 (Hk x 0 ) j This is the Eigenvlue estimtion corresponding to H m x 0 for sufficient lrge m 11 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 11/19 19
12 Eigenvlue problem, probbilistic model Let Q = (q ij), i, j = 1, 2,, n is the probbility mtrix: q ij 0, n q ij = 1 j=1 We construct rndom wlk on the set: {1, 2, n} ccordingly to the below rules: In the t = 0 the prticle is in rndomly chosen stte i 0 ccordingly to binned pdf: p j If in the moment t = n 1 the prticle is in i n 1 stte then in the next moment it goes to the stte i n with the probbility q in 1 j For γ = (i 0, i 1, ) trjectory we define rndom vrible: Now we do: W r (γ) = ( x ) i0 p i0 h i1 i 0 h i2 i 1 h i3 i 2 h ir i r 1 q i1 i 0 q i2 i 1 q i3 i 2 q iri r 1 E [W m (γ)] E [W k (γ)] λm k So to estimte the lrgest Eigenvlue: 12 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 12/19 19 [ ] 1
13 Function interpoltion Lets put f(x 1) = f 1, f(x 2) = f 2, which we know the functions The problem: clculte the f(p) for x 1 < p < x 2 From the interpoltion method we get: f(p) = p x 1 x 2 x 1 f 2 + x 2 p x 2 x 1 f 1 I m jet-lgged writing this so let me put: x 1 = 0 nd x 2 = 1: f(p) = (1 p)f 1 + pf 2 For 2-dim: where: f(p 1, p 2 ) = r 1 r 2 f(δ 1, δ 2 ) δ { 1 pi, δ i = 0 r i = p i, δ i = 1 the sum is over ll pirs (in this cse 4) 13 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 13/19 19
14 Function interpoltion For n-dim we get monstrous: f(p 1, p 2,, p n ) = r 1 r 2 r n f(δ 1,, δ n ) δ the sum is over ll combintions (δ 1,, δ n), where δ i = 0, 1 The bove sum is over 2 n terms nd ech of it hs (n + 1) temrs It s esy to imgine tht for lrge n this is hrd Exmple n = 50 then we hve ingredients There hs to be better wy to do this! From construction: 0 r 1r 2r n 1, r 1r 2r n = 1 We cn tret the r i s probbilities! We define rndom vrible: ξ = (ξ 1,, ξ n ) such tht: δ P(ξ i = 0) = 1 p i, P(ξ i = 1) = p i The extrpoltion vlue is then equl: f(p 1, p 2,, p n) = E [f(ξ 1,, ξ n)] 14 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 14/19 19
15 Trvelling Slesmn Problem Slesmn strting from his bse hs to visit n 1 other loctions nd return to bse hedqurters The problem is to find the shortest wy For lrge n the problem cn t be solver by brutl force s the complexity of the problem is (n 1)! There exist simplified numericl solutions ssuming fctoriztions Unfortuntely even those require nonymous computing power Cn MC help? YES :) The minimum distnce l hs to depend on 2 fctors: P the re of the city the Slesmn is trvelling nd the density of plces he n wnts to visit: P Form this we cn ssume: l P ( n P )b = P b n b 15 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 15/19 19
16 Trveling Slesmn Problem From dimension nlysis: b = 1 2 To get l we need squre root of re From this it s obvious: l P ( n P )b = P 05 n 05 Now we cn multiply the re by lph fctor tht keeps the density constnt then: l α 05 α 05 = α In this cse the distnce between the clients will not chnge, but the number of clients will increse by α so: l α In the end we get: = 1 16 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 16/19 19
17 Trveling Slesmn Problem In totl: l k(np ) 05 Of course the k depends on the shpe of the re nd loctions of client However for lrge n the k strts loosing the dependency It s n symptoticlly free estimtor To use the bove formul we need to somehow clculte k How to estimte this? Well mke TOY MC: tke squre put uniformly n points Then we cn clculte l Then it s trivil: k = l(np ) / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 17/19 19
18 Trveling Slesmn Problem This kind of MC experiment might require lrge CPU power nd time The dventge is tht once we solve the problem we cn use the obtined k for other cses (it s universl constnt!) It turns out tht: k 3 4 Ok, but in this cse we cn clculte l but not the ctul shortest wy! Why the hell we did this exercise?! Turns out tht for most of the problems we re looking for the solution tht is close to smllest l not the exct minimum 18 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 18/19 19
19 Wr Gmes S Andersoon 1966 simulted for Swedish government how would tnk bttle look like Ech of the sides hs 15 tnks tht they llocte on the bttle field The bttle is done in time steps Ech tnk hs 5 sttes: OK Tnk cn only shoot Tnk cn only move Tnk is destroyed Temporry sttes This models mde possible to check different fighting strtegies 19 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 19/19 19
20 Bckup 20 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion 20/19 19
Continuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht
More informationBest Approximation. Chapter The General Case
Chpter 4 Best Approximtion 4.1 The Generl Cse In the previous chpter, we hve seen how n interpolting polynomil cn be used s n pproximtion to given function. We now wnt to find the best pproximtion to given
More informationMonte Carlo method in solving numerical integration and differential equation
Monte Crlo method in solving numericl integrtion nd differentil eqution Ye Jin Chemistry Deprtment Duke University yj66@duke.edu Abstrct: Monte Crlo method is commonly used in rel physics problem. The
More informationOrthogonal Polynomials and Least-Squares Approximations to Functions
Chpter Orthogonl Polynomils nd Lest-Squres Approximtions to Functions **4/5/3 ET. Discrete Lest-Squres Approximtions Given set of dt points (x,y ), (x,y ),..., (x m,y m ), norml nd useful prctice in mny
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationReview of Gaussian Quadrature method
Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge
More informationChapter 3 Polynomials
Dr M DRAIEF As described in the introduction of Chpter 1, pplictions of solving liner equtions rise in number of different settings In prticulr, we will in this chpter focus on the problem of modelling
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More informationPi evaluation. Monte Carlo integration
Pi evlution y 1 1 x Computtionl Physics 2018-19 (Phys Dep IST, Lisbon) Fernndo Bro (311) Monte Crlo integrtion we wnt to evlute the following integrl: F = f (x) dx remember tht the expecttion vlue of the
More informationDefinite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová (Mendel University) Definite integrl MENDELU / Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function
More informationNumerical Integration
Numericl Integrtion Wouter J. Den Hn London School of Economics c 2011 by Wouter J. Den Hn June 3, 2011 Qudrture techniques I = f (x)dx n n w i f (x i ) = w i f i i=1 i=1 Nodes: x i Weights: w i Qudrture
More informationDefinite integral. Mathematics FRDIS MENDELU
Definite integrl Mthemtics FRDIS MENDELU Simon Fišnrová Brno 1 Motivtion - re under curve Suppose, for simplicity, tht y = f(x) is nonnegtive nd continuous function defined on [, b]. Wht is the re of the
More informationThe Wave Equation I. MA 436 Kurt Bryan
1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
More informationa < a+ x < a+2 x < < a+n x = b, n A i n f(x i ) x. i=1 i=1
Mth 33 Volume Stewrt 5.2 Geometry of integrls. In this section, we will lern how to compute volumes using integrls defined by slice nlysis. First, we recll from Clculus I how to compute res. Given the
More informationLecture 17. Integration: Gauss Quadrature. David Semeraro. University of Illinois at Urbana-Champaign. March 20, 2014
Lecture 17 Integrtion: Guss Qudrture Dvid Semerro University of Illinois t Urbn-Chmpign Mrch 0, 014 Dvid Semerro (NCSA) CS 57 Mrch 0, 014 1 / 9 Tody: Objectives identify the most widely used qudrture method
More information7.6 The Use of Definite Integrals in Physics and Engineering
Arknss Tech University MATH 94: Clculus II Dr. Mrcel B. Finn 7.6 The Use of Definite Integrls in Physics nd Engineering It hs been shown how clculus cn be pplied to find solutions to geometric problems
More informationDiscrete Least-squares Approximations
Discrete Lest-squres Approximtions Given set of dt points (x, y ), (x, y ),, (x m, y m ), norml nd useful prctice in mny pplictions in sttistics, engineering nd other pplied sciences is to construct curve
More informationChapter 5 : Continuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 216 Néhémy Lim Chpter 5 : Continuous Rndom Vribles Nottions. N {, 1, 2,...}, set of nturl numbers (i.e. ll nonnegtive integers); N {1, 2,...}, set of ll
More informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution
More informationChapter 5. Numerical Integration
Chpter 5. Numericl Integrtion These re just summries of the lecture notes, nd few detils re included. Most of wht we include here is to be found in more detil in Anton. 5. Remrk. There re two topics with
More informationJDEP 384H: Numerical Methods in Business
BT 3.4: Solving Nonliner Systems Chpter 4: Numericl Integrtion: Deterministic nd Monte Crlo Methods Instructor: Thoms Shores Deprtment of Mthemtics Lecture 20, Februry 29, 2007 110 Kufmnn Center Instructor:
More information1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.
1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt
More informationWe divide the interval [a, b] into subintervals of equal length x = b a n
Arc Length Given curve C defined by function f(x), we wnt to find the length of this curve between nd b. We do this by using process similr to wht we did in defining the Riemnn Sum of definite integrl:
More information20 MATHEMATICS POLYNOMIALS
0 MATHEMATICS POLYNOMIALS.1 Introduction In Clss IX, you hve studied polynomils in one vrible nd their degrees. Recll tht if p(x) is polynomil in x, the highest power of x in p(x) is clled the degree of
More informationAPPLICATIONS OF THE DEFINITE INTEGRAL
APPLICATIONS OF THE DEFINITE INTEGRAL. Volume: Slicing, disks nd wshers.. Volumes by Slicing. Suppose solid object hs boundries extending from x =, to x = b, nd tht its cross-section in plne pssing through
More informationl 2 p2 n 4n 2, the total surface area of the
Week 6 Lectures Sections 7.5, 7.6 Section 7.5: Surfce re of Revolution Surfce re of Cone: Let C be circle of rdius r. Let P n be n n-sided regulr polygon of perimeter p n with vertices on C. Form cone
More informationCS667 Lecture 6: Monte Carlo Integration 02/10/05
CS667 Lecture 6: Monte Crlo Integrtion 02/10/05 Venkt Krishnrj Lecturer: Steve Mrschner 1 Ide The min ide of Monte Crlo Integrtion is tht we cn estimte the vlue of n integrl by looking t lrge number of
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More informationChapter 3 Solving Nonlinear Equations
Chpter 3 Solving Nonliner Equtions 3.1 Introduction The nonliner function of unknown vrible x is in the form of where n could be non-integer. Root is the numericl vlue of x tht stisfies f ( x) 0. Grphiclly,
More informationSUMMER KNOWHOW STUDY AND LEARNING CENTRE
SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18
More informationGoals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite
Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite
More informationJoint distribution. Joint distribution. Marginal distributions. Joint distribution
Joint distribution To specify the joint distribution of n rndom vribles X 1,...,X n tht tke vlues in the smple spces E 1,...,E n we need probbility mesure, P, on E 1... E n = {(x 1,...,x n ) x i E i, i
More informationSturm-Liouville Eigenvalue problem: Let p(x) > 0, q(x) 0, r(x) 0 in I = (a, b). Here we assume b > a. Let X C 2 1
Ch.4. INTEGRAL EQUATIONS AND GREEN S FUNCTIONS Ronld B Guenther nd John W Lee, Prtil Differentil Equtions of Mthemticl Physics nd Integrl Equtions. Hildebrnd, Methods of Applied Mthemtics, second edition
More informationMAA 4212 Improper Integrals
Notes by Dvid Groisser, Copyright c 1995; revised 2002, 2009, 2014 MAA 4212 Improper Integrls The Riemnn integrl, while perfectly well-defined, is too restrictive for mny purposes; there re functions which
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More information1 The Lagrange interpolation formula
Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x
More informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More information7 - Continuous random variables
7-1 Continuous rndom vribles S. Lll, Stnford 2011.01.25.01 7 - Continuous rndom vribles Continuous rndom vribles The cumultive distribution function The uniform rndom vrible Gussin rndom vribles The Gussin
More informationLecture 1. Functional series. Pointwise and uniform convergence.
1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is
More informationRecitation 3: More Applications of the Derivative
Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech
More informationProblem Set 3 Solutions
Chemistry 36 Dr Jen M Stndrd Problem Set 3 Solutions 1 Verify for the prticle in one-dimensionl box by explicit integrtion tht the wvefunction ψ ( x) π x is normlized To verify tht ψ ( x) is normlized,
More informationSummary: Method of Separation of Variables
Physics 246 Electricity nd Mgnetism I, Fll 26, Lecture 22 1 Summry: Method of Seprtion of Vribles 1. Seprtion of Vribles in Crtesin Coordintes 2. Fourier Series Suggested Reding: Griffiths: Chpter 3, Section
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More informationZ b. f(x)dx. Yet in the above two cases we know what f(x) is. Sometimes, engineers want to calculate an area by computing I, but...
Chpter 7 Numericl Methods 7. Introduction In mny cses the integrl f(x)dx cn be found by finding function F (x) such tht F 0 (x) =f(x), nd using f(x)dx = F (b) F () which is known s the nlyticl (exct) solution.
More informationLecture 20: Numerical Integration III
cs4: introduction to numericl nlysis /8/0 Lecture 0: Numericl Integrtion III Instructor: Professor Amos Ron Scribes: Mrk Cowlishw, Yunpeng Li, Nthnel Fillmore For the lst few lectures we hve discussed
More informationINTRODUCTION TO INTEGRATION
INTRODUCTION TO INTEGRATION 5.1 Ares nd Distnces Assume f(x) 0 on the intervl [, b]. Let A be the re under the grph of f(x). b We will obtin n pproximtion of A in the following three steps. STEP 1: Divide
More information4.4 Areas, Integrals and Antiderivatives
. res, integrls nd ntiderivtives 333. Ares, Integrls nd Antiderivtives This section explores properties of functions defined s res nd exmines some connections mong res, integrls nd ntiderivtives. In order
More informationPractice final exam solutions
University of Pennsylvni Deprtment of Mthemtics Mth 26 Honors Clculus II Spring Semester 29 Prof. Grssi, T.A. Asher Auel Prctice finl exm solutions 1. Let F : 2 2 be defined by F (x, y (x + y, x y. If
More informationProblem Set 9. Figure 1: Diagram. This picture is a rough sketch of the 4 parabolas that give us the area that we need to find. The equations are:
(x + y ) = y + (x + y ) = x + Problem Set 9 Discussion: Nov., Nov. 8, Nov. (on probbility nd binomil coefficients) The nme fter the problem is the designted writer of the solution of tht problem. (No one
More informationNUMERICAL INTEGRATION
NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls
More informationUnit #9 : Definite Integral Properties; Fundamental Theorem of Calculus
Unit #9 : Definite Integrl Properties; Fundmentl Theorem of Clculus Gols: Identify properties of definite integrls Define odd nd even functions, nd reltionship to integrl vlues Introduce the Fundmentl
More informationSection 17.2 Line Integrals
Section 7. Line Integrls Integrting Vector Fields nd Functions long urve In this section we consider the problem of integrting functions, both sclr nd vector (vector fields) long curve in the plne. We
More informationRiemann is the Mann! (But Lebesgue may besgue to differ.)
Riemnn is the Mnn! (But Lebesgue my besgue to differ.) Leo Livshits My 2, 2008 1 For finite intervls in R We hve seen in clss tht every continuous function f : [, b] R hs the property tht for every ɛ >
More informationPolynomials and Division Theory
Higher Checklist (Unit ) Higher Checklist (Unit ) Polynomils nd Division Theory Skill Achieved? Know tht polynomil (expression) is of the form: n x + n x n + n x n + + n x + x + 0 where the i R re the
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More informationLine Integrals. Partitioning the Curve. Estimating the Mass
Line Integrls Suppose we hve curve in the xy plne nd ssocite density δ(p ) = δ(x, y) t ech point on the curve. urves, of course, do not hve density or mss, but it my sometimes be convenient or useful to
More informationp(t) dt + i 1 re it ireit dt =
Note: This mteril is contined in Kreyszig, Chpter 13. Complex integrtion We will define integrls of complex functions long curves in C. (This is bit similr to [relvlued] line integrls P dx + Q dy in R2.)
More informationLecture 3 Gaussian Probability Distribution
Introduction Lecture 3 Gussin Probbility Distribution Gussin probbility distribution is perhps the most used distribution in ll of science. lso clled bell shped curve or norml distribution Unlike the binomil
More informationLecture 4: Piecewise Cubic Interpolation
Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 4: Piecewise Cubic Interpoltion Compiled 5 September In this lecture we consider piecewise cubic interpoltion
More informationSection 4.8. D v(t j 1 ) t. (4.8.1) j=1
Difference Equtions to Differentil Equtions Section.8 Distnce, Position, nd the Length of Curves Although we motivted the definition of the definite integrl with the notion of re, there re mny pplictions
More informationChapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1
Tor Kjellsson Stockholm University Chpter 5 5. Strting with the following informtion: R = m r + m r m + m, r = r r we wnt to derive: µ = m m m + m r = R + µ m r, r = R µ m r 3 = µ m R + r, = µ m R r. 4
More informationDOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES
DOIG PHYSICS WITH MATLAB MATHEMATICAL ROUTIES COMPUTATIO OF OE-DIMESIOAL ITEGRALS In Cooper School of Physics, University of Sydney in.cooper@sydney.edu.u DOWLOAD DIRECTORY FOR MATLAB SCRIPTS mth_integrtion_1d.m
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationSolution for Assignment 1 : Intro to Probability and Statistics, PAC learning
Solution for Assignment 1 : Intro to Probbility nd Sttistics, PAC lerning 10-701/15-781: Mchine Lerning (Fll 004) Due: Sept. 30th 004, Thursdy, Strt of clss Question 1. Bsic Probbility ( 18 pts) 1.1 (
More informationHow do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?
XII. LINEAR ALGEBRA: SOLVING SYSTEMS OF EQUATIONS Tody we re going to tlk bout solving systems of liner equtions. These re problems tht give couple of equtions with couple of unknowns, like: 6 2 3 7 4
More informationMath 360: A primitive integral and elementary functions
Mth 360: A primitive integrl nd elementry functions D. DeTurck University of Pennsylvni October 16, 2017 D. DeTurck Mth 360 001 2017C: Integrl/functions 1 / 32 Setup for the integrl prtitions Definition:
More informationLECTURE NOTE #12 PROF. ALAN YUILLE
LECTURE NOTE #12 PROF. ALAN YUILLE 1. Clustering, K-mens, nd EM Tsk: set of unlbeled dt D = {x 1,..., x n } Decompose into clsses w 1,..., w M where M is unknown. Lern clss models p(x w)) Discovery of
More informationB.Sc. in Mathematics (Ordinary)
R48/0 DUBLIN INSTITUTE OF TECHNOLOGY KEVIN STREET, DUBLIN 8 B.Sc. in Mthemtics (Ordinry) SUPPLEMENTAL EXAMINATIONS 01 Numericl Methods Dr. D. Mckey Dr. C. Hills Dr. E.A. Cox Full mrks for complete nswers
More informationChapters 4 & 5 Integrals & Applications
Contents Chpters 4 & 5 Integrls & Applictions Motivtion to Chpters 4 & 5 2 Chpter 4 3 Ares nd Distnces 3. VIDEO - Ares Under Functions............................................ 3.2 VIDEO - Applictions
More informationSection 11.5 Estimation of difference of two proportions
ection.5 Estimtion of difference of two proportions As seen in estimtion of difference of two mens for nonnorml popultion bsed on lrge smple sizes, one cn use CLT in the pproximtion of the distribution
More informationPART 1 MULTIPLE CHOICE Circle the appropriate response to each of the questions below. Each question has a value of 1 point.
PART MULTIPLE CHOICE Circle the pproprite response to ech of the questions below. Ech question hs vlue of point.. If in sequence the second level difference is constnt, thn the sequence is:. rithmetic
More informationNumerical Integration. 1 Introduction. 2 Midpoint Rule, Trapezoid Rule, Simpson Rule. AMSC/CMSC 460/466 T. von Petersdorff 1
AMSC/CMSC 46/466 T. von Petersdorff 1 umericl Integrtion 1 Introduction We wnt to pproximte the integrl I := f xdx where we re given, b nd the function f s subroutine. We evlute f t points x 1,...,x n
More informationImproper Integrals, and Differential Equations
Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted
More informationMath Lecture 23
Mth 8 - Lecture 3 Dyln Zwick Fll 3 In our lst lecture we delt with solutions to the system: x = Ax where A is n n n mtrix with n distinct eigenvlues. As promised, tody we will del with the question of
More information63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1
3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =
More informationEnergy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon
Phys463.nb 49 7 Energy Bnds Ref: textbook, Chpter 7 Q: Why re there insultors nd conductors? Q: Wht will hppen when n electron moves in crystl? In the previous chpter, we discussed free electron gses,
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More informationOrthogonal Polynomials
Mth 4401 Gussin Qudrture Pge 1 Orthogonl Polynomils Orthogonl polynomils rise from series solutions to differentil equtions, lthough they cn be rrived t in vriety of different mnners. Orthogonl polynomils
More informationMORE FUNCTION GRAPHING; OPTIMIZATION. (Last edited October 28, 2013 at 11:09pm.)
MORE FUNCTION GRAPHING; OPTIMIZATION FRI, OCT 25, 203 (Lst edited October 28, 203 t :09pm.) Exercise. Let n be n rbitrry positive integer. Give n exmple of function with exctly n verticl symptotes. Give
More informationLecture 21: Order statistics
Lecture : Order sttistics Suppose we hve N mesurements of sclr, x i =, N Tke ll mesurements nd sort them into scending order x x x 3 x N Define the mesured running integrl S N (x) = 0 for x < x = i/n for
More informationConservation Law. Chapter Goal. 5.2 Theory
Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very
More informationSection 14.3 Arc Length and Curvature
Section 4.3 Arc Length nd Curvture Clculus on Curves in Spce In this section, we ly the foundtions for describing the movement of n object in spce.. Vector Function Bsics In Clc, formul for rc length in
More informationThe Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.
Clculus Li Vs The Fundmentl Theorem of Clculus. The Totl Chnge Theorem nd the Are Under Curve. Recll the following fct from Clculus course. If continuous function f(x) represents the rte of chnge of F
More information13: Diffusion in 2 Energy Groups
3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationAdvanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004
Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when
More informationThe Riemann Integral
Deprtment of Mthemtics King Sud University 2017-2018 Tble of contents 1 Anti-derivtive Function nd Indefinite Integrls 2 3 4 5 Indefinite Integrls & Anti-derivtive Function Definition Let f : I R be function
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationMath 120 Answers for Homework 13
Mth 12 Answers for Homework 13 1. In this problem we will use the fct tht if m f(x M on n intervl [, b] (nd if f is integrble on [, b] then (* m(b f dx M(b. ( The function f(x = 1 + x 3 is n incresing
More information38 Riemann sums and existence of the definite integral.
38 Riemnn sums nd existence of the definite integrl. In the clcultion of the re of the region X bounded by the grph of g(x) = x 2, the x-xis nd 0 x b, two sums ppered: ( n (k 1) 2) b 3 n 3 re(x) ( n These
More informationSession 13
780.20 Session 3 (lst revised: Februry 25, 202) 3 3. 780.20 Session 3. Follow-ups to Session 2 Histogrms of Uniform Rndom Number Distributions. Here is typicl figure you might get when histogrmming uniform
More informationCBE 291b - Computation And Optimization For Engineers
The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn
More information7.2 The Definite Integral
7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where
More information