Solutions to Problem Set 6

Size: px
Start display at page:

Download "Solutions to Problem Set 6"

Transcription

1 Solutons to Problem Set 6 Problem 6. (Resdue theory) a) Problem Boas. n ths problem we wll solve ths ntegral: x sn x x + 4x + 5 dx: To solve ths usng the resdue theorem, we study ths complex ntegral: z z + 4z + 5 ez dz = where s a sem-crcle n the upper half complex plane. When the radus of the sem-crcle goes to nnty, the ntegral of the functon along the arc wll go to zero, and the ntegral along the real axs wll become ths ntegral: z + 4z + 5 dz: f we use Euler's dentty e z = cos z + sn z, then take the magnary part of the ntegral, and substtute z by x, we end up wth our orgnal ntegral. f we had taken the real part, we would have ended up wth the correspondng cosne ntegral. Anyway, n the end, we have to nd the resdues of f (z) at the poles n the upper half plane, and multply the sum of those resdues wth to get the value of the contour ntegral. The magnary part of ths s then the value of our ntegral. Note that ths demands that the degree of the polynomal n the denomnator s greater than (not equal to) the degree of the polynomal n the numerator - f not, the ntegral along the arc would not go to zero. Ths whole recpe s summarzed on page 69 n Boas, where there s a parameter m (whch s one n ths case) that allows for more general ntegrals to be solved. Factorng f (z), we locate the poles of the functon: f (z) = z + 4z + 5 = (z + + )( + ) : There are rst order poles at z =, and the pole at + s n the upper half plane. We need the resdue at ths pole: R(f (z); z = + ) = = lm The contour ntegral s then lm (z + )f (z) z! + z! + z + + = ( + )e ( + )e = = ( )e : Our ntegral s the magnary part of ths: x sn x x + 4x + 5 dx = m :

2 = (cos( )e sn( )e ) = (cos + sn ) : e b) Problem n Boas. n ths problem, we wll solve ths ntegral: cos(x) (4x + 9) dx: Ths ntegral has a very smlar form to the ntegral n the prevous problem, so we go drectly to the recpe on page 69 n Boas. The rst thng we notce s that the ntegrals gven there are between and, but our ntegral s from to. Luckly, the ntegrand s an even functon (f (x) = f ( x)), and they have ths nce property that we have used plenty of tme before: cos(x) (4x + 9) dx = (4x + 9) cos(x)dx: Now we can use the recpe. We check the three ponts:. Snce the numerator s a constant, we can regard t as a polynomal of degree zero. The denomnator s dentely a polynomal, of degree 4.. The numerator s postve for all x R, and so t has no real zeros. 3. The degree of the denomnator s 4, whch s greater than the degree of the numerator (whch s zero). m = n ths case. Wth these three condtons met, t s straght forward to solve the ntegral: just construct the correspondng complex ntegrand f (z) = (4z + 9) ez ; nd the resdues at the poles n the upper half plane, sum them and multply by, and then take the real part (snce we have a cosne n our orgnal ntegrand) to get the nal answer. However, even f t s straght forward, ths may nvolve a lot of work, and a lot of traps to fall nto. So let us get to t. We start by makng an equaton of the denomnator beng zero, to locate the poles and factor the denomnator: (4z + 9) = ) 4z + 9 = ) z = 3 ) f (z) = e z ((z 3 )(z + 3 )) : We see that we have two second order poles, at z = 3. Only the pole at z = 3 s n the upper half plane, so we only need to nd the resdue at ths pole. Usng the method for a second order pole, descrbed near the bottom of page 685 n Boas, we nd the resdue: = R f (z); z = 3 = d dz (z + 3 ) e z (z 3 )ez 6(z + 3 )4 z = 3 z 3 f (z)! z = 3 = (3) 3 6(3) 4 e 3 = e 3 54 :

3 Multplyng ths wth, we get the contour ntegral: = 54 e 3 : Then we take the real part, and multply by the factor /, due to our orgnal ntegraton lmt: cos(x) (4x + 9) dx = Re c) Problem n Boas. Here, we wll study ths ntegral: x sn(x) x dx: = 54 e 3 : Ths ntegral may look suspcous at rst sght - we are ntegratng through two ponts (x = ) n whch the ntegrand s not dened. Now, takng the lmts when x goes towards or, the numerator also goes to zero, and the lmts are dened. By lettng the ntegrand havng these values for x =, there s no problem wth the ntegraton doman. We solve ths by the same recpe as n the two prevous problems, wth one addtonal pont to the resdue theorem: f there are poles on the contour we ntegrate along, we only add half of the resdues there. Gong to the complex functon and factorng the denomnator, we have ths contour ntegral: = ( z)( + z) dz; where s a sem-crcle n the upper half plane wth a radus that wll go to nnty, as usual. The two rst order poles of the ntegrand, at z =, are on the contour. We calculate the resdues: R(f (z); z = ) = lm (z )f (z) = z! ; R(f (z); z = ) = : We add the resdues, dvde by, and multply by (or just smply multply by ): = + = : Now, snce we are nterested n the sne ntegral, we take the magnary part: x sn(x) x dx = m d) Problem n Boas. Now we wll study ths ntegral: x 4 dx: = : Agan there are two ponts n the ntegraton doman (x = ) where the ntegrand s not dened. Ths tme, however, the numerator does not go to 3

4 zero when x goes to, and so we can not smply extend the ntegrand wth addtonal values at these ponts. But we can nd a prncple value. Wth the resdue theorem, ths s straght forward. We rst observe that the ntegrand s an even functon, and so we have x 4 dx = x 4 dx We transform t to a contour ntegral, and factor the denomnator: = (z )(z + )(z )(z + ) dz; where s yet another sem-crcle n the upper half plane wth nnte radus. There s a rst order pole at z = nsde our contour, and two rst order poles, at z =, on the contour. We easly compute the resdues: R(f (z); z = ) = 4 ; R(f (z); z = ) = 4 ; R(f (z); z = ) = 4 : Now, the value of the contour ntegral s tmes the sum of the resdues nsde the contour, plus tmes the sum of the resdues on the contour. The sum of the resdues on the contour s, so we get = 4 = : Now, our orgnal ntegral s related to ths contour ntegral by a factor of /, and we only get a prncpal value: P V x dx = 4 4 : e) Here, we wll study ths ntegral: x + dx: We rst multply wth x, the complex conjugate of the denomnator, n the numerator and denomnator. Then we splt the ntegral n two: (x ) (x + )(x ) dx = x x + dx x + dx: These two ntegrals can be solved straght forward wth the recpe at page 69 n Boas: we solve these contour ntegrals by the resdue theorem ( have already factored the denomnators): f (z)dz = z (z )(z + ) ez dz; f (z)dz = (z )(z + ) ez dz; where s the usual sem-crcle, then we take the real part of both, multply the second by, and then sum up. 4

5 Both ntegrands have poles at z =, where the pole at z = s nsde the contour. The resdues are R(f (z); z = ) = lm(z )f (z) = lm z! z! z + = e R(f (z); z = ) = lm(z )f (z) = lm z! z! We then get these two values for the contour ntegrals: e z z + = f (z)dz = e ; f (z)dz = e : e : Now we take the real parts (snce we are dealng wth the cosne ntegrals): x x + dx = Re x + dx = Re f (z)dz f (z)dz = = e : Then we multply the second term by and add up to get the nal answer: Problem 6. (Tensors) x + dx = e : a) Problem.4.4 n Boas. Fndng the nerta tensor s straght forward by usng equaton (4.5). Snce we are workng wth a contnuous mass dstrbuton, we wll use the ntegral form. The most mportant to note here when comparng to example 5 are the ntegraton lmts. We are studyng a quarter of a sphere. n sphercal coordnates, wth x = r sn cos ; y = r sn sn and z = r cos, we have r [; ]; [; ] and [; ]. The xx component then s = xx = (y + z )dm = (r x )dv V r ( sn cos )r sn drdd = 5 ; where the ntegral s straght forward solved by usng onlne servces such as Wolfram Alpha, or by lookng up n tables. The yy and zz components are found smlarly, and both actually have the same value as xx. Now for the o-dagonal elements, rst xy : xy = xydm = (r 4 sn 4 sn cos )r sn drdd = 5 = yx: The last equalty comes from symmetry n the ntegrand. Smlarly, we nd xz = zx = yz = zy =. As we see, these o-dagonal elements are not equal 5

6 to xy - ths s what the warnng s referrng to. Wrtten as a matrx, the nerta tensor looks lke ths: 5 A : The prncpal moments of nerta are the egenvalues of. We can nd them ether by computer (matlab has very smple functons to compute those), or by hand, by use of the characterstc polynomal. They are 5 ; 5 ; + 5 The prncpal axes are the egenvectors of. Agan, we can nd them by usng a computer, or by hand, by settng up the equatons wth the egenvalues: v = v. You wll be left wth some degrees of freedom here, so have chosen the easest egenvectors by (my) eye. f you use matlab, you wll get normalzed egenvectors. got these three vectors (wrtten n artesan coordnates): (; ; ); (; ; ); (; ; ) : b) Problem.4.6 n Boas. As n the prevous problem, calculatng the nerta tensor n ths case s straght forward usng equaton (4.5). Ths s a dscrete case, and so we use the dscrete versons of the equatons. These are just straght forward sums, where m = for the rst mass, and m = for the second mass. x ; y and z are the gven vector components: xx = X m (y + z ) = 9: Smlarly, we get yy = 9 and zz = 6. For the o-dagonal elements, we get xy = X m x y = 3 = yx ; where the last equalty s yet due to symmetry n the sum. For the other components, we get xz = zx = and yz = zy =. So the nerta tensor looks lke ths: The prncpal moments of nerta are the egenvalues of the nerta tensor: they are f6; 6; g, and the prncpal axes are the egenvectors, where got these: : A : f(; ; ); (; ; ); (; ; )g: 6

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

Math1110 (Spring 2009) Prelim 3 - Solutions

Math1110 (Spring 2009) Prelim 3 - Solutions Math 1110 (Sprng 2009) Solutons to Prelm 3 (04/21/2009) 1 Queston 1. (16 ponts) Short answer. Math1110 (Sprng 2009) Prelm 3 - Solutons x a 1 (a) (4 ponts) Please evaluate lm, where a and b are postve numbers.

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

8.6 The Complex Number System

8.6 The Complex Number System 8.6 The Complex Number System Earler n the chapter, we mentoned that we cannot have a negatve under a square root, snce the square of any postve or negatve number s always postve. In ths secton we want

More information

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order: 68A Solutons to Exercses March 05 (a) Usng a Taylor expanson, and notng that n 0 for all n >, ( + ) ( + ( ) + ) We can t nvert / because there s no Taylor expanson around 0 Lets try to calculate the nverse

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Complex Variables. Chapter 18 Integration in the Complex Plane. March 12, 2013 Lecturer: Shih-Yuan Chen

Complex Variables. Chapter 18 Integration in the Complex Plane. March 12, 2013 Lecturer: Shih-Yuan Chen omplex Varables hapter 8 Integraton n the omplex Plane March, Lecturer: Shh-Yuan hen Except where otherwse noted, content s lcensed under a BY-N-SA. TW Lcense. ontents ontour ntegrals auchy-goursat theorem

More information

Section 3.6 Complex Zeros

Section 3.6 Complex Zeros 04 Chapter Secton 6 Comple Zeros When fndng the zeros of polynomals, at some pont you're faced wth the problem Whle there are clearly no real numbers that are solutons to ths equaton, leavng thngs there

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

Quantum Mechanics I - Session 4

Quantum Mechanics I - Session 4 Quantum Mechancs I - Sesson 4 Aprl 3, 05 Contents Operators Change of Bass 4 3 Egenvectors and Egenvalues 5 3. Denton....................................... 5 3. Rotaton n D....................................

More information

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is. Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these

More information

One Dimension Again. Chapter Fourteen

One Dimension Again. Chapter Fourteen hapter Fourteen One Dmenson Agan 4 Scalar Lne Integrals Now we agan consder the dea of the ntegral n one dmenson When we were ntroduced to the ntegral back n elementary school, we consdered only functons

More information

CIVL 8/7117 Chapter 10 - Isoparametric Formulation 42/56

CIVL 8/7117 Chapter 10 - Isoparametric Formulation 42/56 CIVL 8/77 Chapter 0 - Isoparametrc Formulaton 4/56 Newton-Cotes Example Usng the Newton-Cotes method wth = ntervals (n = 3 samplng ponts), evaluate the ntegrals: x x cos dx 3 x x dx 3 x x dx 4.3333333

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

Digital Signal Processing

Digital Signal Processing Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over

More information

LECTURE 9 CANONICAL CORRELATION ANALYSIS

LECTURE 9 CANONICAL CORRELATION ANALYSIS LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of

More information

18.781: Solution to Practice Questions for Final Exam

18.781: Solution to Practice Questions for Final Exam 18.781: Soluton to Practce Questons for Fnal Exam 1. Fnd three solutons n postve ntegers of x 6y = 1 by frst calculatng the contnued fracton expanson of 6. Soluton: We have 1 6=[, ] 6 6+ =[, ] 1 =[,, ]=[,,

More information

Integrals and Invariants of Euler-Lagrange Equations

Integrals and Invariants of Euler-Lagrange Equations Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,

More information

σ τ τ τ σ τ τ τ σ Review Chapter Four States of Stress Part Three Review Review

σ τ τ τ σ τ τ τ σ Review Chapter Four States of Stress Part Three Review Review Chapter Four States of Stress Part Three When makng your choce n lfe, do not neglect to lve. Samuel Johnson Revew When we use matrx notaton to show the stresses on an element The rows represent the axs

More information

PES 1120 Spring 2014, Spendier Lecture 6/Page 1

PES 1120 Spring 2014, Spendier Lecture 6/Page 1 PES 110 Sprng 014, Spender Lecture 6/Page 1 Lecture today: Chapter 1) Electrc feld due to charge dstrbutons -> charged rod -> charged rng We ntroduced the electrc feld, E. I defned t as an nvsble aura

More information

Statistics and Probability Theory in Civil, Surveying and Environmental Engineering

Statistics and Probability Theory in Civil, Surveying and Environmental Engineering Statstcs and Probablty Theory n Cvl, Surveyng and Envronmental Engneerng Pro. Dr. Mchael Havbro Faber ETH Zurch, Swtzerland Contents o Todays Lecture Overvew o Uncertanty Modelng Random Varables - propertes

More information

Integrals and Invariants of

Integrals and Invariants of Lecture 16 Integrals and Invarants of Euler Lagrange Equatons NPTEL Course Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng, Indan Insttute of Scence, Banagalore

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

COMPLEX NUMBERS AND QUADRATIC EQUATIONS COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

332600_08_1.qxp 4/17/08 11:29 AM Page 481

332600_08_1.qxp 4/17/08 11:29 AM Page 481 336_8_.qxp 4/7/8 :9 AM Page 48 8 Complex Vector Spaces 8. Complex Numbers 8. Conjugates and Dvson of Complex Numbers 8.3 Polar Form and DeMovre s Theorem 8.4 Complex Vector Spaces and Inner Products 8.5

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

Lecture Torsion Properties for Line Segments and Computational Scheme for Piecewise Straight Section Calculations

Lecture Torsion Properties for Line Segments and Computational Scheme for Piecewise Straight Section Calculations Lecture - 003 Torson Propertes for Lne Segments and Computatonal Scheme for Pecewse Straght Secton Calculatons ths conssts of four parts (and how we wll treat each) A - dervaton of geometrc algorthms for

More information

Random Walks on Digraphs

Random Walks on Digraphs Random Walks on Dgraphs J. J. P. Veerman October 23, 27 Introducton Let V = {, n} be a vertex set and S a non-negatve row-stochastc matrx (.e. rows sum to ). V and S defne a dgraph G = G(V, S) and a drected

More information

An Explicit Construction of an Expander Family (Margulis-Gaber-Galil)

An Explicit Construction of an Expander Family (Margulis-Gaber-Galil) An Explct Constructon of an Expander Famly Marguls-Gaber-Gall) Orr Paradse July 4th 08 Prepared for a Theorst's Toolkt, a course taught by Irt Dnur at the Wezmann Insttute of Scence. The purpose of these

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

Unit 5: Quadratic Equations & Functions

Unit 5: Quadratic Equations & Functions Date Perod Unt 5: Quadratc Equatons & Functons DAY TOPIC 1 Modelng Data wth Quadratc Functons Factorng Quadratc Epressons 3 Solvng Quadratc Equatons 4 Comple Numbers Smplfcaton, Addton/Subtracton & Multplcaton

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

12. The Hamilton-Jacobi Equation Michael Fowler

12. The Hamilton-Jacobi Equation Michael Fowler 1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and

More information

General Tips on How to Do Well in Physics Exams. 1. Establish a good habit in keeping track of your steps. For example, when you use the equation

General Tips on How to Do Well in Physics Exams. 1. Establish a good habit in keeping track of your steps. For example, when you use the equation General Tps on How to Do Well n Physcs Exams 1. Establsh a good habt n keepng track o your steps. For example when you use the equaton 1 1 1 + = d d to solve or d o you should rst rewrte t as 1 1 1 = d

More information

Norms, Condition Numbers, Eigenvalues and Eigenvectors

Norms, Condition Numbers, Eigenvalues and Eigenvectors Norms, Condton Numbers, Egenvalues and Egenvectors 1 Norms A norm s a measure of the sze of a matrx or a vector For vectors the common norms are: N a 2 = ( x 2 1/2 the Eucldean Norm (1a b 1 = =1 N x (1b

More information

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

More information

= z 20 z n. (k 20) + 4 z k = 4

= z 20 z n. (k 20) + 4 z k = 4 Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5

More information

10. Canonical Transformations Michael Fowler

10. Canonical Transformations Michael Fowler 10. Canoncal Transformatons Mchael Fowler Pont Transformatons It s clear that Lagrange s equatons are correct for any reasonable choce of parameters labelng the system confguraton. Let s call our frst

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

Using TranSIESTA (II): Integration contour and tbtrans

Using TranSIESTA (II): Integration contour and tbtrans Usng TranSIESTA (II): Integraton contour and tbtrans Frederco D. Novaes December 15, 2009 Outlne Usng TranSIESTA The ntegraton contour Electron densty from GF Why go complex? The non-equlbrum stuaton Usng

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system. Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and

More information

Complex Numbers. x = B B 2 4AC 2A. or x = x = 2 ± 4 4 (1) (5) 2 (1)

Complex Numbers. x = B B 2 4AC 2A. or x = x = 2 ± 4 4 (1) (5) 2 (1) Complex Numbers If you have not yet encountered complex numbers, you wll soon do so n the process of solvng quadratc equatons. The general quadratc equaton Ax + Bx + C 0 has solutons x B + B 4AC A For

More information

MEEM 3700 Mechanical Vibrations

MEEM 3700 Mechanical Vibrations MEEM 700 Mechancal Vbratons Mohan D. Rao Chuck Van Karsen Mechancal Engneerng-Engneerng Mechancs Mchgan echnologcal Unversty Copyrght 00 Lecture & MEEM 700 Multple Degree of Freedom Systems (ext: S.S.

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

(A and B must have the same dmensons to be able to add them together.) Addton s commutatve and assocatve, just lke regular addton. A matrx A multpled

(A and B must have the same dmensons to be able to add them together.) Addton s commutatve and assocatve, just lke regular addton. A matrx A multpled CNS 185: A Bref Revew of Lnear Algebra An understandng of lnear algebra s crtcal as a steppng-o pont for understandng neural networks. Ths handout ncludes basc dentons, then quckly progresses to elementary

More information

Lecture 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

More information

(2mn, m 2 n 2, m 2 + n 2 )

(2mn, m 2 n 2, m 2 + n 2 ) MATH 16T Homewk Solutons 1. Recall that a natural number n N s a perfect square f n = m f some m N. a) Let n = p α even f = 1,,..., k. be the prme factzaton of some n. Prove that n s a perfect square f

More information

SELECTED PROOFS. DeMorgan s formulas: The first one is clear from Venn diagram, or the following truth table:

SELECTED PROOFS. DeMorgan s formulas: The first one is clear from Venn diagram, or the following truth table: SELECTED PROOFS DeMorgan s formulas: The frst one s clear from Venn dagram, or the followng truth table: A B A B A B Ā B Ā B T T T F F F F T F T F F T F F T T F T F F F F F T T T T The second one can be

More information

Linear Feature Engineering 11

Linear Feature Engineering 11 Lnear Feature Engneerng 11 2 Least-Squares 2.1 Smple least-squares Consder the followng dataset. We have a bunch of nputs x and correspondng outputs y. The partcular values n ths dataset are x y 0.23 0.19

More information

Chapter Twelve. Integration. We now turn our attention to the idea of an integral in dimensions higher than one. Consider a real-valued function f : D

Chapter Twelve. Integration. We now turn our attention to the idea of an integral in dimensions higher than one. Consider a real-valued function f : D Chapter Twelve Integraton 12.1 Introducton We now turn our attenton to the dea of an ntegral n dmensons hgher than one. Consder a real-valued functon f : R, where the doman s a nce closed subset of Eucldean

More information

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim Causal Damonds M. Aghl, L. Bombell, B. Plgrm Introducton The correcton to volume of a causal nterval due to curvature of spacetme has been done by Myrhem [] and recently by Gbbons & Solodukhn [] and later

More information

Lagrangian Field Theory

Lagrangian Field Theory Lagrangan Feld Theory Adam Lott PHY 391 Aprl 6, 017 1 Introducton Ths paper s a summary of Chapter of Mandl and Shaw s Quantum Feld Theory [1]. The frst thng to do s to fx the notaton. For the most part,

More information

P A = (P P + P )A = P (I P T (P P ))A = P (A P T (P P )A) Hence if we let E = P T (P P A), We have that

P A = (P P + P )A = P (I P T (P P ))A = P (A P T (P P )A) Hence if we let E = P T (P P A), We have that Backward Error Analyss for House holder Reectors We want to show that multplcaton by householder reectors s backward stable. In partcular we wsh to show fl(p A) = P (A) = P (A + E where P = I 2vv T s the

More information

Electricity and Magnetism - Physics 121 Lecture 10 - Sources of Magnetic Fields (Currents) Y&F Chapter 28, Sec. 1-7

Electricity and Magnetism - Physics 121 Lecture 10 - Sources of Magnetic Fields (Currents) Y&F Chapter 28, Sec. 1-7 Electrcty and Magnetsm - Physcs 11 Lecture 10 - Sources of Magnetc Felds (Currents) Y&F Chapter 8, Sec. 1-7 Magnetc felds are due to currents The Bot-Savart Law Calculatng feld at the centers of current

More information

Work is the change in energy of a system (neglecting heat transfer). To examine what could

Work is the change in energy of a system (neglecting heat transfer). To examine what could Work Work s the change n energy o a system (neglectng heat transer). To eamne what could cause work, let s look at the dmensons o energy: L ML E M L F L so T T dmensonally energy s equal to a orce tmes

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

This model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as:

This model contains two bonds per unit cell (one along the x-direction and the other along y). So we can rewrite the Hamiltonian as: 1 Problem set #1 1.1. A one-band model on a square lattce Fg. 1 Consder a square lattce wth only nearest-neghbor hoppngs (as shown n the fgure above): H t, j a a j (1.1) where,j stands for nearest neghbors

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur Module Random Processes Lesson 6 Functons of Random Varables After readng ths lesson, ou wll learn about cdf of functon of a random varable. Formula for determnng the pdf of a random varable. Let, X be

More information

MTH 263 Practice Test #1 Spring 1999

MTH 263 Practice Test #1 Spring 1999 Pat Ross MTH 6 Practce Test # Sprng 999 Name. Fnd the area of the regon bounded by the graph r =acos (θ). Observe: Ths s a crcle of radus a, for r =acos (θ) r =a ³ x r r =ax x + y =ax x ax + y =0 x ax

More information

From Biot-Savart Law to Divergence of B (1)

From Biot-Savart Law to Divergence of B (1) From Bot-Savart Law to Dvergence of B (1) Let s prove that Bot-Savart gves us B (r ) = 0 for an arbtrary current densty. Frst take the dvergence of both sdes of Bot-Savart. The dervatve s wth respect to

More information

Lecture Note 3. Eshelby s Inclusion II

Lecture Note 3. Eshelby s Inclusion II ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

Complex Numbers, Signals, and Circuits

Complex Numbers, Signals, and Circuits Complex Numbers, Sgnals, and Crcuts 3 August, 009 Complex Numbers: a Revew Suppose we have a complex number z = x jy. To convert to polar form, we need to know the magntude of z and the phase of z. z =

More information

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017) Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed

More information

Programming Project 1: Molecular Geometry and Rotational Constants

Programming Project 1: Molecular Geometry and Rotational Constants Programmng Project 1: Molecular Geometry and Rotatonal Constants Center for Computatonal Chemstry Unversty of Georga Athens, Georga 30602 Summer 2012 1 Introducton Ths programmng project s desgned to provde

More information

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department

More information

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt Physcs 543 Quantum Mechancs II Fall 998 Hartree-Fock and the Self-consstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes.

More information

NEWTON S LAWS. These laws only apply when viewed from an inertial coordinate system (unaccelerated system).

NEWTON S LAWS. These laws only apply when viewed from an inertial coordinate system (unaccelerated system). EWTO S LAWS Consder two partcles. 1 1. If 1 0 then 0 wth p 1 m1v. 1 1 2. 1.. 3. 11 These laws only apply when vewed from an nertal coordnate system (unaccelerated system). consder a collecton of partcles

More information

PHYS 705: Classical Mechanics. Hamilton-Jacobi Equation

PHYS 705: Classical Mechanics. Hamilton-Jacobi Equation 1 PHYS 705: Classcal Mechancs Hamlton-Jacob Equaton Hamlton-Jacob Equaton There s also a very elegant relaton between the Hamltonan Formulaton of Mechancs and Quantum Mechancs. To do that, we need to derve

More information

Formal solvers of the RT equation

Formal solvers of the RT equation Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan

More information

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning

Erratum: A Generalized Path Integral Control Approach to Reinforcement Learning Journal of Machne Learnng Research 00-9 Submtted /0; Publshed 7/ Erratum: A Generalzed Path Integral Control Approach to Renforcement Learnng Evangelos ATheodorou Jonas Buchl Stefan Schaal Department of

More information

The Prncpal Component Transform The Prncpal Component Transform s also called Karhunen-Loeve Transform (KLT, Hotellng Transform, oregenvector Transfor

The Prncpal Component Transform The Prncpal Component Transform s also called Karhunen-Loeve Transform (KLT, Hotellng Transform, oregenvector Transfor Prncpal Component Transform Multvarate Random Sgnals A real tme sgnal x(t can be consdered as a random process and ts samples x m (m =0; ;N, 1 a random vector: The mean vector of X s X =[x0; ;x N,1] T

More information

9. Complex Numbers. 1. Numbers revisited. 2. Imaginary number i: General form of complex numbers. 3. Manipulation of complex numbers

9. Complex Numbers. 1. Numbers revisited. 2. Imaginary number i: General form of complex numbers. 3. Manipulation of complex numbers 9. Comple Numbers. Numbers revsted. Imagnar number : General form of comple numbers 3. Manpulaton of comple numbers 4. The Argand dagram 5. The polar form for comple numbers 9.. Numbers revsted We saw

More information

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16 STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus

More information

Lecture 10: Euler s Equations for Multivariable

Lecture 10: Euler s Equations for Multivariable Lecture 0: Euler s Equatons for Multvarable Problems Let s say we re tryng to mnmze an ntegral of the form: {,,,,,, ; } J f y y y y y y d We can start by wrtng each of the y s as we dd before: y (, ) (

More information

MATH Sensitivity of Eigenvalue Problems

MATH Sensitivity of Eigenvalue Problems MATH 537- Senstvty of Egenvalue Problems Prelmnares Let A be an n n matrx, and let λ be an egenvalue of A, correspondngly there are vectors x and y such that Ax = λx and y H A = λy H Then x s called A

More information

Dynamic Systems on Graphs

Dynamic Systems on Graphs Prepared by F.L. Lews Updated: Saturday, February 06, 200 Dynamc Systems on Graphs Control Graphs and Consensus A network s a set of nodes that collaborates to acheve what each cannot acheve alone. A network,

More information

Effects of Ignoring Correlations When Computing Sample Chi-Square. John W. Fowler February 26, 2012

Effects of Ignoring Correlations When Computing Sample Chi-Square. John W. Fowler February 26, 2012 Effects of Ignorng Correlatons When Computng Sample Ch-Square John W. Fowler February 6, 0 It can happen that ch-square must be computed for a sample whose elements are correlated to an unknown extent.

More information

14 The Postulates of Quantum mechanics

14 The Postulates of Quantum mechanics 14 The Postulates of Quantum mechancs Postulate 1: The state of a system s descrbed completely n terms of a state vector Ψ(r, t), whch s quadratcally ntegrable. Postulate 2: To every physcally observable

More information

1 Derivation of Point-to-Plane Minimization

1 Derivation of Point-to-Plane Minimization 1 Dervaton of Pont-to-Plane Mnmzaton Consder the Chen-Medon (pont-to-plane) framework for ICP. Assume we have a collecton of ponts (p, q ) wth normals n. We want to determne the optmal rotaton and translaton

More information

of Nebraska - Lincoln

of Nebraska - Lincoln Unversty of Nebraska - Lncoln DgtalCommons@Unversty of Nebraska - Lncoln MAT Exam Expostory Papers Math n the Mddle Insttute Partnershp 008 The Square Root of Tffany Lothrop Unversty of Nebraska-Lncoln

More information

Poisson brackets and canonical transformations

Poisson brackets and canonical transformations rof O B Wrght Mechancs Notes osson brackets and canoncal transformatons osson Brackets Consder an arbtrary functon f f ( qp t) df f f f q p q p t But q p p where ( qp ) pq q df f f f p q q p t In order

More information

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0 Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of

More information

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1. 7636S ADVANCED QUANTUM MECHANICS Soluton Set 1 Sprng 013 1 Warm-up Show that the egenvalues of a Hermtan operator  are real and that the egenkets correspondng to dfferent egenvalues are orthogonal (b)

More information

Now that we have laws or better postulates we should explore what they imply

Now that we have laws or better postulates we should explore what they imply I-1 Theorems from Postulates: Now that we have laws or better postulates we should explore what they mply about workng q.m. problems -- Theorems (Levne 7.2, 7.4) Thm 1 -- egen values of Hermtan operators

More information

ALGEBRA MID-TERM. 1 Suppose I is a principal ideal of the integral domain R. Prove that the R-module I R I has no non-zero torsion elements.

ALGEBRA MID-TERM. 1 Suppose I is a principal ideal of the integral domain R. Prove that the R-module I R I has no non-zero torsion elements. ALGEBRA MID-TERM CLAY SHONKWILER 1 Suppose I s a prncpal deal of the ntegral doman R. Prove that the R-module I R I has no non-zero torson elements. Proof. Note, frst, that f I R I has no non-zero torson

More information

DETAILED DERIVATION OF THE GENERAL MASTER EQUATION IN QUANTUM OPTICS

DETAILED DERIVATION OF THE GENERAL MASTER EQUATION IN QUANTUM OPTICS DETAILED DERIVATION OF THE GENERAL MASTER EQUATION IN QUANTUM OPTICS TIMOTHY JONES Abstract. An unorgnal but more detaled dervaton that typcal s presented of the master equaton. 1. Assumptons We begn by

More information

Spectral graph theory: Applications of Courant-Fischer

Spectral graph theory: Applications of Courant-Fischer Spectral graph theory: Applcatons of Courant-Fscher Steve Butler September 2006 Abstract In ths second talk we wll ntroduce the Raylegh quotent and the Courant- Fscher Theorem and gve some applcatons for

More information

Computational Biology Lecture 8: Substitution matrices Saad Mneimneh

Computational Biology Lecture 8: Substitution matrices Saad Mneimneh Computatonal Bology Lecture 8: Substtuton matrces Saad Mnemneh As we have ntroduced last tme, smple scorng schemes lke + or a match, - or a msmatch and -2 or a gap are not justable bologcally, especally

More information

Chapter 3 Differentiation and Integration

Chapter 3 Differentiation and Integration MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton

More information

Spring Force and Power

Spring Force and Power Lecture 13 Chapter 9 Sprng Force and Power Yeah, energy s better than orces. What s net? Course webste: http://aculty.uml.edu/andry_danylov/teachng/physcsi IN THIS CHAPTER, you wll learn how to solve problems

More information

NOTES ON SIMPLIFICATION OF MATRICES

NOTES ON SIMPLIFICATION OF MATRICES NOTES ON SIMPLIFICATION OF MATRICES JONATHAN LUK These notes dscuss how to smplfy an (n n) matrx In partcular, we expand on some of the materal from the textbook (wth some repetton) Part of the exposton

More information

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2 P470 Lecture 6/7 (February 10/1, 014) DIRAC EQUATION The non-relatvstc Schrödnger equaton was obtaned by notng that the Hamltonan H = P (1) m can be transformed nto an operator form wth the substtutons

More information

APPROXIMATE PRICES OF BASKET AND ASIAN OPTIONS DUPONT OLIVIER. Premia 14

APPROXIMATE PRICES OF BASKET AND ASIAN OPTIONS DUPONT OLIVIER. Premia 14 APPROXIMAE PRICES OF BASKE AND ASIAN OPIONS DUPON OLIVIER Prema 14 Contents Introducton 1 1. Framewor 1 1.1. Baset optons 1.. Asan optons. Computng the prce 3. Lower bound 3.1. Closed formula for the prce

More information

UNIT 3 EXPRESSIONS AND EQUATIONS Lesson 4: Fundamental Theorem of Algebra. Instruction. Guided Practice Example 1

UNIT 3 EXPRESSIONS AND EQUATIONS Lesson 4: Fundamental Theorem of Algebra. Instruction. Guided Practice Example 1 Guded Practce 3.4. Example 1 Instructon For each equaton, state the number and type of solutons by frst fndng the dscrmnant. x + 3x =.4x x = 3x = x 9x + 1 = 6x 1. Fnd the dscrmnant of x + 3x =. The equaton

More information