Numerical Methods 1. February 20, Joel Franklin, Cambridge University Press, 2013 and is used with the author s permission.
|
|
- Gerard Potter
- 6 years ago
- Views:
Transcription
1 Numerical Methods 1 February 20, Some of this material appears in Computational Methods for Physics by Joel Franklin, Cambridge University Press, 2013 and is used with the author s permission.
2 1.1 Root-Finding Given a function F (x), find some or all of the set { x i } n i=1 such that F ( x i)= 0. This defines the root-finding problem this type of uestion shows up in a variety of settings E&M example For the time-dependent fields in electricity and magnetism, we know that information (in that setting, the magnitude and direction of E and B at a particular point r = x ˆx + y ŷ + z ẑ at time t) travels at a finite and specific speed: c (in vacuum). The immediate implication is that if we would like to know the electric field at our observation point (and time) due to a charge that is moving along some prescribed trajectory: w(t), we need to evaluate the location of the particle not at time t, but at some earlier time t r (the retarded time). The setup is shown in Figure 1.1 the particle of charge is at w(t) at time t, but it is the earlier location w(t r ) that has information traveling at c to the observation point r at time t. The time it takes for the information from w(t r ) to reach r is t t r, and during that interval, the field information travels at c, so the distance travelled is c (t t r ). At the same time, the distance between r and w(t r ) is, geometrically, r w(t r ) puttingthese together, we have the defining euation for the retarded time: c (t t r )= r w(t r ) p (r w(t r )) (r w(t r )). (1.1) Notice that the retarded time t r is a function of t and r (the observation time and location), and is defined implicitly in terms of w(t r ). We may or may not be able to solve this euation for t r given w analytically (cannot be done for all but the simplest trajectories), but we can define a root-finding problem: F (T )=c (t T ) r w(t ). (1.2) The roots of F (T ) will provide the retarded time Method The procedure we will use is called bisection start with a pair of values, x 0` and x0 r with x 0` <x0 r (hence the ` and r designations) and F (x 0`) F (x0 r) < 0 1
3 ẑ w(t) w(t r ) ˆx r at t ŷ Figure 1.1: A particle with charge moves along the prescribed curve w(t). For the point r at time t, field information comes from the earlier time t r, when the charge was at w(t r ). (so that a root is in between the points). Then calculate the midpoint between this pair, x 0 m 1 2 x0` + x0 r, and evaluate the product p F (x 0`) F (x0 m). If p is less than zero, then the root lies between the left and middle points, so we can move the left and right points over by setting x 1` = x0`, x1 r = x 0 m.if p is greater than zero, the root is between the middle and right points, and we update x 1` = x0 m, x 1 r = x 0 r. The iteration is continued until the absolute value of F at the current midpoint is smaller than some tolerance, F (x n m) apple, (1.3) and then x n m is an approximation to the root. The process is shown pictorially in Figure 1.2, and is described with pseudocode in Algorithm Example the electric potential of a moving charge Given a vector describing the location of a charge at time t, w(t), we can use root-finding to find the retarded time associated with a field point r at time t. From that, we can evaluate the full electric field, or, because it is a 2
4 F (x) x x 0 x 0 m x 0 r x 1 x 1 m x 1 r Figure 1.2: Two iterations of the bisection procedure. Algorithm 1.1 Bisection(F, x 0`,x0 r, ) x` x r x 0` x 0 r x m 1/2(x` + x r ) while F (x m ) > do if F (x m ) F (x r ) < 0 then x` else x r x m x m end if x m 1/2(x` + x r ) end while return x m 3
5 scalar, we can focus on the potential: V (r,t)= r w(t r ) (r w(t r )) ẇ(t r )/c, (1.4) which reduces, for slow source-charge motion, to V (r,t)= 4 0 r w(t r ) (1.5) (we ll use this approximation). Given a charge s trajectory through space, w(t), we can use the routine described in Algorithm 1.2 to generate the (approximate) electric potential at an individual field point, and by evaluating the potential at a variety of points (in a grid, for example), we can generate a contour plot of it (that process is described in Problem 1.3 below). Algorithm 1.2 Vfield(x, y, z, t, w, ) r {x, y, z} F (X) c (t X) r w(x) tr p Bisection(F, 1,t, ) R (r w(tr )) (r w(t r )) 4 0 R V return V Problem 1.0 When you make a function in Mathematica (using Module), you cannot assign values to an argument of the function. To see what happens (and gain familiarity with the arcane error messages Mathematica has in store for you), write a function that adds two numbers and assigns the result to one of the inputs as follows: addem[a, b ] := Module[{}, a = a + b; Return[a];]. Try it out: addem[5,3] whathappens? Another potential pitfall in writing functions can occur when you try to return a variable that has no assigned value. Let s rewrite the addem function so that the return value is not assigned: addem[a, b ] := Module[{retval, c}, c = a + b; Return[retval];]. What output does addem[5,3] produce this time? 4
6 Problem 1.1 Implement the bisection procedure, and test it by finding the roots of f(x) = (x )(x+1) and the first two roots of J 0 (x) (obtainable via BesselJ[0,x] ) what are those roots? Use = as the tolerance in all cases. Problem 1.2 A charged particle moves along the trajectory described by w(t) =R cos(!t) ˆx + d!tẑ. (1.6) Use the bisection routine to find the retarded time for the point r =0at t =0 if R = 10 5 m,! = /s, d =1m. Use = as your tolerance. What is w(t r )? How does this compare with w(0), the value we would use if we ignored the retarded time? Problem 1.3 Given w(t) =d cos(!t) ŷ with d =.05 m,! =5s 1,andtaking/(4 0 )=1, c =1for simplicity, implement Algorithm 1.2 to find V and make contour plots of its value for points between 1! 1 in x and.1! 1 in y, withz =0and t going from 0 to 8 /5 in steps of 2 /(25) using: data=table[contourplot[vfield[x,y,0,t,w,eps],{x,-1,1},{y,.1,1}], {t,0,8 Pi/5, (2 Pi/5)/5}] You can animate these as a flip book via ListAnimate[data]. Use = again, and take 1 = 10 for the bisection inside Algorithm
Physics 202 Laboratory 3. Root-Finding 1. Laboratory 3. Physics 202 Laboratory
Physics 202 Laboratory 3 Root-Finding 1 Laboratory 3 Physics 202 Laboratory The fundamental question answered by this week s lab work will be: Given a function F (x), find some/all of the values {x i }
More informationPhysics 200 Lecture 4. Integration. Lecture 4. Physics 200 Laboratory
Physics 2 Lecture 4 Integration Lecture 4 Physics 2 Laboratory Monday, February 21st, 211 Integration is the flip-side of differentiation in fact, it is often possible to write a differential equation
More information2 Discrete Dynamical Systems (DDS)
2 Discrete Dynamical Systems (DDS) 2.1 Basics A Discrete Dynamical System (DDS) models the change (or dynamics) of single or multiple populations or quantities in which the change occurs deterministically
More informationGoals for This Lecture:
Goals for This Lecture: Learn the Newton-Raphson method for finding real roots of real functions Learn the Bisection method for finding real roots of a real function Look at efficient implementations of
More informationLinear models: the perceptron and closest centroid algorithms. D = {(x i,y i )} n i=1. x i 2 R d 9/3/13. Preliminaries. Chapter 1, 7.
Preliminaries Linear models: the perceptron and closest centroid algorithms Chapter 1, 7 Definition: The Euclidean dot product beteen to vectors is the expression d T x = i x i The dot product is also
More informationPART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435
PART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435 Professor Biswa Nath Datta Department of Mathematical Sciences Northern Illinois University DeKalb, IL. 60115 USA E mail: dattab@math.niu.edu
More informationgap> G:=ReducedGroebnerBasis(B,MonomialGrlexOrdering(y,x)); [ x^2-x, x*y-x, y^2-x ]
Mathematics 667 Gröbner bases in GAP A. Hulpke To start using it, one has to define indeterminates (variables) and polynomials. Indeterminates are displayed either by their internal numbers or you can
More informationSubstitutions and by Parts, Area Between Curves. Goals: The Method of Substitution Areas Integration by Parts
Week #7: Substitutions and by Parts, Area Between Curves Goals: The Method of Substitution Areas Integration by Parts 1 Week 7 The Indefinite Integral The Fundamental Theorem of Calculus, b a f(x) dx =
More informationPhysics 141 Energy 1 Page 1. Energy 1
Physics 4 Energy Page Energy What I tell you three times is true. Lewis Carroll The interplay of mathematics and physics The mathematization of physics in ancient times is attributed to the Pythagoreans,
More informationModeling and Predicting Chaotic Time Series
Chapter 14 Modeling and Predicting Chaotic Time Series To understand the behavior of a dynamical system in terms of some meaningful parameters we seek the appropriate mathematical model that captures the
More information1.1: The bisection method. September 2017
(1/11) 1.1: The bisection method Solving nonlinear equations MA385/530 Numerical Analysis September 2017 3 2 f(x)= x 2 2 x axis 1 0 1 x [0] =a x [2] =1 x [3] =1.5 x [1] =b 2 0.5 0 0.5 1 1.5 2 2.5 1 Solving
More informationNumerical Methods 5633
Numerical Methods 5633 Lecture 4 Michaelmas Term 2017 Marina Krstic Marinkovic mmarina@maths.tcd.ie School of Mathematics Trinity College Dublin Marina Krstic Marinkovic 1 / 17 5633-Numerical Methods Root
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationSolution of Nonlinear Equations
Solution of Nonlinear Equations (Com S 477/577 Notes) Yan-Bin Jia Sep 14, 017 One of the most frequently occurring problems in scientific work is to find the roots of equations of the form f(x) = 0. (1)
More information1 The best of all possible worlds
Notes for 2017-03-18 1 The best of all possible worlds Last time, we discussed three methods of solving f(x) = 0: Newton, modified Newton, and bisection. Newton is potentially faster than bisection; bisection
More informationRoot Finding: Close Methods. Bisection and False Position Dr. Marco A. Arocha Aug, 2014
Root Finding: Close Methods Bisection and False Position Dr. Marco A. Arocha Aug, 2014 1 Roots Given function f(x), we seek x values for which f(x)=0 Solution x is the root of the equation or zero of the
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationGraphics (INFOGR), , Block IV, lecture 8 Deb Panja. Today: Matrices. Welcome
Graphics (INFOGR), 2017-18, Block IV, lecture 8 Deb Panja Today: Matrices Welcome 1 Today Matrices: why and what? Matrix operations Determinants Adjoint/adjugate and inverse of matrices Geometric interpretation
More informationBisected Direct Quadratic Regula Falsi
Applied Mathematical Sciences, Vol. 4, 2010, no. 15, 709-718 Bisected Direct Quadratic Regula Falsi Robert G. Gottlieb 1 and Blair F. Thompson Odyssey Space Research 1120 NASA Parkway, Suite 505 Houston,
More informationREVIEW: Waves on a String
Lecture 14: Solution to the Wave Equation (Chapter 6) and Random Walks (Chapter 7) 1 Description of Wave Motion REVIEW: Waves on a String We are all familiar with the motion of a transverse wave pulse
More informationPOLARIZATION OF LIGHT
POLARIZATION OF LIGHT OVERALL GOALS The Polarization of Light lab strongly emphasizes connecting mathematical formalism with measurable results. It is not your job to understand every aspect of the theory,
More informationNumerical Methods Dr. Sanjeev Kumar Department of Mathematics Indian Institute of Technology Roorkee Lecture No 7 Regula Falsi and Secant Methods
Numerical Methods Dr. Sanjeev Kumar Department of Mathematics Indian Institute of Technology Roorkee Lecture No 7 Regula Falsi and Secant Methods So welcome to the next lecture of the 2 nd unit of this
More informationAnalysis Methods in Atmospheric and Oceanic Science
Analysis Methods in Atmospheric and Oceanic Science 1 AOSC 652 Week 7, Day 1 10 Oct 2016 Student projects: 20% of the final grade: you will receive a numerical score for the project and final grade will
More informationHow to Use Calculus Like a Physicist
How to Use Calculus Like a Physicist Physics A300 Fall 2004 The purpose of these notes is to make contact between the abstract descriptions you may have seen in your calculus classes and the applications
More informationIntroduction to Machine Learning Spring 2018 Note Neural Networks
CS 189 Introduction to Machine Learning Spring 2018 Note 14 1 Neural Networks Neural networks are a class of compositional function approximators. They come in a variety of shapes and sizes. In this class,
More informationAn Improved Hybrid Algorithm to Bisection Method and Newton-Raphson Method
Applied Mathematical Sciences, Vol. 11, 2017, no. 56, 2789-2797 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.710302 An Improved Hybrid Algorithm to Bisection Method and Newton-Raphson
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 5. Nonlinear Equations
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T Heath Chapter 5 Nonlinear Equations Copyright c 2001 Reproduction permitted only for noncommercial, educational
More informationM01M.1 Massive Spring
Part I Mechanics M01M.1 Massive Spring M01M.1 Massive Spring A spring has spring constant K, unstretched length L, and mass per unit length ρ. The spring is suspended vertically from one end in a constant
More informationCS 125 Section #12 (More) Probability and Randomized Algorithms 11/24/14. For random numbers X which only take on nonnegative integer values, E(X) =
CS 125 Section #12 (More) Probability and Randomized Algorithms 11/24/14 1 Probability First, recall a couple useful facts from last time about probability: Linearity of expectation: E(aX + by ) = ae(x)
More informationIterative methods: introduction and background. Vrije Universiteit Amsterdam
Iterative methods: introduction and background André Ran Vrije Universiteit Amsterdam 1 Consider the set of linear equations (I + αt)u = v, where v is given, α > 0 and T is the symmetric tridiagonal Toeplitz
More informationMath Numerical Analysis
Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center
More informationSolving Non-Linear Equations (Root Finding)
Solving Non-Linear Equations (Root Finding) Root finding Methods What are root finding methods? Methods for determining a solution of an equation. Essentially finding a root of a function, that is, a zero
More informationOutline. Math Numerical Analysis. Intermediate Value Theorem. Lecture Notes Zeros and Roots. Joseph M. Mahaffy,
Outline Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research
More informationIntroduction to Finite Di erence Methods
Introduction to Finite Di erence Methods ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu ME 448/548: Introduction to Finite Di erence Approximation
More informationWaves, the Wave Equation, and Phase Velocity
Waves, the Wave Equation, and Phase Velocity What is a wave? The one-dimensional wave equation Wavelength, frequency, period, etc. Phase velocity Complex numbers and exponentials Plane waves, laser beams,
More informationChapter 1. Root Finding Methods. 1.1 Bisection method
Chapter 1 Root Finding Methods We begin by considering numerical solutions to the problem f(x) = 0 (1.1) Although the problem above is simple to state it is not always easy to solve analytically. This
More informationThe Inside-Outside Algorithm
The Inside-Outside Algorithm Karl Stratos 1 The Setup Our data consists of Q observations O 1,..., O Q where each observation is a sequence of symbols in some set X = {1,..., n}. We suspect a Probabilistic
More information1-D Optimization. Lab 16. Overview of Line Search Algorithms. Derivative versus Derivative-Free Methods
Lab 16 1-D Optimization Lab Objective: Many high-dimensional optimization algorithms rely on onedimensional optimization methods. In this lab, we implement four line search algorithms for optimizing scalar-valued
More informationNonlinear Equations. Chapter The Bisection Method
Chapter 6 Nonlinear Equations Given a nonlinear function f(), a value r such that f(r) = 0, is called a root or a zero of f() For eample, for f() = e 016064, Fig?? gives the set of points satisfying y
More informationMaths III - Numerical Methods
Maths III - Numerical Methods Matt Probert matt.probert@york.ac.uk 4 Solving Differential Equations 4.1 Introduction Many physical problems can be expressed as differential s, but solving them is not always
More informationcha1873x_p02.qxd 3/21/05 1:01 PM Page 104 PART TWO
cha1873x_p02.qxd 3/21/05 1:01 PM Page 104 PART TWO ROOTS OF EQUATIONS PT2.1 MOTIVATION Years ago, you learned to use the quadratic formula x = b ± b 2 4ac 2a to solve f(x) = ax 2 + bx + c = 0 (PT2.1) (PT2.2)
More informationElements of Vector Calculus : Scalar Field & its Gradient
Elements of Vector Calculus : Scalar Field & its Gradient Lecture 1 : Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Introduction : In this set of approximately 40 lectures
More informationLecture 2: Linear Algebra Review
EE 227A: Convex Optimization and Applications January 19 Lecture 2: Linear Algebra Review Lecturer: Mert Pilanci Reading assignment: Appendix C of BV. Sections 2-6 of the web textbook 1 2.1 Vectors 2.1.1
More informationLinear Regression. S. Sumitra
Linear Regression S Sumitra Notations: x i : ith data point; x T : transpose of x; x ij : ith data point s jth attribute Let {(x 1, y 1 ), (x, y )(x N, y N )} be the given data, x i D and y i Y Here D
More informationPhysics 20 Homework 2 SIMS 2016
Physics 20 Homework 2 SIMS 2016 Due: Saturday, August 20 th 1. In class, we ignored air resistance in our discussion of projectile motion. Now, let s derive the relevant equation of motion in the case
More informationChapter 3: Root Finding. September 26, 2005
Chapter 3: Root Finding September 26, 2005 Outline 1 Root Finding 2 3.1 The Bisection Method 3 3.2 Newton s Method: Derivation and Examples 4 3.3 How To Stop Newton s Method 5 3.4 Application: Division
More informationthe EL equation for the x coordinate is easily seen to be (exercise)
Physics 6010, Fall 2016 Relevant Sections in Text: 1.3 1.6 Examples After all this formalism it is a good idea to spend some time developing a number of illustrative examples. These examples represent
More informationProjectile Motion and 2-D Dynamics
Projectile Motion and 2-D Dynamics Vector Notation Vectors vs. Scalars In Physics 11, you learned the difference between vectors and scalars. A vector is a quantity that includes both direction and magnitude
More informationCLASS NOTES Models, Algorithms and Data: Introduction to computing 2018
CLASS NOTES Models, Algorithms and Data: Introduction to computing 2018 Petros Koumoutsakos, Jens Honore Walther (Last update: April 16, 2018) IMPORTANT DISCLAIMERS 1. REFERENCES: Much of the material
More informationLecture 21 Reminder/Introduction to Wave Optics
Lecture 1 Reminder/Introduction to Wave Optics Program: 1. Maxwell s Equations.. Magnetic induction and electric displacement. 3. Origins of the electric permittivity and magnetic permeability. 4. Wave
More informationPhysics 342 Lecture 2. Linear Algebra I. Lecture 2. Physics 342 Quantum Mechanics I
Physics 342 Lecture 2 Linear Algebra I Lecture 2 Physics 342 Quantum Mechanics I Wednesday, January 27th, 21 From separation of variables, we move to linear algebra Roughly speaking, this is the study
More information7.2 Steepest Descent and Preconditioning
7.2 Steepest Descent and Preconditioning Descent methods are a broad class of iterative methods for finding solutions of the linear system Ax = b for symmetric positive definite matrix A R n n. Consider
More informationNUMERICAL AND STATISTICAL COMPUTING (MCA-202-CR)
NUMERICAL AND STATISTICAL COMPUTING (MCA-202-CR) Autumn Session UNIT 1 Numerical analysis is the study of algorithms that uses, creates and implements algorithms for obtaining numerical solutions to problems
More informationLab I. 2D Motion. 1 Introduction. 2 Theory. 2.1 scalars and vectors LAB I. 2D MOTION 15
LAB I. 2D MOTION 15 Lab I 2D Motion 1 Introduction In this lab we will examine simple two-dimensional motion without acceleration. Motion in two dimensions can often be broken up into two separate one-dimensional
More information22 Approximations - the method of least squares (1)
22 Approximations - the method of least squares () Suppose that for some y, the equation Ax = y has no solutions It may happpen that this is an important problem and we can t just forget about it If we
More informationAnalysis Methods in Atmospheric and Oceanic Science
Analysis Methods in Atmospheric and Oceanic Science AOSC 652 Week 7, Day 1 13 Oct 2014 1 Student projects: 20% of the final grade: you will receive a numerical score for the project and final grade will
More informationPhysics 452 Lecture 33: A Particle in an E&M Field
Physics 452 Lecture 33: A Particle in an E&M Field J. Peatross In lectures 31 and 32, we considered the Klein-Gordon equation for a free particle. We would like to add a potential to the equation (since
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 5 Nonlinear Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction
More informationLab I. 2D Motion. 1 Introduction. 2 Theory. 2.1 scalars and vectors LAB I. 2D MOTION 15
LAB I. 2D MOTION 15 Lab I 2D Motion 1 Introduction In this lab we will examine simple two-dimensional motion without acceleration. Motion in two dimensions can often be broken up into two separate one-dimensional
More informationOutline. Scientific Computing: An Introductory Survey. Nonlinear Equations. Nonlinear Equations. Examples: Nonlinear Equations
Methods for Systems of Methods for Systems of Outline Scientific Computing: An Introductory Survey Chapter 5 1 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationToday s class. Numerical differentiation Roots of equation Bracketing methods. Numerical Methods, Fall 2011 Lecture 4. Prof. Jinbo Bi CSE, UConn
Today s class Numerical differentiation Roots of equation Bracketing methods 1 Numerical Differentiation Finite divided difference First forward difference First backward difference Lecture 3 2 Numerical
More informationPhysics 225 Relativity and Math Applications. Fall Unit 8 Differential Calculus: now available in 3D
Physics 225 Relativity and Math Applications Fall 2012 Unit 8 Differential Calculus: now available in 3D N.C.R. Makins University of Illinois at Urbana-Champaign 2010 Physics 225 8.2 8.2 Physics 225 8.3
More informationNumerical Methods in Informatics
Numerical Methods in Informatics Lecture 2, 30.09.2016: Nonlinear Equations in One Variable http://www.math.uzh.ch/binf4232 Tulin Kaman Institute of Mathematics, University of Zurich E-mail: tulin.kaman@math.uzh.ch
More informationModule II: Relativity and Electrodynamics
Module II: Relativity and Electrodynamics Lecture 2: Lorentz transformations of observables Amol Dighe TIFR, Mumbai Outline Length, time, velocity, acceleration Transformations of electric and magnetic
More informationScientific Computing. Roots of Equations
ECE257 Numerical Methods and Scientific Computing Roots of Equations Today s s class: Roots of Equations Bracketing Methods Roots of Equations Given a function f(x), the roots are those values of x that
More informationPhysically Based Modeling Differential Equation Basics
Physically Based Modeling Differential Equation Basics Andrew Witkin and David Baraff Pixar Animation Studios Please note: This document is 2001 by Andrew Witkin and David Baraff. This chapter may be freely
More informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More informationNearly Free Electron Gas model - I
Nearly Free Electron Gas model - I Contents 1 Free electron gas model summary 1 2 Electron effective mass 3 2.1 FEG model for sodium...................... 4 3 Nearly free electron model 5 3.1 Primitive
More informationSolving nonlinear equations (See online notes and lecture notes for full details) 1.3: Newton s Method
Solving nonlinear equations (See online notes and lecture notes for full details) 1.3: Newton s Method MA385 Numerical Analysis September 2018 (1/16) Sir Isaac Newton, 1643-1727, England. Easily one of
More informationIntroduction to Linear Algebra
Introduction to Linear Algebra Linear algebra is the algebra of vectors. In a course on linear algebra you will also learn about the machinery (matrices and reduction of matrices) for solving systems of
More informationMathematical Constraint on Functions with Continuous Second Partial Derivatives
1 Mathematical Constraint on Functions with Continuous Second Partial Derivatives J.D. Franson Physics Department, University of Maryland, Baltimore County, Baltimore, MD 15 Abstract A new integral identity
More informationModule 2: Reflecting on One s Problems
MATH55 Module : Reflecting on One s Problems Main Math concepts: Translations, Reflections, Graphs of Equations, Symmetry Auxiliary ideas: Working with quadratics, Mobius maps, Calculus, Inverses I. Transformations
More informationNumerical Methods Lecture 3
Numerical Methods Lecture 3 Nonlinear Equations by Pavel Ludvík Introduction Definition (Root or zero of a function) A root (or a zero) of a function f is a solution of an equation f (x) = 0. We learn
More informationPhysics 212E Spring 2004 Classical and Modern Physics. Computer Exercise #2
Physics 212E Spring 2004 Classical and Modern Physics Chowdary Computer Exercise #2 Launch Mathematica by clicking on the Start menu (lower left hand corner of the screen); from there go up to Science
More informationOrthogonality. 6.1 Orthogonal Vectors and Subspaces. Chapter 6
Chapter 6 Orthogonality 6.1 Orthogonal Vectors and Subspaces Recall that if nonzero vectors x, y R n are linearly independent then the subspace of all vectors αx + βy, α, β R (the space spanned by x and
More informationSection 1.1: Functions Given by Formulas
Section 1.1: Functions Given by Formulas Topics: Using function notation Domain Answering questions when given a formula functions Using TI to compute function values and using the Ans feature of the calculator
More informationSequences and infinite series
Sequences and infinite series D. DeTurck University of Pennsylvania March 29, 208 D. DeTurck Math 04 002 208A: Sequence and series / 54 Sequences The lists of numbers you generate using a numerical method
More informationDesigning Information Devices and Systems I Spring 2018 Lecture Notes Note Introduction to Linear Algebra the EECS Way
EECS 16A Designing Information Devices and Systems I Spring 018 Lecture Notes Note 1 1.1 Introduction to Linear Algebra the EECS Way In this note, we will teach the basics of linear algebra and relate
More informationElectromagnetic Waves Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space
Electromagnetic Waves 1 1. Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space 1 Retarded Potentials For volume charge & current = 1 4πε
More informationCodes for Partially Stuck-at Memory Cells
1 Codes for Partially Stuck-at Memory Cells Antonia Wachter-Zeh and Eitan Yaakobi Department of Computer Science Technion Israel Institute of Technology, Haifa, Israel Email: {antonia, yaakobi@cs.technion.ac.il
More informationCS123 INTRODUCTION TO COMPUTER GRAPHICS. Linear Algebra /34
Linear Algebra /34 Vectors A vector is a magnitude and a direction Magnitude = v Direction Also known as norm, length Represented by unit vectors (vectors with a length of 1 that point along distinct axes)
More informationPhysics 351 Monday, February 26, 2018
Physics 351 Monday, February 26, 2018 You just read the first half ( 10.1 10.7) of Chapter 10, which we ll probably start to discuss this Friday. The midterm exam (March 26) will cover (only!) chapters
More informationPhysics 342 Lecture 2. Linear Algebra I. Lecture 2. Physics 342 Quantum Mechanics I
Physics 342 Lecture 2 Linear Algebra I Lecture 2 Physics 342 Quantum Mechanics I Wednesday, January 3th, 28 From separation of variables, we move to linear algebra Roughly speaking, this is the study of
More information5 Finding roots of equations
Lecture notes for Numerical Analysis 5 Finding roots of equations Topics:. Problem statement. Bisection Method 3. Newton s Method 4. Fixed Point Iterations 5. Systems of equations 6. Notes and further
More informationInterpolation. 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter
Key References: Interpolation 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter 6. 2. Press, W. et. al. Numerical Recipes in C, Cambridge: Cambridge University Press. Chapter 3
More information- HH Why Can Light Propagate in Vacuum? Hsiu-Hau Lin (Apr 1, 2014)
- HH0124 - Why Can Light Propagate in Vacuum? Hsiu-Hau Lin hsiuhau@phys.nthu.edu.tw (Apr 1, 2014) Sunlight is vital for most creatures in this warm and humid planet. From interference experiments, scientists
More informationProperties of Real Numbers
Properties of Real Numbers Essential Understanding. Relationships that are always true for real numbers are called properties, which are rules used to rewrite and compare expressions. Two algebraic expressions
More informationPHYS 328 HOMEWORK # 5
PHYS 38 HOMEWORK # 5 Solutions 1. This first uestion is one of several that will employ Stirling' s approximation to obtain analytic expressions that will help us understand various thermodynamic systems.
More informationMath 56 Homework 1 Michael Downs. ne n 10 + ne n (1)
. Problem (a) Yes. The following equation: ne n + ne n () holds for all n R but, since we re only concerned with the asymptotic behavior as n, let us only consider n >. Dividing both sides by n( + ne n
More information(One Dimension) Problem: for a function f(x), find x 0 such that f(x 0 ) = 0. f(x)
Solving Nonlinear Equations & Optimization One Dimension Problem: or a unction, ind 0 such that 0 = 0. 0 One Root: The Bisection Method This one s guaranteed to converge at least to a singularity, i not
More informationPhysics 111. Thursday, October 07, Conservation of Momentum. Kinetic Energy
ics Thursday, ober 07, 2004 Ch 9: Ch 7: Momentum Conservation of Momentum Work Kinetic Energy Announcements Help this week: Wednesday, 8-9 pm in NSC 118/119 Sunday, 6:30-8 pm in CCLIR 468 Announcements
More informationA Greedy Framework for First-Order Optimization
A Greedy Framework for First-Order Optimization Jacob Steinhardt Department of Computer Science Stanford University Stanford, CA 94305 jsteinhardt@cs.stanford.edu Jonathan Huggins Department of EECS Massachusetts
More informationAugust Progress Report
PATH PREDICTION FOR AN EARTH-BASED DEMONSTRATION BALLOON FLIGHT DANIEL BEYLKIN Mentor: Jerrold Marsden Co-Mentors: Claire Newman and Philip Du Toit August Progress Report. Progress.. Discrete Mechanics
More informationSOLUTIONS to Exercises from Optimization
SOLUTIONS to Exercises from Optimization. Use the bisection method to find the root correct to 6 decimal places: 3x 3 + x 2 = x + 5 SOLUTION: For the root finding algorithm, we need to rewrite the equation
More informationSolutions. Part I Logistic regression backpropagation with a single training example
Solutions Part I Logistic regression backpropagation with a single training example In this part, you are using the Stochastic Gradient Optimizer to train your Logistic Regression. Consequently, the gradients
More informationPreliminaries. Definition: The Euclidean dot product between two vectors is the expression. i=1
90 8 80 7 70 6 60 0 8/7/ Preliminaries Preliminaries Linear models and the perceptron algorithm Chapters, T x + b < 0 T x + b > 0 Definition: The Euclidean dot product beteen to vectors is the expression
More informationHidden Markov Models (recap BNs)
Probabilistic reasoning over time - Hidden Markov Models (recap BNs) Applied artificial intelligence (EDA132) Lecture 10 2016-02-17 Elin A. Topp Material based on course book, chapter 15 1 A robot s view
More informationMultistep Methods for IVPs. t 0 < t < T
Multistep Methods for IVPs We are still considering the IVP dy dt = f(t,y) t 0 < t < T y(t 0 ) = y 0 So far we have looked at Euler s method, which was a first order method and Runge Kutta (RK) methods
More informationRoot Finding Convergence Analysis
Root Finding Convergence Analysis Justin Ross & Matthew Kwitowski November 5, 2012 There are many different ways to calculate the root of a function. Some methods are direct and can be done by simply solving
More informationRobotics I. April 1, the motion starts and ends with zero Cartesian velocity and acceleration;
Robotics I April, 6 Consider a planar R robot with links of length l = and l =.5. he end-effector should move smoothly from an initial point p in to a final point p fin in the robot workspace so that the
More information