3.1: 1, 3, 5, 9, 10, 12, 14, 18
|
|
- Gerald Heath
- 5 years ago
- Views:
Transcription
1 3.:, 3, 5, 9,,, 4, 8 ) We want to solve d d c() d = f() with c() = c = constant and f() = for different boundary conditions to get w() and u(). dw d = dw d d = ( )d w() w() = w() = w() ( ) c d d = u() u() = w() c 6c (a) Boundary conditions: u() = and w() =. u() u() = w() c 6c u() = 6c w() ( ) ] d ( ) 3 + ] ( ) 3 + ] w() = w() ( ) ( ) 3 ] w() = ( ) (b) Boundary conditions: u() = and u() =. u() u() = w() ( ) 3 + ] c 6c u() = = w() c 6c w() = 6 w() = w() ( ) u() = 6c ( ) 3 ( ) ] w() = 6 ( ) 3) Looking at the case of a bar free at both ends, dw d = f() with w() = w() =, we want to know under which conditions a solution eists. Integrating both sides of the differential equation from to gives dw d d = w() + w() = f()d f()d So, for there to be a solution to this problem, f() must satisfy the condition f()d =. Physically, this means that the net force on the bar must be zero. 5) Again we have dw d = f() and c() d = w, with boundary conditions u() = and w() =. This time, f() = f = constant and c() jumps from c = for to c = for >. So, c() = + S( ). dw d d = fd w() + w() = ( )f w() = ( )f /7/ Page of 9
2 d = w() c() = f d = + S( ) f ( )d for ( )d for < f ( )d + f f u() u() = ] for 3f 8 + f 3 8 for < u() = f ] for f ] for < 9) The solution to cu +u = is u() = d +d ep( c )+. With the boundary conditions u() = u() =, d and d can be found as follows. u() = d + d = d = d u() = d + d ep( c ) = d ( ep( c ) ) = d = d = ep( c ) Now, we can find the limit of d, d, and u() as c from the positive side. lim d = lim d = lim c + c + c + ep( c ) = lim ep( c ) c + ep( c ) = = ep( ( ) c ) c + lim u() = + lim c + d + lim d ep( c + c + c ) = + lim ep( c ) = for = lim u() = c + for + for = + for When is away from = the limit of u() as c + clearly goes to. However, near = the limit rapidly diverges from and heads to. And if much faster than c, then the limit of u() will go to. The boundary condition u() = is always satisfied but the end = has a boundary layer. ) We want to solve u = with u() = u() = using finite elements with three hat functions (φ i ()) and h = 4. So, c() = and f() =, and we need to calculate the components of the matri K and vector F in the equation KU = F where U = (U, U, U 3 ) and K ij = c() dφ j d dv i d F i = f()φ i ()d U() = U φ () + U φ () + U 3 φ 3 () The test functions v i () are taken to be the trial functions φ i (), or v i () = φ i (). This makes the K matri symmetric. The components of K and F are given below. F = F = F 3 = h = ( ) K = K = K 33 = h h = 8 K 3 = K 3 = K = K = h ( h ) = 4 K 3 = K 3 = h ( h ) = 4 Now, we can solve the matri equation KU = F to get the vector U and the finite element solution U(). 8 4 U KU = F U = U U = U 3 U 3 3 /7/ Page of 9
3 Finally, we can compare to the finite element solution to the eact solution u() =. Clearly U() cannot equal u() since u() is quadratic and U() is piecewise linear. However, we can compare the value at the mesh points: = 4, =, and = 3 4. It is shown below that the values of U() and u() are identical at the meshpoints. ( ) u(h) = u = 3 ( ) 4 6 = U u(h) = u = 4 ( ) 6 = U 3 u(3h) = u = = U 3 ) The mass matri M is defined as M ij = V i V j d. We would like to know the mass matri for the three hat functions. So V i () would be a hat function centered about = hi, where h is the spacing between mesh points. It is clear that M is symmetric, so M ij = M ji. Also, the integral of the proct of two hat functions is non-zero only if the two hat functions are the same or are adjacent; in other words, M ij = if i j >. So, M is a tri-diagonal, symmetric matri when hat functions are used. Furthermore, since the integral of the square of a hat function is identical regardless of the which hat function is used, the values along the main diagonal are all equal to the same value (calculated to be h 3 ). Also, since the integral of the proct of two adjacent hat functions is identical regardless of which two adjacent pairs of hat functions is selected, the values along the diagonal adjacent to the main diagonal are all equal to the same value (calculated to be h 6 ). For the specific case of three hat functions dividing the interval to with h = 4, we get the following 3 3 matri for M: M = ) We can use Simpson s 6, 4 6, 6 rule to integrate and find the area underneath the bubble function φ 5(). h h φ 5 ()d = h 6 φ 5(h) + 4h 6 φ 5( 3h ) + h 6 φ 5(h) = h 6 () + 4h 6 () + h 6 () = h 3 8) To find the finite element solution to a fied-free hanging bar with c() = (boundary conditions u() = and u () = ), we must add an additional half-hat function, φ N+ () at the end in addition to the N interior hat functions φ i (). The spacing between meshpoints is h = N+ and the test functions are taken to be equal to the trial functions (V i () = φ i ()). (a) The stiffness matri K would now be (N + ) (N + ) instead of N N as in the fied-fied case. The etra row and column mostly have zero entries ecept near the lower right corner. The entry in row N, column (N + ) (which is identical to the entry in row (N + ), column N because K is symmetric) would be equal to h. The entry in row (N + ), column (N + ) would be equal to h. So multiplying the column vector U with the (N + ) row of K gives the epression U N+ U N h which is a first difference of U; specifically, it is a backward difference of U at the point =, which is analogous to u (). This last row of the stiffness matri is representing the boundary condition u () =. (b) With a constant load f() = f, the new last component of the column vector F would be F N+ = hf. The 3 3 matri K and 3 column vector F from problem can be generalized for the N N case. But in addition to that, we must add the etra row and column to get the larger (N + ) (N + ) stiffness matri K and the (N +) column vector F for this fied-free case. The form for K and F (for a fied-free hanging bar with c() = and f() = f ) are given by K = h... F = hf. /7/ Page 3 of 9
4 Now, we can solve KU = F to get the column vector U and the finite element solution U() = Σ (N+) i= U i φ i (), and then compare it with the true solution of u() = f ( ). Particularly, we would like to compare the values evaluated at the mesh points: for the finite element solution this is simply given by U i ; for the true solution it is given by u(hi) = h f ( h i)(i) = h f (N + ) i](i). U = K F U.. U N U N = h f N N... N (N + ) Through some careful counting we can see that for i =,..., (N + ): U i = h f (N + i)i ] i + Σi j=j = h f (N + )i i ] (i)(i + ) i + = h f (N + )i i i + i + i ] U i = h f (N + ) i] (i) This equation is identical to the equation for u(hi), so u(hi) = U i. 3.:, 3, 5, 6. ) We want to solve the equation for a cantilevered beam with a unit force at the midpoint, u () = δ( ) with boundary conditions u() = u () = and M() = M () =. The differential equation can be split into M () = δ( ) and u () = M(), and then solved as follows. M ()d = δ( )d M () M () = S( ) M ()d = (S( ) d M() M() = ( ) (R( ) ) u ()d = R( ) + + ] d u () u () = Q( ) + + u ()d = Q( )d + ] ( + ) d u() u() = C( ) + ( ) u() = 6 (3 + 3 ) for 6 (3 + 3 ) + ( )3 for < /7/ Page 4 of 9
5 3) We want to solve u () = δ() with boundary conditions u( ) = u ( ) = and u() = u () =. The general solution to u () = δ() is given by A + B + C + D 3 u() = 6 3 for A + B + C + D 3 for We can use the boundary conditions to get four equations that can be used to solve for the four unknowns A, B, C, and D. u( ) = A B + C D + 6 = u() = A + B + C + D = u ( ) = B C + 3D = u () = B + C + 3D = A A B 3 C = 6 B C = 3 D D u() = for for 5) Which of the 8 cubic finite elements, φ d (), φ s (),..., φ d 3(), φ s 3() based at the meshpoints =, 3, 3,, are dropped because of the essential boundary conditions below for u = f? (a) Fied-fied: u() = u() =. φ d () and φ d 3() would be dropped. (b) Fied-free: u() = u () =. Only u() = is an essential boundary condition. Thus, φ d () would be dropped. 6) Which of the 8 cubic finite elements, φ d (), φ s (),..., φ d 3(), φ s 3() based at the meshpoints =, 3, 3,, are dropped because of the essential boundary conditions below for u = f? (a) Built-in beam: u = u = at both ends. φ d (), φ s (), φ d 3(), and φ s 3() would be dropped. (b) Simply supported beam: u = u = at both ends. Only u() = u() = are the essential boundary conditions. Thus, φ d () and φ d 3() would be dropped. (c) Cantilevered beam: u() = u () = u () = u () =. Only u() = u () = are the essential boundary conditions. φ d () and φ s () would be dropped. 4 8 /7/ Page 5 of 9
6 MATLAB Assignment We want to convert between a geocentric rectangular coordinate system and a geodetic coordinate system. This conversion is useful in GPS applications. The conversion from geodetic coordinates (latitude φ, longitude λ, altitude h) to geocentric rectangular coordinates (X, Y, Z) is straightforward and is given by X = (N + h) cos φ cos λ () Y = (N + h) cos φ sin λ () Z = ( f) N + h ] sin φ (3) where N = a f( f) sin φ ] (4) and a = m is the length of the semi-major ais of the Earth and f = 97 is the polar flattening of the Earth. The conversion from geocentric rectangular coordinates to geodetic coordinates is more difficult because it requires solving a system of non-linear equations. The number of equations to solve can be reced from three to two by solving for λ analytically. By substituting equations and into the quotient Y X we get tan λ = Y X (5) By squaring both sides of equations and, adding the two resulting equations, and finally taking the square root of both sides of the summed equation, we get X + Y = (N + h) cos φ cos λ + sin λ = (N + h) cos φ The above equation is used to define the function g (φ, h), and equation 3 is used to define the function g (φ, h): g (φ, h) (N + h) cos φ X + Y = (6) g (φ, h) ( f) N + h ] sin φ Z = (7) where N is given in equation 4. Now, we need to calculate the Jacobian, g J = φ g φ g h g h The partial derivatives with respect to altitude h are easy to calculate. They are simply g h = cos φ and N = sin φ. To simplify calculating the partial derivatives with respect to φ, we will first calculate g h N φ = a f( f) sin φ ] 3 ( f( f) sin φ cos φ) = N 3 f( f) sin φ cos φ a φ. g φ = N 3 N cos φ (N + h) sin φ = φ a f( f) sin φ cos φ N sin φ h sin φ = N 3 ] f( a sin φ f) cos φ a N h sin φ = N 3 a = N 3 a sin φ f( f) cos φ + f( f) sin φ ] h sin φ sin φ + f f ] N 3 h sin φ = a ( f) + h sin φ /7/ Page 6 of 9
7 g φ N = ( f) φ sin φ + ( f) N + h ] cos φ N = ( f) 3 a f( f) sin φ cos φ + ( f) N cos φ + h cos φ = ( f) N 3 ] f( cos φ f) sin φ + a + h cos φ a = ( f) N 3 a cos φ f( f) sin φ + f( f) sin φ ] + h cos φ = ( f) N 3 N 3 cos φ + h cos φ = a a ( f) + h cos φ If we define D N 3 a ( f) + h, then the Jacobian can be simply written as D sin φ cos φ J = D cos φ sin φ and the inverse of the Jacobian is J sin φ cos φ sin φ cos φ = = det(j) D cos φ D sin φ D sin φ D cos φ D cos φ D sin φ = sin φ cos φ = D sin φ D cos φ ] D D cos φ D sin φ cos φ sin φ Now, we can find the solution vector u = (φ, h) using Newton s Method. We start with an initial guess vector u init and set the current vector u cur equal to it. For each iteration in the loop, we compute a new vector u new through the computation u new = u cur J g or φnew h new ] = φcur h cur N ] D sin φ cur cos φ cur D cos φ ] cur g (φ cur, h cur ) sin φ cur g (φ cur, h cur ) (8) where D = N 3 a ( f) + h (9) and the functions g and g are defined in equations 6 and 7 respectively, ecept the functions would of course not equal zero when evaluated at the current u cur unless u cur happened to eactly be the true solution. Every time the loop is repeated, u new of the prior loop becomes u cur of the current loop. The new calculated vector u new should be closer to the actual solution. Eventually, u new should be close enough to the actual solution, within some tolerance, that we can break out of the loop and treat u new as the solution. The MATLAB function that does this conversion essentially uses equations to 4 to convert the current estimate for the geodetic coordinates into geocentric rectangular coordinates (X cur, Y cur, Z cur ), and then it calculates the Euclidean distance between this point and the point (X, Y, Z) specified by the user. If this distance is smaller than some threshold distance, which can be specified by the user through an optional argument in the MATLAB function, it will break out of the loop. Newton s method requires an initial guess to start from; however, testing showed that regardless of where the initial guess was, the unique correct solution was quickly obtained in a few iterations. Therefore, the MATLAB function makes specifying the initial guess of the latitude and altitude optional, and by default it sets both to zero (at the equator with zero altitude). Finally, the longitude is eplicitly calculated from X and Y alone, using equation 5. The code for the MATLAB function, gps, that converts from geocentric rectangular to geodetic coordinates is displayed in Listing. /7/ Page 7 of 9
8 We can look at an eample where we wish to convert the geocentric rectangular coordinates (X, Y, Z) to geodetic coordinates latitude, longitude, and altitude. X, Y, and Z are specified as follows: X = Y = Z = m m m Then the following MATLAB code can be used to call the function and get the converted coordinates: phi, lambda, h] = gps( , , ) The variable h has the altitude (h) in meters, the variable phi has the latitude (φ) in degrees, and the variable lambda has the longitude (λ) in degrees. The values for these three variables, for this eample, are φ 57.3 λ 9.95 h m We can verify the results of the function by plugging these values into equations to 4 and getting the original specified values for X, Y, and Z. Listing : MATLAB function that converts from geocentral rectangular to geodetic coordinates. function l a t, long, a l t ] = gps ( X, Y, Z, v a r a r g i n ) %GPS Converts from g e o c e n t r i c r e c t a n g u l a r to g e o d e t i c c o o r d i n a t e s. 3 % l a t, long, a l t ] = gps (X, Y, Z, i n i t l a t =, i n i t a l t =,... 4 % t o l = e 6, ma iter = ) 5 % c o n v e r t s from g e o c e n t r i c r e c t a n g u l a r c o o r d i n a t e s (X, Y, Z) in meters to 6 % g e o d e t i c c o o r d i n a t e s : l a t i t u d e ( l a t ) and l o n g i t u d e ( long ) in d e g r e e s 7 % and a l t i t u d e ( a l t ) in meters. The conversion i s done by s o l v i n g a 8 % system o f non l i n e a r e q u a t i o n s using Newton s method, so an i n i t i a l 9 % guess f o r t he l a t i t u d e ( i n i t l a t ) in d e g r e e s and the a l t i t u d e % ( i n i t a l t ) in meters i s necessary. By d e f a u l t the i n i t i a l l a t i t u d e i s % s e t to d e g r e e s ( equator ) and the i n i t i a l a l t i t u d e i s s e t to m. % 3 % The t o l e r a n c e ( t o l ), which s p e c i f i e s the maimum a c c e p t i b l e d i s t a n c e in 4 % meters between the s p e c i f i e d (X, Y, Z) p o i n t and the e s t i m a t e d p o i n t 5 % g iven by the c a l c u l a t e d l a t i t u d e, l o n g i t u d e, a l t i t u d e, i s s e t to e 6 m 6 % by d e f a u l t. The maimum number o f i t e r a t i o n s ( ma iter ) b e f o r e a b o r t i n g 7 % Newton s method and p r i n t i n g an e r r o r message i s s e t to by d e f a u l t. 8 9 % Constants f o r Earth a = ; % m] l e n g t h o f semi major a i s o f Earth f = /97; % p o l a r f l a t t e n i n g 3 % Set d e f a u l t v a l u e s f o r t o l e r a n c e and ma iter i f not s p e c i f i e d. 4 numvarargs = length ( v a r a r g i n ) ; 5 i f ( numvarargs > 4) 6 error ( Requires at most 4 o p t i o n a l inputs ) ; 7 end 8 optargs =,, e 6, }; % Set d e f a u l t s 9 newvals = c e l l f u n (@( ) isempty ( ), v a r a r g i n ) ; 3 optargs ( newvals ) = v a r a r g i n ( newvals ) ; 3 i n i t l a t, i n i t a l t, t o l, ma iter ] = optargs : } ; 3 /7/ Page 8 of 9
9 33 % Newton s Method 34 u = degtorad ( i n i t l a t ) ; i n i t a l t ] ; 35 g r i g h t = sqrt (X X + Y Y) ; Z ] ; 36 g right norm = norm( g r i g h t ) ; % g r i g h t n o r m = s q r t (X X + Y Y + Z Z ) ; 37 for i =: ma iter ] 38 phi = u ( ) ; h = u ( ) ; 39 N = a ( f ( f ) ( sin ( phi )ˆ))ˆ( /); 4 D = h + ( f )ˆ Nˆ3 / a ˆ ; 4 g l e f t = (N+h ) cos ( phi ) ; (N ( f )ˆ + h ) sin ( phi ) ] ; 4 g = g l e f t g r i g h t ; 43 u = u sin ( phi )/D, cos ( phi )/D ; cos ( phi ), sin ( phi ) ] g ; % C a l c u l a t e c u r r e n t d i s t a n c e from d e s i r e d p o i n t ( e r r o r ) and check i f i t 46 % i s w i t h i n t he s p e c i f i e d t o l e r a n c e. 47 i f ( abs (norm( g l e f t ) g right norm ) < t o l ) 48 break ; % Tolerance s a t i s f i e d, break out o f loop. 49 end 5 end 5 5 % Check i f loop aborted because o f maimum i t e r a t i o n s 53 i f ( i >= ma iter ) 54 error ( gps : maiteration, Eceeded maimum number o f i t e r a t i o n s (%d ) without f i n d i n g, a s o l u t i o n within the t o l e r a n c e ], ma iter ) ; 57 end % Set output v a r i a b l e s 6 long = radtodeg ( atan (Y, X) ) ; 6 l a t = radtodeg ( phi ) ; 6 a l t = h ; 63 end /7/ Page 9 of 9
Numerical Methods. Root Finding
Numerical Methods Solving Non Linear 1-Dimensional Equations Root Finding Given a real valued function f of one variable (say ), the idea is to find an such that: f() 0 1 Root Finding Eamples Find real
More informationChapter 6. Nonlinear Equations. 6.1 The Problem of Nonlinear Root-finding. 6.2 Rate of Convergence
Chapter 6 Nonlinear Equations 6. The Problem of Nonlinear Root-finding In this module we consider the problem of using numerical techniques to find the roots of nonlinear equations, f () =. Initially we
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationNumerical Integration (Quadrature) Another application for our interpolation tools!
Numerical Integration (Quadrature) Another application for our interpolation tools! Integration: Area under a curve Curve = data or function Integrating data Finite number of data points spacing specified
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 informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
More informationChapter 9. Derivatives. Josef Leydold Mathematical Methods WS 2018/19 9 Derivatives 1 / 51. f x. (x 0, f (x 0 ))
Chapter 9 Derivatives Josef Leydold Mathematical Methods WS 208/9 9 Derivatives / 5 Difference Quotient Let f : R R be some function. The the ratio f = f ( 0 + ) f ( 0 ) = f ( 0) 0 is called difference
More informationDepartment of Applied Mathematics and Theoretical Physics. AMA 204 Numerical analysis. Exam Winter 2004
Department of Applied Mathematics and Theoretical Physics AMA 204 Numerical analysis Exam Winter 2004 The best six answers will be credited All questions carry equal marks Answer all parts of each question
More information18.303: Introduction to Green s functions and operator inverses
8.33: Introduction to Green s functions and operator inverses S. G. Johnson October 9, 2 Abstract In analogy with the inverse A of a matri A, we try to construct an analogous inverse  of differential
More informationIntroduction to Finite Element Method. Dr. Aamer Haque
Introduction to Finite Element Method 4 th Order Beam Equation Dr. Aamer Haque http://math.iit.edu/~ahaque6 ahaque7@iit.edu Illinois Institute of Technology July 1, 009 Outline Euler-Bernoulli Beams Assumptions
More informationLinear Equations in Linear Algebra
1 Linear Equations in Linear Algebra 1.1 SYSTEMS OF LINEAR EQUATIONS LINEAR EQUATION,, 1 n A linear equation in the variables equation that can be written in the form a a a b 1 1 2 2 n n a a is an where
More informationNUMERICAL METHODS FOR SOLVING EQUATIONS
Mathematics Revision Guides Numerical Methods for Solving Equations Page of M.K. HOME TUITION Mathematics Revision Guides Level: AS / A Level AQA : C3 Edecel: C3 OCR: C3 NUMERICAL METHODS FOR SOLVING EQUATIONS
More information3.3 Limits and Infinity
Calculus Maimus. Limits Infinity Infinity is not a concrete number, but an abstract idea. It s not a destination, but a really long, never-ending journey. It s one of those mind-warping ideas that is difficult
More informationTopic 8c Multi Variable Optimization
Course Instructor Dr. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: rcrumpf@utep.edu Topic 8c Multi Variable Optimization EE 4386/5301 Computational Methods in EE Outline Mathematical Preliminaries
More informationIntegration. Section 8: Using partial fractions in integration
Integration Section 8: Using partial fractions in integration Notes and Eamples These notes contain subsections on Using partial fractions in integration Putting all the integration techniques together
More informationMAE 107 Homework 8 Solutions
MAE 107 Homework 8 Solutions 1. Newton s method to solve 3exp( x) = 2x starting at x 0 = 11. With chosen f(x), indicate x n, f(x n ), and f (x n ) at each step stopping at the first n such that f(x n )
More information= V I = Bus Admittance Matrix. Chapter 6: Power Flow. Constructing Ybus. Example. Network Solution. Triangular factorization. Let
Chapter 6: Power Flow Network Matrices Network Solutions Newton-Raphson Method Fast Decoupled Method Bus Admittance Matri Let I = vector of currents injected into nodes V = vector of node voltages Y bus
More information2 3 x = 6 4. (x 1) 6
Solutions to Math 201 Final Eam from spring 2007 p. 1 of 16 (some of these problem solutions are out of order, in the interest of saving paper) 1. given equation: 1 2 ( 1) 1 3 = 4 both sides 6: 6 1 1 (
More information3 Applications of Derivatives Instantaneous Rates of Change Optimization Related Rates... 13
Contents Limits Derivatives 3. Difference Quotients......................................... 3. Average Rate of Change...................................... 4.3 Derivative Rules...........................................
More informationReview sheet Final Exam Math 140 Calculus I Fall 2015 UMass Boston
Review sheet Final Eam Math Calculus I Fall 5 UMass Boston The eam is closed tetbook NO CALCULATORS OR ELECTRONIC DEVICES ARE ALLOWED DURING THE EXAM The final eam will contain problems of types similar
More informationIn the real world, objects don t just move back and forth in 1-D! Projectile
Phys 1110, 3-1 CH. 3: Vectors In the real world, objects don t just move back and forth in 1-D In principle, the world is really 3-dimensional (3-D), but in practice, lots of realistic motion is 2-D (like
More informationHomework 5 - Solutions
Spring Math 54 Homework 5 - Solutions BF 3.4.4. d. The spline interpolation routine below produces the following coefficients: i a i b i c i d i -..869948.75637848.656598 -.5.9589.487644.9847639.887.9863.34456976.489747
More informationACCUPLACER MATH 0310
The University of Teas at El Paso Tutoring and Learning Center ACCUPLACER MATH 00 http://www.academics.utep.edu/tlc MATH 00 Page Linear Equations Linear Equations Eercises 5 Linear Equations Answer to
More informationStatistical Geometry Processing Winter Semester 2011/2012
Statistical Geometry Processing Winter Semester 2011/2012 Linear Algebra, Function Spaces & Inverse Problems Vector and Function Spaces 3 Vectors vectors are arrows in space classically: 2 or 3 dim. Euclidian
More informationUniversity of Alabama Department of Physics and Astronomy. PH 125 / LeClair Spring A Short Math Guide. Cartesian (x, y) Polar (r, θ)
University of Alabama Department of Physics and Astronomy PH 125 / LeClair Spring 2009 A Short Math Guide 1 Definition of coordinates Relationship between 2D cartesian (, y) and polar (r, θ) coordinates.
More informationMathematics Numbers: Absolute Value of Functions I
a place of mind F A C U L T Y O F E D U C A T I O N Department of Curriculum and Pedagogy Mathematics Numbers: Absolute Value of Functions I Science and Mathematics Education Research Group Supported by
More informationChapter 12: Iterative Methods
ES 40: Scientific and Engineering Computation. Uchechukwu Ofoegbu Temple University Chapter : Iterative Methods ES 40: Scientific and Engineering Computation. Gauss-Seidel Method The Gauss-Seidel method
More informationFixed Point Theorem and Sequences in One or Two Dimensions
Fied Point Theorem and Sequences in One or Two Dimensions Dr. Wei-Chi Yang Let us consider a recursive sequence of n+ = n + sin n and the initial value can be an real number. Then we would like to ask
More informationLimits and Continuity
Limits and Continuity Philippe B. Laval Kennesaw State University January 2, 2005 Contents Abstract Notes and practice problems on its and continuity. Limits 2. Introduction... 2.2 Theory:... 2.2. GraphicalMethod...
More informationMath 471 (Numerical methods) Chapter 3 (second half). System of equations
Math 47 (Numerical methods) Chapter 3 (second half). System of equations Overlap 3.5 3.8 of Bradie 3.5 LU factorization w/o pivoting. Motivation: ( ) A I Gaussian Elimination (U L ) where U is upper triangular
More informationCLASS NOTES Computational Methods for Engineering Applications I Spring 2015
CLASS NOTES Computational Methods for Engineering Applications I Spring 2015 Petros Koumoutsakos Gerardo Tauriello (Last update: July 2, 2015) IMPORTANT DISCLAIMERS 1. REFERENCES: Much of the material
More information0.1. Linear transformations
Suggestions for midterm review #3 The repetitoria are usually not complete; I am merely bringing up the points that many people didn t now on the recitations Linear transformations The following mostly
More information(i) find the points where f(x) is discontinuous, and classify each point of discontinuity.
Math Final Eam - Practice Problems. A function f is graphed below. f() 5 4 8 7 5 4 4 5 7 8 4 5 (a) Find f(0), f( ), f(), and f(4) Find the domain and range of f (c) Find the intervals where f () is positive
More information3.3.1 Linear functions yet again and dot product In 2D, a homogenous linear scalar function takes the general form:
3.3 Gradient Vector and Jacobian Matri 3 3.3 Gradient Vector and Jacobian Matri Overview: Differentiable functions have a local linear approimation. Near a given point, local changes are determined by
More informationCOUNCIL ROCK HIGH SCHOOL MATHEMATICS. A Note Guideline of Algebraic Concepts. Designed to assist students in A Summer Review of Algebra
COUNCIL ROCK HIGH SCHOOL MATHEMATICS A Note Guideline of Algebraic Concepts Designed to assist students in A Summer Review of Algebra [A teacher prepared compilation of the 7 Algebraic concepts deemed
More informationHP-33S Calculator Program XYZ 2
Convert XYZ Geocentric Co-ordinates on any ellipsoid to Latitude, Longitude and Height Programmer: Dr. Bill Hazelton Date: May, 2005. Version: 2.0 Line Instruction Display User Instructions Y0001 LBL Y
More informationMA3232 Numerical Analysis Week 9. James Cooley (1926-)
MA umerical Analysis Week 9 James Cooley (96-) James Cooley is an American mathematician. His most significant contribution to the world of mathematics and digital signal processing is the Fast Fourier
More informationBrief Revision Notes and Strategies
Brief Revision Notes and Strategies Straight Line Distance Formula d = ( ) + ( y y ) d is distance between A(, y ) and B(, y ) Mid-point formula +, y + M y M is midpoint of A(, y ) and B(, y ) y y Equation
More informationFinite Elements for Nonlinear Problems
Finite Elements for Nonlinear Problems Computer Lab 2 In this computer lab we apply finite element method to nonlinear model problems and study two of the most common techniques for solving the resulting
More informationNonlinear examples. June 21, Solution by fixed-point iteration 3. 5 Exercises 9
Nonlinear examples June 21, 2012 Contents 1 Bratu s equation 1 1.1 Manufacturing a problem to match a solution.................. 2 1.1.1 A manufactured Bratu s problem in 1D.................. 2 1.2 Weak
More informationPage 1. These are all fairly simple functions in that wherever the variable appears it is by itself. What about functions like the following, ( ) ( )
Chain Rule Page We ve taken a lot of derivatives over the course of the last few sections. However, if you look back they have all been functions similar to the following kinds of functions. 0 w ( ( tan
More informationPolynomial Functions and Their Graphs
Polynomial Functions and Their Graphs Definition of a Polynomial Function Let n be a nonnegative integer and let a n, a n- 1,, a 2, a 1, a 0, be real numbers with a n 0. The function defined by f (x) a
More informationIteration & Fixed Point
Iteration & Fied Point As a method for finding the root of f this method is difficult, but it illustrates some important features of iterstion. We could write f as f g and solve g. Definition.1 (Fied Point)
More informationGeometric Modeling Summer Semester 2010 Mathematical Tools (1)
Geometric Modeling Summer Semester 2010 Mathematical Tools (1) Recap: Linear Algebra Today... Topics: Mathematical Background Linear algebra Analysis & differential geometry Numerical techniques Geometric
More informationMath 261 Final Exam - Practice Problem Solutions. 1. A function f is graphed below.
Math Final Eam - Practice Problem Solutions. A function f is graphed below. f() 8 7 7 8 (a) Find f(), f( ), f(), and f() f() = ;f( ).;f() is undefined; f() = (b) Find the domain and range of f Domain:
More informationNumerical Analysis Fall. Roots: Open Methods
Numerical Analysis 2015 Fall Roots: Open Methods Open Methods Open methods differ from bracketing methods, in that they require only a single starting value or two starting values that do not necessarily
More informationFinding Roots of Equations
Finding Roots of Equations Solution Methods Overview Bisection/Half-interval Search Method of false position/regula Falsi Secant Method Newton Raphson Iteration Method Many more. Open Methods Bracketing
More information2.3 The Fixed-Point Algorithm
.3 The Fied-Point Algorithm 1. Mean Value Theorem: Theorem Rolle stheorem: Suppose that f is continuous on a, b and is differentiable on a, b. If f a f b, then there eists a number c in a, b such that
More informationC. Finding roots of trinomials: 1st Example: x 2 5x = 14 x 2 5x 14 = 0 (x 7)(x + 2) = 0 Answer: x = 7 or x = -2
AP Calculus Students: Welcome to AP Calculus. Class begins in approimately - months. In this packet, you will find numerous topics that were covered in your Algebra and Pre-Calculus courses. These are
More informationCHAPTER 3 Applications of Differentiation
CHAPTER Applications of Differentiation Section. Etrema on an Interval.............. 0 Section. Rolle s Theorem and the Mean Value Theorem. 07 Section. Increasing and Decreasing Functions and the First
More informationMEI Core 2. Sequences and series. Section 1: Definitions and Notation
Notes and Eamples MEI Core Sequences and series Section : Definitions and Notation In this section you will learn definitions and notation involving sequences and series, and some different ways in which
More informationlim x c) lim 7. Using the guidelines discussed in class (domain, intercepts, symmetry, asymptotes, and sign analysis to
Math 7 REVIEW Part I: Problems Using the precise definition of the it, show that [Find the that works for any arbitrarily chosen positive and show that it works] Determine the that will most likely work
More informationMat104 Fall 2002, Improper Integrals From Old Exams
Mat4 Fall 22, Improper Integrals From Old Eams For the following integrals, state whether they are convergent or divergent, and give your reasons. () (2) (3) (4) (5) converges. Break it up as 3 + 2 3 +
More information2.13 Linearization and Differentials
Linearization and Differentials Section Notes Page Sometimes we can approimate more complicated functions with simpler ones These would give us enough accuracy for certain situations and are easier to
More informationi nfoclearinghouse.com Solving Non-linear Equations with SCILAB Gilberto E. Urroz, Ph.D., P.E. Distributed by
Solving Non-linear Equations with SCILAB By Gilberto E. Urroz, Ph.D., P.E. Distributed by i nfoclearinghouse.com 00 Gilberto E. Urroz All Rights Reserved A "zip" file containing all of the programs in
More information1 Exponential Functions Limit Derivative Integral... 5
Contents Eponential Functions 3. Limit................................................. 3. Derivative.............................................. 4.3 Integral................................................
More informationTRIG REVIEW NOTES. Co-terminal Angles: Angles that end at the same spot. (sines, cosines, and tangents will equal)
TRIG REVIEW NOTES Convert from radians to degrees: multiply by 0 180 Convert from degrees to radians: multiply by 0. 180 Co-terminal Angles: Angles that end at the same spot. (sines, cosines, and tangents
More informationPHYSICS 1030 Homework #9
PHYSICS 1030 Homework #9 (Due Dec. 6, 2017) Find the position of the planet Mars at time t D December 6, 2017, 5:00 am EST. You will do this by following the steps shown below. (a) Convert the time t to
More informationExact and heuristic minimization of Boolean functions
Eact and heuristic minimization of Boolean functions Quine a McCluskey Method a) Generation of all prime implicants b) Selection of a minimum subset of prime implicants, which will represent the original
More informationdu(l) 5 dl = U(0) = 1 and 1.) Substitute for U an unspecified trial function into governing equation, i.e. dx + = =
Consider an ODE of te form: Finite Element Metod du fu g d + wit te following Boundary Conditions: U(0) and du(l) 5 dl.) Substitute for U an unspecified trial function into governing equation, i.e. ^ U
More informationSEE and DISCUSS the pictures on pages in your text. Key picture:
Math 6 Notes 1.1 A PREVIEW OF CALCULUS There are main problems in calculus: 1. Finding a tangent line to a curve though a point on the curve.. Finding the area under a curve on some interval. SEE and DISCUSS
More informationis the intuition: the derivative tells us the change in output y (from f(b)) in response to a change of input x at x = b.
Uses of differentials to estimate errors. Recall the derivative notation df d is the intuition: the derivative tells us the change in output y (from f(b)) in response to a change of input at = b. Eamples.
More informationComputer Problems for Taylor Series and Series Convergence
Computer Problems for Taylor Series and Series Convergence The two problems below are a set; the first should be done without a computer and the second is a computer-based follow up. 1. The drawing below
More informationLinear Algebra Review
Chapter 1 Linear Algebra Review It is assumed that you have had a course in linear algebra, and are familiar with matrix multiplication, eigenvectors, etc. I will review some of these terms here, but quite
More informationLinear Algebraic Equations
Linear Algebraic Equations Linear Equations: a + a + a + a +... + a = c 11 1 12 2 13 3 14 4 1n n 1 a + a + a + a +... + a = c 21 2 2 23 3 24 4 2n n 2 a + a + a + a +... + a = c 31 1 32 2 33 3 34 4 3n n
More information1. Introduction. Consider the Navier Stokes equations ηu t ν 2 u + (u grad) u + grad p = f div u = 0 (1.1)
University of Maryland Department of Computer Science TR-C4930 University of Maryland Institute for Advanced Computer Studies TR-009-0 BOUNDARY CONDITIONS IN APPROXIMATE COMMUTATOR PRECONDITIONERS FOR
More informationFunctions with orthogonal Hessian
Functions with orthogonal Hessian B. Dacorogna P. Marcellini E. Paolini Abstract A Dirichlet problem for orthogonal Hessians in two dimensions is eplicitly solved, by characterizing all piecewise C 2 functions
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 4
2.29 Spring 2015 Lecture 4 Review Lecture 3 Truncation Errors, Taylor Series and Error Analysis Taylor series: 2 3 n n i1 i i i i i n f( ) f( ) f '( ) f ''( ) f '''( )... f ( ) R 2! 3! n! n1 ( n1) Rn f
More informationPreface. Computing Definite Integrals. 3 cos( x) dx. x 3
Preface Here are the solutions to the practice problems for my Calculus I notes. Some solutions will have more or less detail than other solutions. The level of detail in each solution will depend up on
More informationSection 3.3 Graphs of Polynomial Functions
3.3 Graphs of Polynomial Functions 179 Section 3.3 Graphs of Polynomial Functions In the previous section we eplored the short run behavior of quadratics, a special case of polynomials. In this section
More informationPractice Midterm Solutions
Practice Midterm Solutions Math 4B: Ordinary Differential Equations Winter 20 University of California, Santa Barbara TA: Victoria Kala DO NOT LOOK AT THESE SOLUTIONS UNTIL YOU HAVE ATTEMPTED EVERY PROBLEM
More informationThe Fundamental Theorem of Calculus Part 3
The Fundamental Theorem of Calculus Part FTC Part Worksheet 5: Basic Rules, Initial Value Problems, Rewriting Integrands A. It s time to find anti-derivatives algebraically. Instead of saying the anti-derivative
More informationADDITIONAL MATHEMATICS
ADDITIONAL MATHEMATICS Paper 0606/ Paper Key messages Candidates should be reminded of the importance of reading the rubric on the eamination paper. Accuracy is of vital importance with final answers to
More informationA BRIEF REVIEW OF ALGEBRA AND TRIGONOMETRY
A BRIEF REVIEW OF ALGEBRA AND TRIGONOMETR Some Key Concepts:. The slope and the equation of a straight line. Functions and functional notation. The average rate of change of a function and the DIFFERENCE-
More informationSECTION 4-3 Approximating Real Zeros of Polynomials Polynomial and Rational Functions
Polynomial and Rational Functions 79. P() 9 9 8. P() 6 6 8 7 8 8. The solutions to the equation are all the cube roots of. (A) How many cube roots of are there? (B) is obviously a cube root of ; find all
More informationSYMBOLIC AND EXACT STRUCTURE PREDICTION FOR SPARSE GAUSSIAN ELIMINATION WITH PARTIAL PIVOTING
SYMBOLIC AND EXACT STRUCTURE PREDICTION FOR SPARSE GAUSSIAN ELIMINATION WITH PARTIAL PIVOTING LAURA GRIGORI, JOHN R. GILBERT, AND MICHEL COSNARD Abstract. In this paper we consider two structure prediction
More informationMOL 410/510: Introduction to Biological Dynamics Fall 2012 Problem Set #4, Nonlinear Dynamical Systems (due 10/19/2012) 6 MUST DO Questions, 1
MOL 410/510: Introduction to Biological Dynamics Fall 2012 Problem Set #4, Nonlinear Dynamical Systems (due 10/19/2012) 6 MUST DO Questions, 1 OPTIONAL question 1. Below, several phase portraits are shown.
More informationAP Calculus AB Summer Assignment
AP Calculus AB Summer Assignment Name: When you come back to school, you will be epected to have attempted every problem. These skills are all different tools that you will pull out of your toolbo this
More informationAPPM 1360 Final Exam Spring 2016
APPM 36 Final Eam Spring 6. 8 points) State whether each of the following quantities converge or diverge. Eplain your reasoning. a) The sequence a, a, a 3,... where a n ln8n) lnn + ) n!) b) ln d c) arctan
More informationMath Review. Daniel B. Rowe, Ph.D. Professor Department of Mathematics, Statistics, and Computer Science. Copyright 2018 by D.B.
Math Review Daniel B. Rowe, Ph.D. Professor Department of Mathematics, Statistics, and Computer Science Copyright 2018 by 1 Outline Differentiation Definition Analytic Approach Numerical Approach Integration
More informationA Basic Course in Real Analysis Prof. P. D. Srivastava Department of Mathematics Indian Institute of Technology, Kharagpur
A Basic Course in Real Analysis Prof. P. D. Srivastava Department of Mathematics Indian Institute of Technology, Kharagpur Lecture - 36 Application of MVT, Darbou Theorem, L Hospital Rule (Refer Slide
More informationEigenvectors and Eigenvalues 1
Ma 2015 page 1 Eigenvectors and Eigenvalues 1 In this handout, we will eplore eigenvectors and eigenvalues. We will begin with an eploration, then provide some direct eplanation and worked eamples, and
More informationChapter 4. Solution of a Single Nonlinear Algebraic Equation
Single Nonlinear Algebraic Equation - 56 Chapter 4. Solution of a Single Nonlinear Algebraic Equation 4.1. Introduction Life, my fris, is nonlinear. As such, in our roles as problem-solvers, we will be
More informationUBC-SFU-UVic-UNBC Calculus Exam Solutions 7 June 2007
This eamination has 15 pages including this cover. UBC-SFU-UVic-UNBC Calculus Eam Solutions 7 June 007 Name: School: Signature: Candidate Number: Rules and Instructions 1. Show all your work! Full marks
More informationContinuous functions. Limits of non-rational functions. Squeeze Theorem. Calculator issues. Applications of limits
Calculus Lia Vas Continuous functions. Limits of non-rational functions. Squeeze Theorem. Calculator issues. Applications of limits Continuous Functions. Recall that we referred to a function f() as a
More informationSolutions to Problem Sheet for Week 6
THE UNIVERSITY OF SYDNEY SCHOOL OF MATHEMATICS AND STATISTICS Solutions to Problem Sheet for Week 6 MATH90: Differential Calculus (Advanced) Semester, 07 Web Page: sydney.edu.au/science/maths/u/ug/jm/math90/
More informationPHYSICS 1030 Homework #9
PHYSICS 1030 Homework #9 (Due Dec. 5, 2018, 6:00 pm) Find the position of the planet Mars at time t D December 5, 2018, 7:50 pm EST. You will do this by following the steps shown below. (a) Convert the
More informationMECh300H Introduction to Finite Element Methods. Finite Element Analysis (F.E.A.) of 1-D Problems
MECh300H Introduction to Finite Element Methods Finite Element Analysis (F.E.A.) of -D Problems Historical Background Hrenikoff, 94 frame work method Courant, 943 piecewise polynomial interpolation Turner,
More informationLimits, Continuity, and Differentiability Solutions
Limits, Continuity, and Differentiability Solutions We have intentionally included more material than can be covered in most Student Study Sessions to account for groups that are able to answer the questions
More informationCHAPTER 3 Applications of Differentiation
CHAPTER Applications of Differentiation Section. Etrema on an Interval................... 0 Section. Rolle s Theorem and the Mean Value Theorem...... 0 Section. Increasing and Decreasing Functions and
More informationMIDTERM EXAM SOLUTIONS
MIDTERM EXAM 1 2012 SOLUTIONS Directions and rules. The eam will last 70 minutes; the last five minutes of class will be used for collecting the eams. No electronic devices of any kind will be allowed,
More informationSolutions to review problems MAT 125, Fall 2004
Solutions to review problems MAT 125, Fall 200 1. For each of the following functions, find the absolute maimum and minimum values for f() in the given intervals. Also state the value where they occur.
More informationThe Gradient. Consider the topography of the Earth s surface.
9/16/5 The Gradient.doc 1/8 The Gradient Consider the topography of the Earth s surface. We use contours of constant elevation called topographic contours to epress on maps (a -dimensional graphic) the
More informationMATH.2720 Introduction to Programming with MATLAB Vector and Matrix Algebra
MATH.2720 Introduction to Programming with MATLAB Vector and Matrix Algebra A. Vectors A vector is a quantity that has both magnitude and direction, like velocity. The location of a vector is irrelevant;
More informationUseful Mathematics. 1. Multivariable Calculus. 1.1 Taylor s Theorem. Monday, 13 May 2013
Useful Mathematics Monday, 13 May 013 Physics 111 In recent years I have observed a reticence among a subpopulation of students to dive into mathematics when the occasion arises in theoretical mechanics
More informationROOT FINDING REVIEW MICHELLE FENG
ROOT FINDING REVIEW MICHELLE FENG 1.1. Bisection Method. 1. Root Finding Methods (1) Very naive approach based on the Intermediate Value Theorem (2) You need to be looking in an interval with only one
More informationSolutions to the Exercises of Chapter 8
8A Domains of Functions Solutions to the Eercises of Chapter 8 1 For 7 to make sense, we need 7 0or7 So the domain of f() is{ 7} For + 5 to make sense, +5 0 So the domain of g() is{ 5} For h() to make
More information(x 3)(x + 5) = (x 3)(x 1) = x + 5
RMT 3 Calculus Test olutions February, 3. Answer: olution: Note that + 5 + 3. Answer: 3 3) + 5) = 3) ) = + 5. + 5 3 = 3 + 5 3 =. olution: We have that f) = b and f ) = ) + b = b + 8. etting these equal
More informationAP Calculus AB Summer Assignment
AP Calculus AB Summer Assignment Name: When you come back to school, it is my epectation that you will have this packet completed. You will be way behind at the beginning of the year if you haven t attempted
More informationAnswer Key 1973 BC 1969 BC 24. A 14. A 24. C 25. A 26. C 27. C 28. D 29. C 30. D 31. C 13. C 12. D 12. E 3. A 32. B 27. E 34. C 14. D 25. B 26.
Answer Key 969 BC 97 BC. C. E. B. D 5. E 6. B 7. D 8. C 9. D. A. B. E. C. D 5. B 6. B 7. B 8. E 9. C. A. B. E. D. C 5. A 6. C 7. C 8. D 9. C. D. C. B. A. D 5. A 6. B 7. D 8. A 9. D. E. D. B. E. E 5. E.
More information