Computational Physics HW2
|
|
- Harvey Burns
- 6 years ago
- Views:
Transcription
1 Computational Physics HW2 Luke Bouma July 27, Plotting experimental data 1.1 Plotting sunspots.txt in Python: The figure above is the output from from numpy import l o a d t x t data = l o a d t x t ( sunspots. txt, f l o a t ) time = data [ :, 0 ] sunspotn = data [ :, 1 ] p l t. x l a b e l ( Time [ month #] ) p l t. y l a b e l ( Number o f sunspots ) 1
2 p l t. p l o t ( time, sunspotn ) 1.2 Modifying the program to display the first 1000 data points: Where the only output lines that were changed are: time = data [ 0 : , 0 ] sunspotn = data [ 0 : , 1 ] 1.3 Calculating running average of the data. Define the running average to be Y k = 1 2r r m= r y k+m, where r = 5 and y k are the sunspot numbers. This wound up being a bit messier than I wanted, but it still makes good sense: #Import from numpy import loadtxt, sum, zeros, s i z e data = l o a d t x t ( sunspots. txt, f l o a t ) time = data [ :, 0 ] sunspotn = data [ :, 1 ] 2
3 #Smoothing r = 5 Y = z e r o s ( s i z e ( sunspotn ) ) for k in range ( s i z e ( sunspotn ) ) : i f k < 5 or k > 3100: Y[ k ] = sunspotn [ k ] else : for m in range ( k 5,k +5,1): Y[ k ] += 1/(2 r ) sunspotn [m] #P l o t t i n g monthn = 1000 p l t. x l a b e l ( Time [ month #] ) p l t. y l a b e l ( Number o f sunspots ) p l t. p l o t ( time [ 0 : monthn ], sunspotn [ 0 : monthn ], b, time [ 0 : monthn ], Y[ 0 : monthn ], r s ) 3
4 Figure 1: Daw, I love you too! 2 Curve plotting 2.1 Deltoid curve plot Plot the parametric equations x = 2 cos θ + cos 2θ, y = 2 sin θ sin2θ where 0 θ < 2π. Take a set of values of θ between zero and 2π, and calculate x and y for each from the equations above, then plot y as a function of x. We get Fig. 1 from the script from math import sin, cos, pi from numpy import l i n s p a c e, s i z e, z e r o s l e n = 100 x, y = z e r o s ( l e n ), z e r o s ( l e n ) theta = l i n s p a c e (0,2 pi, l e n ) for j in range ( 0, s i z e ( theta ) ) : x [ j ] = 2 cos ( theta [ j ] ) + cos (2 theta [ j ] ) y [ j ] = 2 s i n ( theta [ j ] ) + s i n (2 theta [ j ] ) p l t. x l a b e l ( x ) p l t. y l a b e l ( y ) p l t. p l o t ( x, y, bs ) 4
5 2.2 Plotting the Galilean spiral, r = θ 2 for 0 θ 24π. is generated from from math import sin, cos, pi from numpy import l i n s p a c e, s i z e, z e r o s l e n = 1000 x, y, r = z e r o s ( l e n ), z e r o s ( l e n ), z e r o s ( l e n ) theta = l i n s p a c e (0,10 pi, l e n ) for j in range ( 0, s i z e ( theta ) ) : r [ j ] = theta [ j ] 2 x [ j ] = r [ j ] cos ( theta [ j ] ) y [ j ] = r [ j ] s i n ( theta [ j ] ) p l t. x l a b e l ( x ) p l t. y l a b e l ( y ) p l t. p l o t ( x, y, r ) 2.3 Fey s function polar plot Using the script from math import sin, cos, exp, pi from numpy import l i n s p a c e, s i z e, z e r o s 5
6 l e n = x, y, r = z e r o s ( l e n ), z e r o s ( l e n ), z e r o s ( l e n ) theta = l i n s p a c e (0,24 pi, l e n ) for j in range ( 0, s i z e ( theta ) ) : r [ j ] = exp ( cos ( theta [ j ] ) ) 2 cos (4 theta [ j ] ) + \ s i n ( theta [ j ] / 1 2 ) 5 x [ j ] = r [ j ] cos ( theta [ j ] ) y [ j ] = r [ j ] s i n ( theta [ j ] ) p l t. x l a b e l ( x ) p l t. y l a b e l ( y ) p l t. p l o t ( x, y, r ) we generate 3 Plotting atomic lattice (VPython) This is the one that I m skipping/avoiding because I already did it in the notes, to generate 6
7 which is an approximate NaCl lattice. Doing a FCC (face-centered cubic) lattice would be a question of adding another atom 2a/2 offset on the diagonal from the central atom. 4 Deterministic chaos and the Feigenbaum plot Consider the logistic map, defined by x0 = rx(1 x). For a given constant r take a value of x, say x = 1/2, and feed it into the RHS of this equation to get x0. Then take the value you get, feed it back in, and so on. I.e., you iterate. One of three things will happen: 1. The value settles down to a fixed number, and stays there. For instance, x = 0 is one such fixed point of the logistic map. 2. You never settle to a single value, but instead settle into a periodic pattern, rotating around a set number of values, repeating them in sequence forever. This is a limit cycle. 3. It goes crazy. It generates a seemingly random sequence of numbers that appear to have no rhyme or reason to them at all. This is deterministic chaos, in the sense that it really does look chaotic, but deterministic because we know that the values have a rule behind them. The behavior is determined, it just isn t obvious that this is the case. We write a program to calculate and display the behavior of the logistic map: import m a t p l o t l i b. p y p l o t a s p l t from numpy import l i n s p a c e, s i z e #Returns r e s u l t o f n o p e r a t i o n s o f l o g i s t i c map on x, g i v e n r def l o g i s t i c M a p ( r, x, n ) : i, xprime = 0, 0. 7
8 while i < n : xprime = r x (1 x ) x = xprime i += 1 return xprime #Returns l i s t o f l a t t e r 1000 mappings to p l o t def g e t X l i s t ( r, x ) : i = 0 while i < 250: x. append ( l o g i s t i c M a p ( r, x [ s i z e ( x ) 1], 1 ) ) i += 1 return x xmin, xmax = , 3.86 for r in l i n s p a c e ( xmin, xmax, ) : print ( r ) x = [ l o g i s t i c M a p ( r, 1/2, ) ] x = g e t X l i s t ( r, x ) for i in range ( s i z e ( x ) ) : p l t. s c a t t e r ( r, x [ i ], 1) p l t. x l a b e l ( r ) p l t. y l a b e l ( x ( Orbit ) ) p l t. t i t l e ( L o g i s t i c Map ( Feigenbaum Plot ) ) p l t. a x i s ( [ xmin, xmax, 0, 1 ] ) which, once plotted in Figs. 2 and 3 gives pretty nice results. Answering the questions: a) Fixed points look like single orbits (values of x) corresponding to single values of r. Limit cycles are when two or more values of x correspond to a single value of r. Chaos is when a seemingly random large number of x points are given by a single value of r. b) The switch from orderly (fixed points, limit cycles) to chaotic happens roughly at r = (obviously, wiki). Aside: This deterministic chaos is seen in more complex physical systems too, most notably fluid dynamics and the weather. We get an idea through these plots of what chaotic behavior means: change initial conditions by just a little, and you ll get a rapid onset of highly diverse, hard-to-predict results. The classic example is Edward Lorenz s butterfly effect (although scifi writer Ray Bradbury may have suggested the idea in a story 25 years earlier). We can definitely go further than these two figures though. Generating them with roughly 10 6 points took upwards of ten minutes on my laptop, and we should be able to do much better. Noting the hint at the end of the problem to vectorize the code, we rewrite the whole program and find nicer results: from numpy import l i n s p a c e, array, vstack, f u l l, copy #Returns r e s u l t o f n o p e r a t i o n s o f l o g i s t i c map on x, given r 8
9 Figure 2: The logistic map, x n+1 = rx n (1 x n ), for 2.8 r 4.0. def l o g i s t i c M a p ( r, x, n ) : i, xprime = 0, 0. while i < n : xprime = r x (1 x ) x = xprime i += 1 return xprime #Returns f u l l matrix o f x v a l u e s with each map i t e r a t i o n f o r numx def getmatrix ( r, x, numx) : i = 0 while i < numx: i f i == 0 : #obnoxious, but seemed necessary x = vstack ( ( x, l o g i s t i c M a p ( r, x, 1 ) ) ) else : x = vstack ( ( x [ 0 : i, : ], l o g i s t i c M a p ( r, x [ i 1, : ], 1 ) ) ) i += 1 return x numr = 1000 #number o f r v a l u e s numx = 1000 #number o f x v a l u e s f o r each r v a l u e rmin, rmax = 2. 2, 3 #r between min and max alp = 0.02 #from 0 to 1, changes c o n t r a s t o f p l o t t e d p t s 9
10 Figure 3: The logistic map in the first large island of stability, which begins at r = l i n s p a c e ( rmin, rmax, numr) x = array ( [ 1 / 2 ] numr) x = l o g i s t i c M a p ( r, x, 1000) #s t a b i l i z e s t a r t i n g v a l u e s over 1000 i t e r a t i o n s x = getmatrix ( r, x, numx) #g e t the data o f x o r b i t s f o r each r ( v e c t o r i z e d ) print ( Got matrix, now p l o t t i n g. ) rarr = f u l l ( (numx,numr), 1. ) #used to v e c t o r i z e s c a t t e r p l o t for i in range (numx) : rarr [ i, : ] = copy ( r ) p l t. x l a b e l ( r ) p l t. y l a b e l ( x ( Orbit ) ) p l t. t i t l e ( L o g i s t i c Map ( Feigenbaum Plot ) ) p l t. a x i s ( [ rmin, rmax, 0, 1 ] ) p l t. s c a t t e r ( rarr, x, s =4, alpha=alp ) This program runs like, at least a factor of 100 faster than the other. So obviously we push it further, and make more plots: 10
11 Figure 4: More points. Figure 5: Getting artistic. 11
12 5 Least squares and the photoelectric effect This is the method of least squares. If we have N data points with coordinates (x i, y i ), then the sum of the squares of the distances between some guess line y = mx + c and the points is χ 2 = N (mx i + c y i ) 2 = i=1 N m 2 x 2 i + c 2 + yi 2 + 2mx i c 2mx i y i 2cy i. i=1 The least squares fit is the line that minimizes χ 2. Find this minimum by differentiating with respect to both m and c, and setting the derivatives to zero: N N N m x 2 i + c x i x i y i = 0 i=1 i=1 i=1 N N m x i + cn y i = 0. i=1 i=1 For convenience, define E x = 1 N N x i, i=1 E y = 1 N N i=1 y i, E xx = 1 N N x 2 i, i=1 E xy = 1 N N x i y i, i=1 so that we can rewrite our equations Solving simultaneously for m and c gives m = E xy E x E y E xx E 2 x me xx + ce x = E xy me x + c = E y., c = E xxe y E x E xy E xx Ex 2. These are the equations for the least-squares fit of a straight line to N data points. They give the values of m and c for the last line that best fits the given data. 5.1 Millikan: read data points, make graph from numpy import l o a d t x t data = l o a d t x t ( m i l l i k a n. txt, f l o a t ) x = data [ :, 0 ] y = data [ :, 1 ] p l t. p l o t ( x, y, k. ) 5.2 Find E y, E x x, E x y and print m,c of best-fit line 12
13 N = s i z e ( x ) E x = 1/N sum( x ) E y = 1/N sum( y ) E xx = 1/N sum( x 2) E xy = 1/N sum( x y ) m = ( E xy E x E y ) / ( E xx E x 2) c = ( E xx E y E x E xy ) / ( E xx E x 2) print ( m=%.3 f % m, c=%.3 f % c ) 5.3 Draw straight line Draw the straight line by evaluating mx i + c using the values of m, c from the previous part. Store these values to a list or array, then graph this array as a solid line. This gives a program: from numpy import loadtxt, s i z e data = l o a d t x t ( m i l l i k a n. txt, f l o a t ) x = data [ :, 0]/1 e15 y = data [ :, 1 ] N = s i z e ( x ) E x = sum( x ) / N E y = sum( y ) / N E xx = sum( x 2) / N E xy = sum( x y ) / N m = ( E xy E x E y ) / ( E xx E x E x ) c = ( E xx E y E x E xy ) / ( E xx E x 2) b e s t F i t = m x + c print ( m=%.3 f % m, c=%.3 f % c ) p l t. p l o t ( x, y, k. ) p l t. p l o t ( x, b e s t F i t ) 13
14 5.4 Calculating Planck s constant Given that this data is actually of a photoelectric experiment, where the voltage required to stop ejected electrons, V, is measured as a function of ν, the frequency of incoming light, so our equation y = mx + c actually is V = h e ν φ, so that y = V, m = h/e, c = φ. We are given that e = C, and since we have m = VHz 1, we can calculate h = m e = = VCs = Js which is indeed within a few percent of the literature value. 6 Trapezoidal rule and Romberg rule for integrals Evaluate I = 1 0 sin 2 100x dx 6.1 Use the trapezoidal rule to calculate it to accuracy ε = 10 6 The approximation used in the trapezoidal rule integrating f over [a, b] is ( 1 I h 2 f(a) + 1 N ) 2 f(b) + f(a + kh), for ε = 1 12 h2 (f (a) f (b)). For us, f(x) = sin 2 100x, so f (x) = d dx (sin 100x sin 100x) = 10x 1/2 cos 100x sin 100x, which simplified is 5 sin (20 x)/ x and to avoid dividing by zero, 5 sin(20 x)/ x = 5 ( 20 x (20 x) 3 x 3! k=1 + (20 x) 5 5! (20 x) 7 7! + (20 x) 9 9! ( = x x x x x 5 + O(x 5 ) 3! 5! 7! 9! 11! ) + O(x 11/2 ) which we implement to be accurate up to O(x 6 ) to avoid the division by zero. The program we write is from numpy import sin, s q r t def f a c t o r i a l ( x ) : i f x == 1 : return 1 else : return x f a c t o r i a l ( x 1) ) 14
15 def f ( x ) : return s i n ( s q r t (100 x )) 2 def fprime ( x ) : return 5 ( x/ f a c t o r i a l ( 3 ) x 2/ f a c t o r i a l ( 5 ) 20 7 x 3/ f a c t o r i a l ( 7 ) x 4/ f a c t o r i a l ( 9 ) x 5/ f a c t o r i a l ( 1 1 ) ) def eps (h, a, b ) : return h 2 ( fprime ( a ) fprime ( b ) ) / 12 a, b, I = 0, 1, 0. for i in range ( 5 0 ) : N = 2 i h = (b a )/N e r r o r = eps (h, a, b ) print ( i, h, e r r o r ) i f e r r o r < 1e 6: I = h ( f ( a )/2 + f ( b )/2) for k in range ( 1,N) : I += h f ( a+k h ) break print ( I ) which results I = , good to Mathematica s result to 10 decimals! 6.2 Romberg technique Honestly I can t be bothered here. Gaussian stuff in any case. Looks pretty boring, and we can just do the higher-order 7 Heat capacity of a solid Debye s theory of solids gives the heat capacity of a solid at temperature T to be C V = 9V ρk B θd /T 0 x 4 e x (e x 1) 2 dx, for V the volume of the solid, ρ the number density of atoms, and θ D the Debye temperature, a property of solids that depends on their density and the speed of sound. a) Write a function cv(t) that calculates C V for a given value of the temperature, for a sample of 1000 cm 3 = 10 3 m 3 of solid Al, which has number density ρ = m 3, and θ D = 428K. Use Gaussian quadrature to evaluate the integral, with N = 50 sample points. 15
16 from numpy import exp, l i n s p a c e from gaussxw import gaussxw, gaussxwab def f ( x ) : return x 4 exp ( x ) / ( exp ( x) 1) 2 def cv (T) : V = 1e 3 #m 3 rho = e28 #1/m 3 kb = e 23 #J/K thetad = 428 p r e f a c = 9 V rho kb (T/ thetad ) 3 #Get p o i n t s and w e i g h t s f o r i n t e g r a l N = 50 a, b = 0., thetad /T xp, wp = gaussxwab (N, a, b ) #Perform i n t e g r a t i o n s = 0. for k in range (N) : s += wp [ k ] f ( xp [ k ] ) #p r i n t ( s ) return p r e f a c s T = l i n s p a c e ( 5, 5 0 0, ) c v L i s t = [ ] for i in range ( ) : c v L i s t. append ( cv (T[ i ] ) ) p l t. p l o t (T, c v L i s t ) p l t. x l a b e l ( Temperature $ [K] $ ) p l t. y l a b e l ( Heat c a p a c i t y $C V$ $ [ J /( kg\ cdot K) ] $ ) p l t. t i t l e ( $C v (T) $ o f Aluminum, Debye P r e d i c t i o n ) gives the figure 16
17 17
PHYSICS 3266 SPRING 2016
PHYSICS 3266 SPRIG 2016 Each problem is worth 5 points as discussed in the syllabus. For full credit you must include in your solution a copy of your program (well commented and listed any students that
More informationComputational Physics HW3
Computational Physics HW3 Luke Bouma August 9, 215 1 The Stefan-Boltzmann constant In the angular frequency interval ω to ω + dω, black bodies of unit area radiate a thermal energy per second equal to
More informationPhysics 411: Homework 3
Physics 411: Homework 3 Because of the cancellation of class on January 28, this homework is a double-length homework covering two week s material, and you have two weeks to do it. It is due in class onthursday,
More informationSolutions to homework assignment #7 Math 119B UC Davis, Spring for 1 r 4. Furthermore, the derivative of the logistic map is. L r(x) = r(1 2x).
Solutions to homework assignment #7 Math 9B UC Davis, Spring 0. A fixed point x of an interval map T is called superstable if T (x ) = 0. Find the value of 0 < r 4 for which the logistic map L r has a
More informationLESSON 21: DIFFERENTIALS OF MULTIVARIABLE FUNCTIONS MATH FALL 2018
LESSON 21: DIFFERENTIALS OF MULTIVARIABLE FUNCTIONS MATH 16020 FALL 2018 ELLEN WELD 1. Quick Review of Differentials Ex 1. Consider the function f(x) x. We know that f(9) 9 3, but what is f(9.1) 9.1? Obviously,
More informationLESSON 21: DIFFERENTIALS OF MULTIVARIABLE FUNCTIONS OCTOBER 18, 2017
LESSON 21: DIFFERENTIALS OF MULTIVARIABLE FUNCTIONS OCTOBER 18, 2017 Today we do a quick review of differentials for functions of a single variable and then discuss how to extend this notion to functions
More informationElectro Magnetic Field Dr. Harishankar Ramachandran Department of Electrical Engineering Indian Institute of Technology Madras
Electro Magnetic Field Dr. Harishankar Ramachandran Department of Electrical Engineering Indian Institute of Technology Madras Lecture - 7 Gauss s Law Good morning. Today, I want to discuss two or three
More informationIntroduction. Prediction MATH February 2017
21 February 2017 Predicting the future is very difficult, especially if it s about the future. Niels Bohr Can we say what is going to happen: in the next minute? tomorrow? next year? Predicting the future
More informationMath Problem Set #3 Solution 19 February 2001
Math 203-04 Problem Set #3 Solution 19 February 2001 Exercises: 1. B & D, Section 2.3, problem #3. In your answer, give both exact values and decimal approximations for the amount of salt in the tank at
More informationMATH ASSIGNMENT 03 SOLUTIONS
MATH444.0 ASSIGNMENT 03 SOLUTIONS 4.3 Newton s method can be used to compute reciprocals, without division. To compute /R, let fx) = x R so that fx) = 0 when x = /R. Write down the Newton iteration for
More informationWhy are Discrete Maps Sufficient?
Why are Discrete Maps Sufficient? Why do dynamical systems specialists study maps of the form x n+ 1 = f ( xn), (time is discrete) when much of the world around us evolves continuously, and is thus well
More informationIntroduction to Dynamical Systems Basic Concepts of Dynamics
Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic
More informationMidterm Exam I for Math 110, Spring Solutions
Name: Midterm Exam I for Math 110, Spring 2015 Solutions Problem Points Score 1 12 2 12 3 12 4 12 5 12 6 12 7 12 8 12 9 12 10 12 Total 120 You have ninety minutes for this exam. Please show ALL your work
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationUsing Matlab to integrate Ordinary Differential Equations (ODEs)
Using Matlab to integrate Ordinary Differential Equations (ODEs) Erica McEvoy (Dated: June 17, 2009) 1. INTRODUCTION Ordinary differential equations tend to arise whenever you need to model changing quantities
More informationTropical Polynomials
1 Tropical Arithmetic Tropical Polynomials Los Angeles Math Circle, May 15, 2016 Bryant Mathews, Azusa Pacific University In tropical arithmetic, we define new addition and multiplication operations on
More informationIterated Functions. Tom Davis November 5, 2009
Iterated Functions Tom Davis tomrdavis@earthlink.net http://www.geometer.org/mathcircles November 5, 2009 Abstract In this article we will examine various properties of iterated functions. If f(x) is a
More informationExtracting beauty from chaos
1997 2004, Millennium Mathematics Project, University of Cambridge. Permission is granted to print and copy this page on paper for non commercial use. For other uses, including electronic redistribution,
More informationMaps and differential equations
Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number
More informationR x n. 2 R We simplify this algebraically, obtaining 2x n x n 1 x n x n
Math 42 Homework 4. page 3, #9 This is a modification of the bisection method. Write a MATLAB function similar to bisect.m. Here, given the points P a a,f a and P b b,f b with f a f b,we compute the point
More informationMath 12: Discrete Dynamical Systems Homework
Math 12: Discrete Dynamical Systems Homework Department of Mathematics, Harvey Mudd College Several of these problems will require computational software to help build our insight about discrete dynamical
More informationf[x_, μ_] := 4. μ... Nest[f[#,...] &,...,...] Plot[{x, f[...]}, {x, 0, 1}, AspectRatio Equal]
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 In this exercise, we use renormalization-group and scaling methods to study the onset of chaos. There are several routes by which a dynamical system can start exhibiting chaotic
More informationMA 1128: Lecture 19 4/20/2018. Quadratic Formula Solving Equations with Graphs
MA 1128: Lecture 19 4/20/2018 Quadratic Formula Solving Equations with Graphs 1 Completing-the-Square Formula One thing you may have noticed when you were completing the square was that you followed the
More informationNUMERICAL ANALYSIS WEEKLY OVERVIEW
NUMERICAL ANALYSIS WEEKLY OVERVIEW M. AUTH 1. Monday 28 August Students are encouraged to download Anaconda Python. Anaconda is a version of Python that comes with some numerical packages (numpy and matplotlib)
More informationGoals: Second-order Linear Equations Linear Independence of Solutions and the Wronskian Homogeneous DEs with Constant Coefficients
Week #3 : Higher-Order Homogeneous DEs Goals: Second-order Linear Equations Linear Independence of Solutions and the Wronskian Homogeneous DEs with Constant Coefficients 1 Second-Order Linear Equations
More informationMATH 2053 Calculus I Review for the Final Exam
MATH 05 Calculus I Review for the Final Exam (x+ x) 9 x 9 1. Find the limit: lim x 0. x. Find the limit: lim x + x x (x ).. Find lim x (x 5) = L, find such that f(x) L < 0.01 whenever 0 < x
More informationLectures about Python, useful both for beginners and experts, can be found at (http://scipy-lectures.github.io).
Random Matrix Theory (Sethna, "Entropy, Order Parameters, and Complexity", ex. 1.6, developed with Piet Brouwer) 2016, James Sethna, all rights reserved. This is an ipython notebook. This hints file is
More informationSolving Quadratic & Higher Degree Equations
Chapter 7 Solving Quadratic & Higher Degree Equations Sec 1. Zero Product Property Back in the third grade students were taught when they multiplied a number by zero, the product would be zero. In algebra,
More informationLecture 1: Period Three Implies Chaos
Math 7h Professor: Padraic Bartlett Lecture 1: Period Three Implies Chaos Week 1 UCSB 2014 (Source materials: Period three implies chaos, by Li and Yorke, and From Intermediate Value Theorem To Chaos,
More informationCOSMOS: Making Robots and Making Robots Intelligent Lecture 3: Introduction to discrete-time dynamics
COSMOS: Making Robots and Making Robots Intelligent Lecture 3: Introduction to discrete-time dynamics Jorge Cortés and William B. Dunbar June 3, 25 Abstract In this and the coming lecture, we will introduce
More informationNumerical Methods - Preliminaries
Numerical Methods - Preliminaries Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Preliminaries 2013 1 / 58 Table of Contents 1 Introduction to Numerical Methods Numerical
More informationMath 113 HW #10 Solutions
Math HW #0 Solutions 4.5 4. Use the guidelines of this section to sketch the curve Answer: Using the quotient rule, y = x x + 9. y = (x + 9)(x) x (x) (x + 9) = 8x (x + 9). Since the denominator is always
More informationMath 1310 Section 4.1: Polynomial Functions and Their Graphs. A polynomial function is a function of the form ...
Math 1310 Section 4.1: Polynomial Functions and Their Graphs A polynomial function is a function of the form... where 0,,,, are real numbers and n is a whole number. The degree of the polynomial function
More informationSection 2.8: The Power Chain Rule
calculus sin frontera Section 2.8: The Power Chain Rule I want to see what happens when I take a function I know, like g(x) = x 2, and raise it to a power, f (x) = ( x 2 ) a : f (x) = ( x 2 ) a = ( x 2
More informationMATH2321, Calculus III for Science and Engineering, Fall Name (Printed) Print your name, the date, and then sign the exam on the line
MATH2321, Calculus III for Science and Engineering, Fall 218 1 Exam 2 Name (Printed) Date Signature Instructions STOP. above. Print your name, the date, and then sign the exam on the line This exam consists
More informationLab 1: Iterative Methods for Solving Linear Systems
Lab 1: Iterative Methods for Solving Linear Systems January 22, 2017 Introduction Many real world applications require the solution to very large and sparse linear systems where direct methods such as
More informationLesson 6-1: Relations and Functions
I ll bet you think numbers are pretty boring, don t you? I ll bet you think numbers have no life. For instance, numbers don t have relationships do they? And if you had no relationships, life would be
More information3 Algebraic Methods. we can differentiate both sides implicitly to obtain a differential equation involving x and y:
3 Algebraic Methods b The first appearance of the equation E Mc 2 in Einstein s handwritten notes. So far, the only general class of differential equations that we know how to solve are directly integrable
More information7.2 Linear equation systems. 7.3 Linear least square fit
72 Linear equation systems In the following sections, we will spend some time to solve linear systems of equations This is a tool that will come in handy in many di erent places during this course For
More informationBasic Theory of Dynamical Systems
1 Basic Theory of Dynamical Systems Page 1 1.1 Introduction and Basic Examples Dynamical systems is concerned with both quantitative and qualitative properties of evolution equations, which are often ordinary
More informationMS 3011 Exercises. December 11, 2013
MS 3011 Exercises December 11, 2013 The exercises are divided into (A) easy (B) medium and (C) hard. If you are particularly interested I also have some projects at the end which will deepen your understanding
More informationKey Point. The nth order linear homogeneous equation with constant coefficients
General Solutions of Higher-Order Linear Equations In section 3.1, we saw the following fact: Key Point. The nth order linear homogeneous equation with constant coefficients a n y (n) +... + a 2 y + a
More informationLecture 3: Sizes of Infinity
Math/CS 20: Intro. to Math Professor: Padraic Bartlett Lecture 3: Sizes of Infinity Week 2 UCSB 204 Sizes of Infinity On one hand, we know that the real numbers contain more elements than the rational
More informationPractice Questions From Calculus II. 0. State the following calculus rules (these are many of the key rules from Test 1 topics).
Math 132. Practice Questions From Calculus II I. Topics Covered in Test I 0. State the following calculus rules (these are many of the key rules from Test 1 topics). (Trapezoidal Rule) b a f(x) dx (Fundamental
More informationPHYS 7411 Spring 2015 Computational Physics Homework 3
PHYS 7411 Spring 2015 Computational Physics Homework 3 Due by 3:00pm in Nicholson 447 on 13 March 2015, Friday the 13th (Any late assignments will be penalized in the amount of 25% per day late. Any copying
More informationMath 231E, Lecture 13. Area & Riemann Sums
Math 23E, Lecture 3. Area & Riemann Sums Motivation for Integrals Question. What is an integral, and why do we care? Answer. A tool to compute a complicated expression made up of smaller pieces. Example.
More informationLorenz Equations. Lab 1. The Lorenz System
Lab 1 Lorenz Equations Chaos: When the present determines the future, but the approximate present does not approximately determine the future. Edward Lorenz Lab Objective: Investigate the behavior of a
More informationChapter 1: January 26 January 30
Chapter : January 26 January 30 Section.7: Inequalities As a diagnostic quiz, I want you to go through the first ten problems of the Chapter Test on page 32. These will test your knowledge of Sections.
More informationMath Numerical Analysis Mid-Term Test Solutions
Math 400 - Numerical Analysis Mid-Term Test Solutions. Short Answers (a) A sufficient and necessary condition for the bisection method to find a root of f(x) on the interval [a,b] is f(a)f(b) < 0 or f(a)
More informationWhere Is Newton Taking Us? And How Fast?
Name: Where Is Newton Taking Us? And How Fast? In this activity, you ll use a computer applet to investigate patterns in the way the approximations of Newton s Methods settle down to a solution of the
More information5.2 Polynomial Operations
5.2 Polynomial Operations At times we ll need to perform operations with polynomials. At this level we ll just be adding, subtracting, or multiplying polynomials. Dividing polynomials will happen in future
More informationPolynomial functions right- and left-hand behavior (end behavior):
Lesson 2.2 Polynomial Functions For each function: a.) Graph the function on your calculator Find an appropriate window. Draw a sketch of the graph on your paper and indicate your window. b.) Identify
More informationIdeas from Vector Calculus Kurt Bryan
Ideas from Vector Calculus Kurt Bryan Most of the facts I state below are for functions of two or three variables, but with noted exceptions all are true for functions of n variables..1 Tangent Line Approximation
More informationCHAPTER 2 POLYNOMIALS KEY POINTS
CHAPTER POLYNOMIALS KEY POINTS 1. Polynomials of degrees 1, and 3 are called linear, quadratic and cubic polynomials respectively.. A quadratic polynomial in x with real coefficient is of the form a x
More informationAn Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations
An Overly Simplified and Brief Review of Differential Equation Solution Methods We will be dealing with initial or boundary value problems. A typical initial value problem has the form y y 0 y(0) 1 A typical
More informationPartial Fractions. June 27, In this section, we will learn to integrate another class of functions: the rational functions.
Partial Fractions June 7, 04 In this section, we will learn to integrate another class of functions: the rational functions. Definition. A rational function is a fraction of two polynomials. For example,
More informationInfinity Unit 2: Chaos! Dynamical Systems
Infinity Unit 2: Chaos! Dynamical Systems Iterating Linear Functions These questions are about iterating f(x) = mx + b. Seed: x 1. Orbit: x 1, x 2, x 3, For each question, give examples and a symbolic
More informationMain topics for the First Midterm Exam
Main topics for the First Midterm Exam The final will cover Sections.-.0, 2.-2.5, and 4.. This is roughly the material from first three homeworks and three quizzes, in addition to the lecture on Monday,
More informationStudy Unit 2 : Linear functions Chapter 2 : Sections and 2.6
1 Study Unit 2 : Linear functions Chapter 2 : Sections 2.1 2.4 and 2.6 1. Function Humans = relationships Function = mathematical form of a relationship Temperature and number of ice cream sold Independent
More informationSlope Fields: Graphing Solutions Without the Solutions
8 Slope Fields: Graphing Solutions Without the Solutions Up to now, our efforts have been directed mainly towards finding formulas or equations describing solutions to given differential equations. Then,
More informationWhat if the characteristic equation has a double root?
MA 360 Lecture 17 - Summary of Recurrence Relations Friday, November 30, 018. Objectives: Prove basic facts about basic recurrence relations. Last time, we looked at the relational formula for a sequence
More informationPractice problems from old exams for math 132 William H. Meeks III
Practice problems from old exams for math 32 William H. Meeks III Disclaimer: Your instructor covers far more materials that we can possibly fit into a four/five questions exams. These practice tests are
More informationChaos and Liapunov exponents
PHYS347 INTRODUCTION TO NONLINEAR PHYSICS - 2/22 Chaos and Liapunov exponents Definition of chaos In the lectures we followed Strogatz and defined chaos as aperiodic long-term behaviour in a deterministic
More informationMath 234 Exam 3 Review Sheet
Math 234 Exam 3 Review Sheet Jim Brunner LIST OF TOPIS TO KNOW Vector Fields lairaut s Theorem & onservative Vector Fields url Divergence Area & Volume Integrals Using oordinate Transforms hanging the
More informationUnit 2 Rational Functionals Exercises MHF 4UI Page 1
Unit 2 Rational Functionals Exercises MHF 4UI Page Exercises 2.: Division of Polynomials. Divide, assuming the divisor is not equal to zero. a) x 3 + 2x 2 7x + 4 ) x + ) b) 3x 4 4x 2 2x + 3 ) x 4) 7. *)
More informationIterated Functions. Tom Davis August 2, 2017
Iterated Functions Tom Davis tomrdavis@earthlink.net http://www.geometer.org/mathcircles August 2, 2017 Abstract In this article we will examine various properties of iterated functions. If f(x) is a function,
More informationUNIVERSITETET I OSLO
UNIVERSITETET I OSLO Det matematisk-naturvitenskapelige fakultet Examination in: INF1100 Introduction to programming with scientific applications Day of examination: Thursday, October 13, 2011 Examination
More informationNumerical Methods Lecture 2 Simultaneous Equations
Numerical Methods Lecture 2 Simultaneous Equations Topics: matrix operations solving systems of equations pages 58-62 are a repeat of matrix notes. New material begins on page 63. Matrix operations: Mathcad
More informationLine Search Algorithms
Lab 1 Line Search Algorithms Investigate various Line-Search algorithms for numerical opti- Lab Objective: mization. Overview of Line Search Algorithms Imagine you are out hiking on a mountain, and you
More informationNumerical Integration
Chapter 1 Numerical Integration In this chapter we examine a few basic numerical techniques to approximate a definite integral. You may recall some of this from Calculus I where we discussed the left,
More informationASTRO 114 Lecture Okay. We re going to continue now with our discussion of stars. We re moving
ASTRO 114 Lecture 40 1 Okay. We re going to continue now with our discussion of stars. We re moving into the next chapter and the next chapter not only talks about the properties of stars but the entire
More informationNUMERICAL METHODS. x n+1 = 2x n x 2 n. In particular: which of them gives faster convergence, and why? [Work to four decimal places.
NUMERICAL METHODS 1. Rearranging the equation x 3 =.5 gives the iterative formula x n+1 = g(x n ), where g(x) = (2x 2 ) 1. (a) Starting with x = 1, compute the x n up to n = 6, and describe what is happening.
More informationMAT335H1F Lec0101 Burbulla
Fall 2011 Q 2 (x) = x 2 2 Q 2 has two repelling fixed points, p = 1 and p + = 2. Moreover, if I = [ p +, p + ] = [ 2, 2], it is easy to check that p I and Q 2 : I I. So for any seed x 0 I, the orbit of
More informationHow to work out really complicated motion. Iteration and Problem Solving Strategies. Let s go. Vertical spring-mass.
Iteration and Problem Solving Strategies How to solve anything! How to work out really complicated motion Break it up into little tiny steps. Use an approximate method for each step. Add them all up. Vertical
More informationMore Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n.
More Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n. If there are points which, after many iterations of map then fixed point called an attractor. fixed point, If λ
More information() Chapter 8 November 9, / 1
Example 1: An easy area problem Find the area of the region in the xy-plane bounded above by the graph of f(x) = 2, below by the x-axis, on the left by the line x = 1 and on the right by the line x = 5.
More informationLesson 3-1: Solving Linear Systems by Graphing
For the past several weeks we ve been working with linear equations. We ve learned how to graph them and the three main forms they can take. Today we re going to begin considering what happens when we
More informationDynamical Systems: Lecture 1 Naima Hammoud
Dynamical Systems: Lecture 1 Naima Hammoud Feb 21, 2017 What is dynamics? Dynamics is the study of systems that evolve in time What is dynamics? Dynamics is the study of systems that evolve in time a system
More informationMath 122 Fall Handout 15: Review Problems for the Cumulative Final Exam
Math 122 Fall 2008 Handout 15: Review Problems for the Cumulative Final Exam The topics that will be covered on Final Exam are as follows. Integration formulas. U-substitution. Integration by parts. Integration
More informationCalculus Favorite: Stirling s Approximation, Approximately
Calculus Favorite: Stirling s Approximation, Approximately Robert Sachs Department of Mathematical Sciences George Mason University Fairfax, Virginia 22030 rsachs@gmu.edu August 6, 2011 Introduction Stirling
More informationONE DIMENSIONAL CHAOTIC DYNAMICAL SYSTEMS
Journal of Pure and Applied Mathematics: Advances and Applications Volume 0 Number 0 Pages 69-0 ONE DIMENSIONAL CHAOTIC DYNAMICAL SYSTEMS HENA RANI BISWAS Department of Mathematics University of Barisal
More informationY1 Double Maths Assignment λ (lambda) Exam Paper to do Core 1 Solomon C on the VLE. Drill
α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω Nature is an infinite sphere of which the centre is everywhere and the circumference nowhere Blaise Pascal Y Double Maths Assignment λ (lambda) Tracking
More informationIntroduction to Nonlinear Dynamics and Chaos
Introduction to Nonlinear Dynamics and Chaos Sean Carney Department of Mathematics University of Texas at Austin Sean Carney (University of Texas at Austin) Introduction to Nonlinear Dynamics and Chaos
More informationMITOCW MITRES_18-007_Part5_lec3_300k.mp4
MITOCW MITRES_18-007_Part5_lec3_300k.mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high-quality educational resources
More informationEngg. Math. II (Unit-IV) Numerical Analysis
Dr. Satish Shukla of 33 Engg. Math. II (Unit-IV) Numerical Analysis Syllabus. Interpolation and Curve Fitting: Introduction to Interpolation; Calculus of Finite Differences; Finite Difference and Divided
More informationCalculus with Analytic Geometry I Exam 8 Take Home Part.
Calculus with Analytic Geometry I Exam 8 Take Home Part. INSTRUCTIONS: SHOW ALL WORK. Write clearly, using full sentences. Use equal signs appropriately; don t use them between quantities that are not
More informationLIMITS AT INFINITY MR. VELAZQUEZ AP CALCULUS
LIMITS AT INFINITY MR. VELAZQUEZ AP CALCULUS RECALL: VERTICAL ASYMPTOTES Remember that for a rational function, vertical asymptotes occur at values of x = a which have infinite its (either positive or
More informationWhat is a random variable
OKAN UNIVERSITY FACULTY OF ENGINEERING AND ARCHITECTURE MATH 256 Probability and Random Processes 04 Random Variables Fall 20 Yrd. Doç. Dr. Didem Kivanc Tureli didemk@ieee.org didem.kivanc@okan.edu.tr
More informationLesson 21 Not So Dramatic Quadratics
STUDENT MANUAL ALGEBRA II / LESSON 21 Lesson 21 Not So Dramatic Quadratics Quadratic equations are probably one of the most popular types of equations that you ll see in algebra. A quadratic equation has
More informationMath 165 Final Exam worksheet solutions
C Roettger, Fall 17 Math 165 Final Exam worksheet solutions Problem 1 Use the Fundamental Theorem of Calculus to compute f(4), where x f(t) dt = x cos(πx). Solution. From the FTC, the derivative of the
More informationand in each case give the range of values of x for which the expansion is valid.
α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω Mathematics is indeed dangerous in that it absorbs students to such a degree that it dulls their senses to everything else P Kraft Further Maths A (MFPD)
More informationConnectedness. Proposition 2.2. The following are equivalent for a topological space (X, T ).
Connectedness 1 Motivation Connectedness is the sort of topological property that students love. Its definition is intuitive and easy to understand, and it is a powerful tool in proofs of well-known results.
More informationProject 3: Pendulum. Physics 2300 Spring 2018 Lab partner
Physics 2300 Spring 2018 Name Lab partner Project 3: Pendulum In this project you will explore the behavior of a pendulum. There is no better example of a system that seems simple at first but turns out
More informationSolving Quadratic & Higher Degree Equations
Chapter 9 Solving Quadratic & Higher Degree Equations Sec 1. Zero Product Property Back in the third grade students were taught when they multiplied a number by zero, the product would be zero. In algebra,
More informationPhysics 351 Wednesday, January 10, 2018
Physics 351 Wednesday, January 10, 2018 Chapers 1 5 mostly review freshman physics, so we ll go through them very quickly in the first few days of class. Read Chapters 1+2 for Friday. Read Chapter 3 (momentum
More informationPractice Final Exam Solutions
Important Notice: To prepare for the final exam, one should study the past exams and practice midterms (and homeworks, quizzes, and worksheets), not just this practice final. A topic not being on the practice
More informationPhysics 509: Bootstrap and Robust Parameter Estimation
Physics 509: Bootstrap and Robust Parameter Estimation Scott Oser Lecture #20 Physics 509 1 Nonparametric parameter estimation Question: what error estimate should you assign to the slope and intercept
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationLecture 10: Powers of Matrices, Difference Equations
Lecture 10: Powers of Matrices, Difference Equations Difference Equations A difference equation, also sometimes called a recurrence equation is an equation that defines a sequence recursively, i.e. each
More informationMidterm Review. Igor Yanovsky (Math 151A TA)
Midterm Review Igor Yanovsky (Math 5A TA) Root-Finding Methods Rootfinding methods are designed to find a zero of a function f, that is, to find a value of x such that f(x) =0 Bisection Method To apply
More information