Session 13
|
|
- Arron Bishop
- 5 years ago
- Views:
Transcription
1 Session 3 (lst revised: Februry 25, 202) Session 3. Follow-ups to Session 2 Histogrms of Uniform Rndom Number Distributions. Here is typicl figure you might get when histogrmming uniform distributions using gnuplot (using plot... with steps to get the histogrm look): Two sets of numbers hve n verge of 2000 counts per Uniform Rndom Number Distributions # in bin uniform uniform 2 uniform 3 uniform Thu My 0 :38: bin nd two sets hve n verge of 500 per bin. The size of the fluctutions should scle s # of counts in bin, which is evident in the figure. Tht is, the fluctutions of the top sets look roughly twice s lrge s those of the bottom set nd stndrd devition of looks bout right. The Length of Rndom Wlks. In Session you explored how the finl distnce R from the origin of rndom wlk scles with the number of steps N. Remember tht there is no reltionship for ny individul wlk, but only sttisticl reltionship tht rises by verging R for given N over mny trils. Here is grph tht you might obtin: It shows good fit to the expected N dependence! x
2 Session 3 (lst revised: Februry 25, 202) rndom wlks best-fit line sqrt(# of steps) Rndom Wlks fit slope = finl distnce from origin # of steps Thu My 0 0:22: b. Recp of Session 2: Monte Crlo Integrtion nd Led-in to Session 3 We explored in Session 2 how we could pproximte n integrl f(x) dx by picking rndom sequence of N vlues of x, which we lbel {x i }, distributed uniformly on [, b]. Tht is, the probbility of choosing number between x nd x + dx is P uniform (x) dx = dx. (3.) b Note tht this PDF (probbility distribution function) is normlized, i.e., P uniform (x) dx = b s you d expect since the totl probbility must be unity. dx =, (3.2) The Monte Crlo expression for the integrl rises by pproximting the verge of f(x) in [, b]: We cn write this s f = b f(x) dx f(x i ). (3.3) N i f(x) P uniform (x) dx f(x i ). (3.4) N i
3 Session 3 (lst revised: Februry 25, 202) 3 3 You should hve found tht the ccurcy of the integrl improved sttisticlly with the number of points. Here is grph showing the results of 6 trils for ech vlue of N: Monte Crlo Integrtion monte crlo estimte best-fit line fit slope = bsolute error Thu My 0 0:2: # of integrtion points The expected behvior is reveled by the fit (error improves like / N), but with lrge fluctutions. We should verge over mny more trils if we wnt clener demonstrtion! Now wht if we used different probbility distribution P (x)? Then if our N vlues of x, which we ll lbel { x i }, re distributed ccording to P (x), the verge of the f( x i ) s yields n pproximtion to the integrl over f(x)p (x): N f( x i ) i f(x) P (x) dx. (3.5) This mens tht we cn lwys divide given integrnd into prt tht we tret s probbility distribution (clled P (x) here) nd the rest (which is clled f(x) here). If we choose P (x) so tht f(x) vries only reltively slowly with x, we will get n ccurte estimte of the integrl. This is clled importnce smpling. The generliztion to mny dimensions is immedite: we simply replce one-dimensionl x by multi-dimensionl vector x everywhere. In Session 2 we sw tht GSL cn generte rndom numbers ccording to vrious one-dimensionl PDF s, such s uniform nd gussin distributions. But the multi-dimensionl integrls we ll wnt
4 Session 3 (lst revised: Februry 25, 202) 3 4 to do for physicl problems will generlly hve nonstndrd P (x). We will use the Metropolis lgorithm to generte n pproprite distribution of x i s. c. Thermodynmic Properties of the Ising Model If we wnt to clculte the thermodynmic properties of system t temperture T (tht is, in equilibrium) tht hs Hmiltonin H(x), we hve the problem of clculting high dimensionl integrl with rpidly vrying integrnd. Here H(x) is the energy of the system for given configurtion x. A configurtion is specifiction of the stte of the system in some relevnt vribles. For quntum mechnicl system of N spinless prticles, configurtion could be the N position vectors. For system of spins on lttice, configurtion could be specifiction of the spin projection t ech lttice point. We ll use the Ising model in one nd two dimensions s our first exmple. The system consists of set of lttice sites with spins, which tke on only two possible vlues, + or. The spins re rrnged in chin (one-dimensionl Ising model) or squre lttice (two-dimensionl Ising model). The interction between spins is short-rnged, which is represented in our model by interctions only between nerest neighbors. Tht is, there is potentil energy only from djcent spins. Ech spin hs two nerest neighbors in one dimension nd four nerest neighbors in two dimensions (digonls don t count here!). We hve to specify wht hppens t the edge of the system. E.g., do the ends of the onedimensionl chin hve only one nerest neighbor, or do we imgine the chin wrpped into continuous bnd? In the ltter cse, we hve periodic boundry conditions, becuse it is equivlent to imgining n infinite chin tht repets itself periodiclly. We lso llow for n externl constnt mgnetic field, H, which intercts with ech of the spins individully to contribute Zeemn energy. (Note: the nottion H for the mgnetic field is somewht unfortunte becuse of the confusion with the Hmiltonin. In other references such s Lndu nd Pez, H µb is used insted.) We choose units so tht the Hmiltonin tkes the simple form: H(x) = J i,j S i S j i HS i, (3.6) where S i = ± is the spin. Here x is just shorthnd for specifiction of the spin S i on ech site. For one-dimensionl Ising model with N sites, you could store this s n rry of length N with ech element either + or.
5 Session 3 (lst revised: Februry 25, 202) 3 5 The nottion i, j stnds for nerest neighbor pirs of spins. Obviously the product S i S j cn be either + or. In mny pplictions we ll set the externl field H to zero. The constnt J (clled the exchnge energy ) specifies the strength of the spin-spin interction. For ferromgnetic interction, J > 0, while for n nti-feromgnetic interction J < 0. Wht do you think is the physicl origin of J for ferromgnetic interction? The mgnetiztion of the system is the verge vlue of the spin, S i, where the verge is tken over the entire lttice. If the spins re eqully likely to be up (+) s down ( ), then the net mgnetiztion is zero. In n externl mgnetic field, up or down will be fvored energeticlly (if H > 0 in H(x) bove, which is fvored?), nd there will generlly be net mgnetiztion. In ferromgnet t sufficiently low temperture, there will be spontneous mgnetiztion even in the bsence of n externl field. This my occur over regions ( domins ) or over the entire lttice. A configurtion in the Ising model is microstte of the system. If we hve N sites in onedimensionl Ising model, how mny possible configurtions (or microsttes) re there? How mny different vlues of the energy (with H = 0)? The therml verge of quntity A tht depends on the configurtion x (exmples of A re the energy per degree of freedom E = H T /N or the mgnetiztion per degree of freedom M = i S i T /N) is given by the cnonicl ensemble (the denomintor is the prtition function Z) dx A(x) e H(x)/kT A(x) T = dx = dx A(x) P e H(x )/kt eq (H(x)), (3.7) where we hve pulled out the Boltzmnn fctors to identify the probbility distribution function P eq (H(x)): e H(x)/kT P eq (H(x)) = dx e H(x )/kt. (3.8) (Is this PDF normlized?) In the Ising model exmple, the integrtion over x is just sum over the spin configurtions. If our configurtions re chosen s eigensttes of energy (which they re in the Ising model exmple), then we cn just sum over energies rther thn every configurtion: A T = (# of sttes with energy E) A(E) e E/kT E E (# of sttes with energy E ) e E /kt = E A(E) P (E). (3.9) If we cn construct set of N configurtions { x i } ( sttisticl smple) tht re distributed ccording to P, then we cn pply our Monte Crlo method (i.e., we will be doing importnce smpling) nd A T = (/N) i A( x i). This is wht the Metropolis lgorithm does for us!
6 Session 3 (lst revised: Februry 25, 202) 3 6 d. Metropolis Algorithm The first thing to get stright is tht the Metropolis lgorithm hs nothing to do with Supermn or Fritz Lng :). It is nmed for the first uthor on the pper tht described the lgorithm. (A bit of trivi: nother uthor on the pper is Edwrd Teller.) The Metropolis lgorithm genertes Mrkov chin; this is sequence of configurtions x i tht will be distributed ccording to the cnonicl distribution (i.e., the Boltzmnn fctor will tell us the reltive probbility of ny configurtion). The Mrkov process constructs stte x l+ from the previous stte x l ccording to trnsition probbility W (x l x l+ ). The ide is to construct W such tht in the limit of lrge number of configurtions, the distribution of sttes x i pproches the equilibrium Boltzmn distribution P eq from Eq. (3.8). A sufficient condition for this to hppen is tht the principle of detiled blnce should hold. In equilibrium, the rte of x x should equl the rte of x x, or else it is not equilibrium! In clculting these rtes, there re two terms to multiply: ( ) ( ) probbility probbility of x x (# of x s) = of x x (# of x s) (3.0) time time or which is detiled blnce. This mens tht W (x x ) N P eq (x) = W (x x) N P eq (x ), (3.) W (x x ) W (x x) = P eq(x) P eq (x ). (3.2) But we know tht P eq (x) = Z e H(x)/kT. (3.3) Any choice of W stisfying Eq. (3.2) should work, in principle. For exmple, we cn tke W (x x ) = { e δe/kt if δe > 0 if δe 0 (3.4) with rbitrry for now (set it equl to unity for convenience) nd where δe E x E x. (3.5) We cn check tht it works by considering ll possible outcomes, s shown in this chrt: δe W (x x ) e E/kT δe W (x x) e E/kT E x > E x > 0 e (E x Ex)/kT e Ex/kT < 0 E x < E x < 0 e Ex/kT > 0 e (Ex E x )/kt e E x /kt e E x /kt
7 Session 3 (lst revised: Februry 25, 202) 3 7 If you do the multiplictions, you ll see tht detiled blnce is stisfied in every cse. Section 2.2 of the Binder/Heermn excerpt describes n implementtion of Metropolis for the Ising model. The bsic ide is tht from strting configurtion, one genertes cndidte for new configurtion somehow (for exmple, by flipping one rndomly selected spin). After clculting the energy chnge between new nd old, one hs the trnsition probbility W from Eq. (3.4). Generte rndom number between 0 nd ; if the number is less thn W keep the new configurtion but otherwise flip the spin bck nd keep this s the new configurtion. We ll go over vrious implementtion issues in this Session nd the next. These include: Equilibrtion. If we re going to use configurtions generted by Metropolis to clculte therml verges, we wnt to mke sure we hve reched equilibrium. How long does this tke nd wht is the signture? Efficiency. The initil two-dimensionl Ising model code hs severl inefficiencies tht re generic to mny types of Monte Crlo simultions. In Session 4 we ll look t how to optimize to gin fctor of six in speed. Autocorreltion. When evluting n verge over Monte Crlo configurtions, we wnt to skip the first n steps nd then use dt tken every n 0 Monte Crlo steps. How do we determine how lrge to tke n 0 nd n? The utocorreltion function cn be used to determine resonble choice for n 0. e. References [] R.H. Lndu nd M.J. Pez, Computtionl Physics: Problem Solving with Computers (Wiley- Interscience, 997). [See the info webpge for detils on new version.] [2] M. Hjorth-Jensen, Lecture Notes on Computtionl Physics (2009). These re notes from course offered t the University of Oslo. See the webpge for links to excerpts. [3] W. Press et l., Numericl Recipes in C++, 3rd ed. (Cmbridge, 2007). Chpters from the 2nd edition re vilble online from There re lso Fortrn versions.
6810 Session 13 (last revised: April 3, 2017) 13 1
680 Session 3 (lst revised: April 3, 207) 3 3. 680 Session 3. Follow-ups to Session 2 Histogrms of Uniform Rndom Number Distributions. Here is typicl figure you might get when histogrmming uniform distributions
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 informationTP 10:Importance Sampling-The Metropolis Algorithm-The Ising Model-The Jackknife Method
TP 0:Importnce Smpling-The Metropoli Algorithm-The Iing Model-The Jckknife Method June, 200 The Cnonicl Enemble We conider phyicl ytem which re in therml contct with n environment. The environment i uully
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 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 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 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 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 informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties
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 informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct
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 informationWeek 10: Line Integrals
Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.
More informationn f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1
The Fundmentl Theorem of Clculus As we continue to study the re problem, let s think bck to wht we know bout computing res of regions enclosed by curves. If we wnt to find the re of the region below the
More informationMA 124 January 18, Derivatives are. Integrals are.
MA 124 Jnury 18, 2018 Prof PB s one-minute introduction to clculus Derivtives re. Integrls re. In Clculus 1, we lern limits, derivtives, some pplictions of derivtives, indefinite integrls, definite integrls,
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 information1 Probability Density Functions
Lis Yn CS 9 Continuous Distributions Lecture Notes #9 July 6, 28 Bsed on chpter by Chris Piech So fr, ll rndom vribles we hve seen hve been discrete. In ll the cses we hve seen in CS 9, this ment tht our
More informationContinuous 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 informationAPPROXIMATE INTEGRATION
APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be
More informationReview of basic calculus
Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below
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 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 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 informationMassachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 6
Msschusetts Institute of Technology Quntum Mechnics I (8.) Spring 5 Solutions to Problem Set 6 By Kit Mtn. Prctice with delt functions ( points) The Dirc delt function my be defined s such tht () (b) 3
More informationIntegral equations, eigenvalue, function interpolation
Integrl equtions, eigenvlue, function interpoltion Mrcin Chrząszcz mchrzsz@cernch Monte Crlo methods, 26 My, 2016 1 / Mrcin Chrząszcz (Universität Zürich) Integrl equtions, eigenvlue, function interpoltion
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 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 informationand that at t = 0 the object is at position 5. Find the position of the object at t = 2.
7.2 The Fundmentl Theorem of Clculus 49 re mny, mny problems tht pper much different on the surfce but tht turn out to be the sme s these problems, in the sense tht when we try to pproimte solutions we
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 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 informationMATH , Calculus 2, Fall 2018
MATH 36-2, 36-3 Clculus 2, Fll 28 The FUNdmentl Theorem of Clculus Sections 5.4 nd 5.5 This worksheet focuses on the most importnt theorem in clculus. In fct, the Fundmentl Theorem of Clculus (FTC is rgubly
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 informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationOverview of Calculus I
Overview of Clculus I Prof. Jim Swift Northern Arizon University There re three key concepts in clculus: The limit, the derivtive, nd the integrl. You need to understnd the definitions of these three things,
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 informationChapter 0. What is the Lebesgue integral about?
Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous
More informationP 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)
1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this
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 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 information5.7 Improper Integrals
458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the
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 informationf(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral
Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one
More informationACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019
ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS MATH00030 SEMESTER 208/209 DR. ANTHONY BROWN 7.. Introduction to Integrtion. 7. Integrl Clculus As ws the cse with the chpter on differentil
More informationHow can we approximate the area of a region in the plane? What is an interpretation of the area under the graph of a velocity function?
Mth 125 Summry Here re some thoughts I ws hving while considering wht to put on the first midterm. The core of your studying should be the ssigned homework problems: mke sure you relly understnd those
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 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 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 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 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 informationThe First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).
The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples
More informationCHM Physical Chemistry I Chapter 1 - Supplementary Material
CHM 3410 - Physicl Chemistry I Chpter 1 - Supplementry Mteril For review of some bsic concepts in mth, see Atkins "Mthemticl Bckground 1 (pp 59-6), nd "Mthemticl Bckground " (pp 109-111). 1. Derivtion
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 information221B Lecture Notes WKB Method
Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using
More informationARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac
REVIEW OF ALGEBRA Here we review the bsic rules nd procedures of lgebr tht you need to know in order to be successful in clculus. ARITHMETIC OPERATIONS The rel numbers hve the following properties: b b
More informationMath 426: Probability Final Exam Practice
Mth 46: Probbility Finl Exm Prctice. Computtionl problems 4. Let T k (n) denote the number of prtitions of the set {,..., n} into k nonempty subsets, where k n. Argue tht T k (n) kt k (n ) + T k (n ) by
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 information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
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 informationConsequently, the temperature must be the same at each point in the cross section at x. Let:
HW 2 Comments: L1-3. Derive the het eqution for n inhomogeneous rod where the therml coefficients used in the derivtion of the het eqution for homogeneous rod now become functions of position x in the
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 informationOperations with Polynomials
38 Chpter P Prerequisites P.4 Opertions with Polynomils Wht you should lern: How to identify the leding coefficients nd degrees of polynomils How to dd nd subtrct polynomils How to multiply polynomils
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 informationWe partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.
Mth 255 - Vector lculus II Notes 4.2 Pth nd Line Integrls We begin with discussion of pth integrls (the book clls them sclr line integrls). We will do this for function of two vribles, but these ides cn
More informationClassical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011
Clssicl Mechnics From Moleculr to Con/nuum Physics I WS 11/12 Emilino Ippoli/ October, 2011 Wednesdy, October 12, 2011 Review Mthemtics... Physics Bsic thermodynmics Temperture, idel gs, kinetic gs theory,
More informationReinforcement learning II
CS 1675 Introduction to Mchine Lerning Lecture 26 Reinforcement lerning II Milos Huskrecht milos@cs.pitt.edu 5329 Sennott Squre Reinforcement lerning Bsics: Input x Lerner Output Reinforcement r Critic
More informationName Solutions to Test 3 November 8, 2017
Nme Solutions to Test 3 November 8, 07 This test consists of three prts. Plese note tht in prts II nd III, you cn skip one question of those offered. Some possibly useful formuls cn be found below. Brrier
More informationAbstract inner product spaces
WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the
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 informationPhysics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2011
Physics 74 Grdute Quntum Mechnics Solutions to Finl Exm, Fll 0 You my use () clss notes, () former homeworks nd solutions (vilble online), (3) online routines, such s Clebsch, provided by me, or (4) ny
More information13.4. Integration by Parts. Introduction. Prerequisites. Learning Outcomes
Integrtion by Prts 13.4 Introduction Integrtion by Prts is technique for integrting products of functions. In this Section you will lern to recognise when it is pproprite to use the technique nd hve the
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 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 informationDiscrete Mathematics and Probability Theory Spring 2013 Anant Sahai Lecture 17
EECS 70 Discrete Mthemtics nd Proility Theory Spring 2013 Annt Shi Lecture 17 I.I.D. Rndom Vriles Estimting the is of coin Question: We wnt to estimte the proportion p of Democrts in the US popultion,
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 informationPhysics 202H - Introductory Quantum Physics I Homework #08 - Solutions Fall 2004 Due 5:01 PM, Monday 2004/11/15
Physics H - Introductory Quntum Physics I Homework #8 - Solutions Fll 4 Due 5:1 PM, Mondy 4/11/15 [55 points totl] Journl questions. Briefly shre your thoughts on the following questions: Of the mteril
More informationNumerical Analysis. 10th ed. R L Burden, J D Faires, and A M Burden
Numericl Anlysis 10th ed R L Burden, J D Fires, nd A M Burden Bemer Presenttion Slides Prepred by Dr. Annette M. Burden Youngstown Stte University July 9, 2015 Chpter 4.1: Numericl Differentition 1 Three-Point
More informationLecture 1: Introduction to integration theory and bounded variation
Lecture 1: Introduction to integrtion theory nd bounded vrition Wht is this course bout? Integrtion theory. The first question you might hve is why there is nything you need to lern bout integrtion. You
More informationMain topics for the First Midterm
Min topics for the First Midterm The Midterm will cover Section 1.8, Chpters 2-3, Sections 4.1-4.8, nd Sections 5.1-5.3 (essentilly ll of the mteril covered in clss). Be sure to know the results of the
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 informationMath Calculus with Analytic Geometry II
orem of definite Mth 5.0 with Anlytic Geometry II Jnury 4, 0 orem of definite If < b then b f (x) dx = ( under f bove x-xis) ( bove f under x-xis) Exmple 8 0 3 9 x dx = π 3 4 = 9π 4 orem of definite Problem
More informationMath 113 Exam 2 Practice
Mth Em Prctice Februry, 8 Em will cover sections 6.5, 7.-7.5 nd 7.8. This sheet hs three sections. The first section will remind you bout techniques nd formuls tht you should know. The second gives number
More informationx = b a n x 2 e x dx. cdx = c(b a), where c is any constant. a b
CHAPTER 5. INTEGRALS 61 where nd x = b n x i = 1 (x i 1 + x i ) = midpoint of [x i 1, x i ]. Problem 168 (Exercise 1, pge 377). Use the Midpoint Rule with the n = 4 to pproximte 5 1 x e x dx. Some quick
More informationDefinition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim
Mth 9 Course Summry/Study Guide Fll, 2005 [1] Limits Definition of Limit: We sy tht L is the limit of f(x) s x pproches if f(x) gets closer nd closer to L s x gets closer nd closer to. We write lim f(x)
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 informationModule 6: LINEAR TRANSFORMATIONS
Module 6: LINEAR TRANSFORMATIONS. Trnsformtions nd mtrices Trnsformtions re generliztions of functions. A vector x in some set S n is mpped into m nother vector y T( x). A trnsformtion is liner if, for
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 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 informationLesson 1: Quadratic Equations
Lesson 1: Qudrtic Equtions Qudrtic Eqution: The qudrtic eqution in form is. In this section, we will review 4 methods of qudrtic equtions, nd when it is most to use ech method. 1. 3.. 4. Method 1: Fctoring
More informationLecture 3. Limits of Functions and Continuity
Lecture 3 Limits of Functions nd Continuity Audrey Terrs April 26, 21 1 Limits of Functions Notes I m skipping the lst section of Chpter 6 of Lng; the section bout open nd closed sets We cn probbly live
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 informationMath 113 Fall Final Exam Review. 2. Applications of Integration Chapter 6 including sections and section 6.8
Mth 3 Fll 0 The scope of the finl exm will include: Finl Exm Review. Integrls Chpter 5 including sections 5. 5.7, 5.0. Applictions of Integrtion Chpter 6 including sections 6. 6.5 nd section 6.8 3. Infinite
More information221A Lecture Notes WKB Method
A Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using ψ x, t = e
More informationQuantum Physics I (8.04) Spring 2016 Assignment 8
Quntum Physics I (8.04) Spring 206 Assignment 8 MIT Physics Deprtment Due Fridy, April 22, 206 April 3, 206 2:00 noon Problem Set 8 Reding: Griffiths, pges 73-76, 8-82 (on scttering sttes). Ohnin, Chpter
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 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 informationWe will see what is meant by standard form very shortly
THEOREM: For fesible liner progrm in its stndrd form, the optimum vlue of the objective over its nonempty fesible region is () either unbounded or (b) is chievble t lest t one extreme point of the fesible
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 informationSection 6.1 Definite Integral
Section 6.1 Definite Integrl Suppose we wnt to find the re of region tht is not so nicely shped. For exmple, consider the function shown elow. The re elow the curve nd ove the x xis cnnot e determined
More information