Quadrature approaches to the solution of two point boundary value problems
|
|
- Moris Eaton
- 5 years ago
- Views:
Transcription
1 Quadrature approaches to the solution of two point boundary value problems Seth F. Oppenheimer Mohsen Razzaghi Department of Mathematics and Statistics Mississippi State University Drawer MA MSU, MS Abstract In a split level numerical analysis course or introductory applied mathematics course the topics of Gaussian Quadrature and the solution of two point boundary value problems are often covered. It seems to the students that these topics are unrelated. However, by using Green s functions, one can develop very nice easily implemented numerical schemes to approimate the solution to two point boundary value problems via Gaussian quadrature. 1 Introduction The purpose of this note is to present a technique for the numerical solution of two point boundary value problems which, while well known in the research community [4], is not taught at the undergraduate level. This is unfortunate as the technique is easy to implement and remarkably accurate. In a first numerical analysis course, Gaussian quadrature is usually introduced. However, beyond approimating the integrals of some known functions, little is done with it. We will demonstrate a very nice, nontrivial application for Gaussian quadrature to the numerical solution of differential equations. For completeness, we will outline all of the necessary background topics. In section 2, we review orthogonal polynomials and Gaussian quadrature. In section 3, we will state the general two point boundary value problem we wish to solve and give an eample of transforming such a problem into an equivalent Fredholm integral equation via the use of a Green s function. Finally, in the last section we will show how to use Gaussian quadrature to transform the integral representation of the differential equation into a system of algebraic equations that will yield an approimate solution of the problem. Two eamples will be given. 1
2 2 Orthogonal polynomials and Gaussian quadrature In this section we will develop the set of orthogonal polynomials on the interval [- 1,1] called the Legendre polynomials and their application to the approimation of definite integrals, Gaussian quadrature. Treatments of both of these subjects may be found in standard numerical analysis books such as [2] and [5]. An inepensive monograph on orthogonal polynomials is [6] and the more general application of the techniques we are discussing in this paper may be found in [3]. We consider the space of continuous functions on [, 1], C ([, 1]). We know from the Stone-Weierstrass theorem that the set of polynomials is dense in C ([, 1]), that is given a continuous function f and a positive number ε, there is a polynomial p so that for all in [, 1] p () f () < ε. The question is, given an f, howdowefind p? The answer is interpolation. That is, we choose N +1 points 0,..., N and find the unique polynomial of degree N +1 that agrees with f at those points. That is, we want the polynomial to take on the value f ( i ) at i. However, for some functions and some choices of points, increasing the number of interpolation points will not improve the polynomial approimation. To find sets of interpolation points so that as the number of interpolating points increases the difference between f and the interpolating polynomial decreases requires a detour into a discussion of orthogonal polynomials. Given a positive integrable function ρ on [, 1], called a weight function, we can define an inner product on C [(, 1)] by (f,g) = f()g()ρ()d This allows us to define orthogonal functions by saying that f and g are orthogonal if (f,g) =0. Since a polynomial of finite degree has only a finite set of roots, it is easy to see that any set 1,, 2,..., Mª is linearly independent in C [(, 1)]. We can therefore use the Gram-Schmidt orthogonalization process q 0 () = 1 q k () = k X k,q i () (q i (),q i ()) q i (),k=2,..., M i=0 to obtain a set of orthogonal polynomials {q 0,q 1,..., q M } which spans the same subspace in C [(, 1)] as 1,, 2,..., Mª. When ρ () 1, we obtain the Legendre polynomials, L 0,L 1,..., L M This set of orthogonal polynomials has many nice properties. One of these nice properties is the fact that the Legendre polynomials satisfy the following recursion 2
3 relation L 0 () = 1,L 1 () = L m+1 () = 2m +1 m +1 L m () m m +1 L m (). This allows us to easily generate the Legendre polynomials. The property of the Legendre polynomials we are most interested in, concerns polynomial interpolation. Let I N map C ([, 1]) onto the polynomials of degree N or less wherei N (f) is the interpolating polynomial for f where the interpolating points are the zeros of the N th Legendre polynomial. Then we have the following. lim I N (f)() f () 2 d =0. N In other words, under rather weak conditions, the sequence of interpolating polynomials will converge mean square to the function f. Our approach to approimating will be to approimate it by f () d I N+1 (f)() d an easily computed integral of a polynomial. (The increase in interpolation points is for notational convenience below.) We can make the computation of the integral still easier by bypassing the computation of the interpolating polynomial I N+1 (f)(). We do that with the theorem below. The quadrature method below [3] is equivalent to simply integrating I N+1 (f)(). Theorem 1 Let 0,..., N be the roots of the (N +1)-th orthogonal polynomial, p N+1, on (a, b) with respect to the weight function ρ and let w 0,..., w N be the solution of the linear system Then ( j ) k w j = Z b a k ρ () d, k =0,...N w j > 0 for j =0,..., N and p ( j ) w j = Z b a p()ρ () d, for all p P 2N+1. The w j are called weights. We will call the roots, j, the collocation points. 3
4 We will find the roots and weights for N =7which we will use below. In this case the approimate roots, calculated on Mathematica, are , , , , , , , and We set up the 8 8 linear system in Matlab and solved for the weights Just to see how well this works, observe that i w i e d = e 1 e and Not too bad. In general we have, 7X ep ( j ) w j lim I N+1 (f)() d = f () d N 3 Two point boundary value problems We will consider the following nonlinear two point boundary value problem L[u] = p 0 () u 00 + p 1 ()u 0 + p 2 () u = g(, u), <<1 B l [u] = α l u() + β l u 0 () = 0 (NLBVP) B r [u] = α r u(1) + β r u 0 (1) = 0 First we will eamine the following simpler boundary value problem. L[u] = p 0 () u 00 + p 1 ()u 0 + p 2 () u = f(), <<1 (BVP) B l [u] = α l u() + β l u 0 () = 0 B r [u] = α r u(1) + β r u 0 (1) = 0 4
5 where p 0 is nonzero and continuously differentiable. It is known, [1], that if the operator L does not have 0 as an eigenvalue, then (BVP) has a unique solution for each continuous f. Thechallengeisfinding the promised solution. There are a wide variety of methods for approimating a solution to (BVP) such as shooting methods and finite difference methods. However, classically, one sought a Green s function and found the solution to (BVP) as an integral. Here we give the integral equation form. u() = K[f] (IEF) = G (, s) f(s)ds. We will now do an eample in which we find the Green s function G. Green s function shall be used in our numerical eamples later on. Consider the following problem This u 00 = f (EX1) u 0 () = 0 u (1) = 0 We may interpret u physically as the temperature in a bar with an insulated left hand end, a right hand end held a fied temperature of 0, and a distributed heat source of intensity f. We will now do various manipulations to put (EX1) into integral equation form (IEF). First we integrate from to to obtain u 0 () u 0 () = Using the boundary condition at yields We now integrate from to 1. u 0 () = u (1) u () = Z Using the boundary condition at 1 yields u () = We now integrate by parts setting w = Z f (t) dt Z s Z s Z s and dv = ds f (t) dt f (t) dtds f (t) dtds f (t) dt 5
6 Then " Z s u () = s = = = f (t) dt 1 f (t) dt f (s) ds + G (, s) f (s) ds Z Z # sf (s) ds f (t) dt f (s) ds + sf (s) ds sf (s) ds where ½ 1 s< G (, s) = s 1 s>. We are now ready to return to (NLBVP). The idea at this point is very simple. If u is a solution to (NLBVP), then setting f () =g (, u) will yield that u also satisfies (BPV ) and consequently or u() = u() = G (, s) f(s)ds. G (, s) g (s, u (s)) ds, (NLIEF) which is a Fredholm integral equation of the first kind. This is a nonlinear integral equation. Since we will be working in situations where the solution to (NLBVP) is unique, a solution to (NLIEF) is the unique solution to (NLBVP). In the net section we put together what we have done so far and give an algorithm for the numerical solution of (NLIEF). 4 The numerical approach Let 0,..., N be the roots of the (N +1)-th Legendre polynomial and w 0,..., w N be the corresponding weights so that for a continuous function f Z b a f()d f ( j ) w j We consider (NLIEF) and apply the approimation u() = G (, s) g (s, u (s)) ds G (, j ) g ( j,u( j )) w j. 6
7 Since this approimation is valid for every in [, 1], itistruefor i, i = 0, 1,...N. Thus u( i ) G ( i, j ) g ( j,u( j )) w j,i=0, 1,...N (AIEF) Let us denote the approimation of u( i ) by u i and define u i by requiring (AIEF) to hold eactly for the approimation. That is u i = G ( i, j ) g ( j,u j ) w j,i=0, 1,...N (AE) This a system of N +1 algebraic equations in the N +1 unknowns u 0,u 1,..., u N. For our first eample, we will study a linear problem. u 00 = g(, u) u 0 () = 0 u(1) = 0. where g(u, ) =u + e e The eact solution to this problem is u() =e + 1 e2. e We will use N =7, for eight points. The system turns out to be an 8 8 linear system which we can solve by our favorite method. The results may be 7
8 seen in Figure 1 and Table 1. 0 Figure 1. -The true solution, * the approimate solution y i u( i ) u i error Table 1. Comparison of the true solution with the quadrature approimation For our second eample, we will merely change g to a nonlinear function of u. g(, u) =u In this case the solution is u () = We will use N =7, for eight points. The system turns out to be an 8 8 nonlinear system. The results may be seen in Figure 2 and Table 2. The 8
9 solution to (AE) in this case was found using the built in Matlab function fmins. However, in a numerical analysis class, this could be used as an opportunity to apply a quasi-newton method. 0 Figure 2. -The true solution, *the approimate solution y i u( i ) u i error Table 2. Comparison of the true solution with the quadrature approimation While not as nice as the results for the linear problem, it is still very good for eight points. To see this, let us compare the results obtained for a finite 9
10 difference approimation. The results may be seen in Table 3 and Figure 3. Figure 3. -The true solution, +the approimate solution via finite difference, * via quadrature y i u( i ) u i error Table 3. Comparison of the true solution with the finite difference approimation 5 References References [1] G. Birkhoff and G. C. Rota, Ordinary Differential Equations Second Edition, Blaisdell Publishing Company,
11 [2] R. L. Burden and J. D. Faires, Numerical Analysis, Sith Edition, Brooks/Cole, 1997 [3] C. Canuto, M.Y. Hussaini, A. Quarteroni and T.A. Zang, Spectral Methods in Fluid Dynamic, Springer-Verlag, 1987 [4] G.N.ElnagarandM.Razzaghi,A Pseudospectral Method for Hammerstein Equations, Journal of Mathematical Analysis and Applications 199, (1996). [5] F. B. Hildebrand, Introduction to Numerical Analysis, Dover, [6] G. Sasone, Orthogonal Functions, Dover,
Recall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationMATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors.
MATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors. Orthogonal sets Let V be a vector space with an inner product. Definition. Nonzero vectors v 1,v
More informationarxiv:gr-qc/ v1 6 Sep 2006
Introduction to spectral methods Philippe Grandclément Laboratoire Univers et ses Théories, Observatoire de Paris, 5 place J. Janssen, 995 Meudon Cede, France This proceeding is intended to be a first
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 informationNumerical Solution of Linear Volterra-Fredholm Integral Equations Using Lagrange Polynomials
Numerical of Linear Volterra-Fredholm Integral Equations Using Polynomials Muna M. Mustafa 1 and Iman N. Ghanim 2* Department of Mathematics, College of Science for Women, University of Baghdad, Baghdad-Iraq
More informationConjugate Gradient Method
Conjugate Gradient Method Hung M Phan UMass Lowell April 13, 2017 Throughout, A R n n is symmetric and positive definite, and b R n 1 Steepest Descent Method We present the steepest descent method for
More informationLinear Algebra. and
Instructions Please answer the six problems on your own paper. These are essay questions: you should write in complete sentences. 1. Are the two matrices 1 2 2 1 3 5 2 7 and 1 1 1 4 4 2 5 5 2 row equivalent?
More informationwhich arises when we compute the orthogonal projection of a vector y in a subspace with an orthogonal basis. Hence assume that P y = A ij = x j, x i
MODULE 6 Topics: Gram-Schmidt orthogonalization process We begin by observing that if the vectors {x j } N are mutually orthogonal in an inner product space V then they are necessarily linearly independent.
More informationMath 5630: Conjugate Gradient Method Hung M. Phan, UMass Lowell March 29, 2019
Math 563: Conjugate Gradient Method Hung M. Phan, UMass Lowell March 29, 219 hroughout, A R n n is symmetric and positive definite, and b R n. 1 Steepest Descent Method We present the steepest descent
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
More information(Mathematical Operations with Arrays) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Mathematical Operations with Arrays) Contents Getting Started Matrices Creating Arrays Linear equations Mathematical Operations with Arrays Using Script
More informationMATH39001 Generating functions. 1 Ordinary power series generating functions
MATH3900 Generating functions The reference for this part of the course is generatingfunctionology by Herbert Wilf. The 2nd edition is downloadable free from http://www.math.upenn. edu/~wilf/downldgf.html,
More informationThere are two things that are particularly nice about the first basis
Orthogonality and the Gram-Schmidt Process In Chapter 4, we spent a great deal of time studying the problem of finding a basis for a vector space We know that a basis for a vector space can potentially
More informationNUMERICAL SIMULATION OF LONG WAVE RUNUP ON A SLOPING BEACH*
NUMERICAL SIMULATION OF LONG WAVE RUNUP ON A SLOPING BEACH* * presented at Long Waves Symposium (in parallel with the XXX IAHR Congress) August 5-7, 003, AUTh, Thessaloniki, Greece. by HAKAN I. TARMAN
More informationPositive Definite Kernels: Opportunities and Challenges
Positive Definite Kernels: Opportunities and Challenges Michael McCourt Department of Mathematical and Statistical Sciences University of Colorado, Denver CUNY Mathematics Seminar CUNY Graduate College
More informationMATH 304 Linear Algebra Lecture 34: Review for Test 2.
MATH 304 Linear Algebra Lecture 34: Review for Test 2. Topics for Test 2 Linear transformations (Leon 4.1 4.3) Matrix transformations Matrix of a linear mapping Similar matrices Orthogonality (Leon 5.1
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationHilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.
Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,
More informationSolutions: Problem Set 3 Math 201B, Winter 2007
Solutions: Problem Set 3 Math 201B, Winter 2007 Problem 1. Prove that an infinite-dimensional Hilbert space is a separable metric space if and only if it has a countable orthonormal basis. Solution. If
More informationStudent name: Student ID: TA s name and/or section: MATH 3B (Butler) Midterm II, 20 February 2009
Student name: Student ID: TA s name and/or section: MATH 3B (Butler) Midterm II, 0 February 009 This test is closed book and closed notes. No calculator is allowed for this test. For full credit show all
More informationMATH 235: Inner Product Spaces, Assignment 7
MATH 235: Inner Product Spaces, Assignment 7 Hand in questions 3,4,5,6,9, by 9:3 am on Wednesday March 26, 28. Contents Orthogonal Basis for Inner Product Space 2 2 Inner-Product Function Space 2 3 Weighted
More informationDiagonalizing Matrices
Diagonalizing Matrices Massoud Malek A A Let A = A k be an n n non-singular matrix and let B = A = [B, B,, B k,, B n ] Then A n A B = A A 0 0 A k [B, B,, B k,, B n ] = 0 0 = I n 0 A n Notice that A i B
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationOn compact operators
On compact operators Alen Aleanderian * Abstract In this basic note, we consider some basic properties of compact operators. We also consider the spectrum of compact operators on Hilbert spaces. A basic
More informationA Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations
Mathematics A Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations Mohamed Meabed KHADER * and Ahmed Saied HENDY Department of Mathematics, Faculty of Science,
More informationPreface. Figures Figures appearing in the text were prepared using MATLAB R. For product information, please contact:
Linear algebra forms the basis for much of modern mathematics theoretical, applied, and computational. The purpose of this book is to provide a broad and solid foundation for the study of advanced mathematics.
More informationSOLUTION KEY TO THE LINEAR ALGEBRA FINAL EXAM 1 2 ( 2) ( 1) c a = 1 0
SOLUTION KEY TO THE LINEAR ALGEBRA FINAL EXAM () We find a least squares solution to ( ) ( ) A x = y or 0 0 a b = c 4 0 0. 0 The normal equation is A T A x = A T y = y or 5 0 0 0 0 0 a b = 5 9. 0 0 4 7
More informationConceptual Questions for Review
Conceptual Questions for Review Chapter 1 1.1 Which vectors are linear combinations of v = (3, 1) and w = (4, 3)? 1.2 Compare the dot product of v = (3, 1) and w = (4, 3) to the product of their lengths.
More information1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,
More informationReview and problem list for Applied Math I
Review and problem list for Applied Math I (This is a first version of a serious review sheet; it may contain errors and it certainly omits a number of topic which were covered in the course. Let me know
More informationFunctional Analysis. James Emery. Edit: 8/7/15
Functional Analysis James Emery Edit: 8/7/15 Contents 0.1 Green s functions in Ordinary Differential Equations...... 2 0.2 Integral Equations........................ 2 0.2.1 Fredholm Equations...................
More informationAlgebra Exam Topics. Updated August 2017
Algebra Exam Topics Updated August 2017 Starting Fall 2017, the Masters Algebra Exam will have 14 questions. Of these students will answer the first 8 questions from Topics 1, 2, and 3. They then have
More informationReduction to the associated homogeneous system via a particular solution
June PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 5) Linear Algebra This study guide describes briefly the course materials to be covered in MA 5. In order to be qualified for the credit, one
More informationResearch Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
Applied Mathematics Volume 013, Article ID 591636, 5 pages http://dx.doi.org/10.1155/013/591636 Research Article A Legendre Wavelet Spectral Collocation Method for Solving Oscillatory Initial Value Problems
More informationA Modification in Successive Approximation Method for Solving Nonlinear Volterra Hammerstein Integral Equations of the Second Kind
Journal of Mathematical Extension Vol. 8, No. 1, (214), 69-86 A Modification in Successive Approximation Method for Solving Nonlinear Volterra Hammerstein Integral Equations of the Second Kind Sh. Javadi
More informationhttp://www.springer.com/3-540-30725-7 Erratum Spectral Methods Fundamentals in Single Domains C. Canuto M.Y. Hussaini A. Quarteroni T.A. Zang Springer-Verlag Berlin Heidelberg 2006 Due to a technical error
More information+ y = 1 : the polynomial
Notes on Basic Ideas of Spherical Harmonics In the representation of wavefields (solutions of the wave equation) one of the natural considerations that arise along the lines of Huygens Principle is the
More informationMATH 304 Linear Algebra Lecture 23: Diagonalization. Review for Test 2.
MATH 304 Linear Algebra Lecture 23: Diagonalization. Review for Test 2. Diagonalization Let L be a linear operator on a finite-dimensional vector space V. Then the following conditions are equivalent:
More informationSuperconvergence Results for the Iterated Discrete Legendre Galerkin Method for Hammerstein Integral Equations
Journal of Computer Science & Computational athematics, Volume 5, Issue, December 05 DOI: 0.0967/jcscm.05.0.003 Superconvergence Results for the Iterated Discrete Legendre Galerkin ethod for Hammerstein
More informationHOSTOS COMMUNITY COLLEGE DEPARTMENT OF MATHEMATICS
HOSTOS COMMUNITY COLLEGE DEPARTMENT OF MATHEMATICS MAT 217 Linear Algebra CREDIT HOURS: 4.0 EQUATED HOURS: 4.0 CLASS HOURS: 4.0 PREREQUISITE: PRE/COREQUISITE: MAT 210 Calculus I MAT 220 Calculus II RECOMMENDED
More informationFINITE-DIMENSIONAL LINEAR ALGEBRA
DISCRETE MATHEMATICS AND ITS APPLICATIONS Series Editor KENNETH H ROSEN FINITE-DIMENSIONAL LINEAR ALGEBRA Mark S Gockenbach Michigan Technological University Houghton, USA CRC Press Taylor & Francis Croup
More informationA Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions
Applied Mathematical Sciences, Vol. 5, 2011, no. 23, 1145-1152 A Comparison between Solving Two Dimensional Integral Equations by the Traditional Collocation Method and Radial Basis Functions Z. Avazzadeh
More informationSyllabus for Applied Mathematics Graduate Student Qualifying Exams, Dartmouth Mathematics Department
Syllabus for Applied Mathematics Graduate Student Qualifying Exams, Dartmouth Mathematics Department Alex Barnett, Scott Pauls, Dan Rockmore August 12, 2011 We aim to touch upon many topics that a professional
More information(Refer Slide Time: 2:04)
Linear Algebra By Professor K. C. Sivakumar Department of Mathematics Indian Institute of Technology, Madras Module 1 Lecture 1 Introduction to the Course Contents Good morning, let me welcome you to this
More informationFamily Feud Review. Linear Algebra. October 22, 2013
Review Linear Algebra October 22, 2013 Question 1 Let A and B be matrices. If AB is a 4 7 matrix, then determine the dimensions of A and B if A has 19 columns. Answer 1 Answer A is a 4 19 matrix, while
More information(Linear equations) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Linear equations) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationA numerical study of a technique for shifting eigenvalues of radial basis function differentiation matrices
A numerical study of a technique for shifting eigenvalues of radial basis function differentiation matrices Scott A. Sarra Marshall University and Alfa R.H. Heryudono University of Massachusetts Dartmouth
More informationInterpolating Accuracy without underlying f (x)
Example: Tabulated Data The following table x 1.0 1.3 1.6 1.9 2.2 f (x) 0.7651977 0.6200860 0.4554022 0.2818186 0.1103623 lists values of a function f at various points. The approximations to f (1.5) obtained
More informationApplied Linear Algebra
Applied Linear Algebra Peter J. Olver School of Mathematics University of Minnesota Minneapolis, MN 55455 olver@math.umn.edu http://www.math.umn.edu/ olver Chehrzad Shakiban Department of Mathematics University
More informationChapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
More informationModule name Calculus and Linear Algebra (Maths 2) To provide additional mathematical tools required for core statistical and actuarial modules
MODULE SPECIFICATION UNDERGRADUATE PROGRAMMES KEY FACTS Module name and Linear Algebra (Maths 2) Module code AS2051 School Cass Business School Department or equivalent UG Programme UK credits 20 ECTS
More informationNumerical 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 informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationDiscrete Projection Methods for Integral Equations
SUB Gttttingen 7 208 427 244 98 A 5141 Discrete Projection Methods for Integral Equations M.A. Golberg & C.S. Chen TM Computational Mechanics Publications Southampton UK and Boston USA Contents Sources
More informationON THE CONVERGENCE OF THE FOURIER APPROXIMATION FOR EIGENVALUES AND EIGENFUNCTIONS OF DISCONTINUOUS PROBLEMS
SIAM J. UMER. AAL. Vol. 40, o. 6, pp. 54 69 c 003 Society for Industrial and Applied Mathematics O THE COVERGECE OF THE FOURIER APPROXIMATIO FOR EIGEVALUES AD EIGEFUCTIOS OF DISCOTIUOUS PROBLEMS M. S.
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationF = m a. t 2. stress = k(x) strain
The Wave Equation Consider a bar made of an elastic material. The bar hangs down vertically from an attachment point = and can vibrate vertically but not horizontally. Since chapter 5 is the chapter on
More informationarxiv: v1 [physics.comp-ph] 22 Jul 2010
Gaussian integration with rescaling of abscissas and weights arxiv:007.38v [physics.comp-ph] 22 Jul 200 A. Odrzywolek M. Smoluchowski Institute of Physics, Jagiellonian University, Cracov, Poland Abstract
More informationEFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL EQUATIONS BASED ON THE GENERALIZED LAGUERRE POLYNOMIALS
Journal of Fractional Calculus and Applications, Vol. 3. July 212, No.13, pp. 1-14. ISSN: 29-5858. http://www.fcaj.webs.com/ EFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL
More informationMATH 240 Spring, Chapter 1: Linear Equations and Matrices
MATH 240 Spring, 2006 Chapter Summaries for Kolman / Hill, Elementary Linear Algebra, 8th Ed. Sections 1.1 1.6, 2.1 2.2, 3.2 3.8, 4.3 4.5, 5.1 5.3, 5.5, 6.1 6.5, 7.1 7.2, 7.4 DEFINITIONS Chapter 1: Linear
More informationLINEAR ALGEBRA 1, 2012-I PARTIAL EXAM 3 SOLUTIONS TO PRACTICE PROBLEMS
LINEAR ALGEBRA, -I PARTIAL EXAM SOLUTIONS TO PRACTICE PROBLEMS Problem (a) For each of the two matrices below, (i) determine whether it is diagonalizable, (ii) determine whether it is orthogonally diagonalizable,
More informationMATH 31 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL
MATH 3 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL MAIN TOPICS FOR THE FINAL EXAM:. Vectors. Dot product. Cross product. Geometric applications. 2. Row reduction. Null space, column space, row space, left
More informationA Gauss Lobatto quadrature method for solving optimal control problems
ANZIAM J. 47 (EMAC2005) pp.c101 C115, 2006 C101 A Gauss Lobatto quadrature method for solving optimal control problems P. Williams (Received 29 August 2005; revised 13 July 2006) Abstract This paper proposes
More informationMAT Linear Algebra Collection of sample exams
MAT 342 - Linear Algebra Collection of sample exams A-x. (0 pts Give the precise definition of the row echelon form. 2. ( 0 pts After performing row reductions on the augmented matrix for a certain system
More informationNumerical Integration Schemes for Unequal Data Spacing
American Journal of Applied Mathematics 017; () 48- http//www.sciencepublishinggroup.com/j/ajam doi 10.1148/j.ajam.01700.1 ISSN 330-0043 (Print); ISSN 330-00X (Online) Numerical Integration Schemes for
More informationPRACTICE FINAL EXAM. why. If they are dependent, exhibit a linear dependence relation among them.
Prof A Suciu MTH U37 LINEAR ALGEBRA Spring 2005 PRACTICE FINAL EXAM Are the following vectors independent or dependent? If they are independent, say why If they are dependent, exhibit a linear dependence
More information33A Linear Algebra and Applications: Practice Final Exam - Solutions
33A Linear Algebra and Applications: Practice Final Eam - Solutions Question Consider a plane V in R 3 with a basis given by v = and v =. Suppose, y are both in V. (a) [3 points] If [ ] B =, find. (b)
More informationJune 2011 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations
June 20 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations The topics covered in this exam can be found in An introduction to differential equations
More informationThe 'linear algebra way' of talking about "angle" and "similarity" between two vectors is called "inner product". We'll define this next.
Orthogonality and QR The 'linear algebra way' of talking about "angle" and "similarity" between two vectors is called "inner product". We'll define this next. So, what is an inner product? An inner product
More informationWe wish to solve a system of N simultaneous linear algebraic equations for the N unknowns x 1, x 2,...,x N, that are expressed in the general form
Linear algebra This chapter discusses the solution of sets of linear algebraic equations and defines basic vector/matrix operations The focus is upon elimination methods such as Gaussian elimination, and
More informationMODULE 8 Topics: Null space, range, column space, row space and rank of a matrix
MODULE 8 Topics: Null space, range, column space, row space and rank of a matrix Definition: Let L : V 1 V 2 be a linear operator. The null space N (L) of L is the subspace of V 1 defined by N (L) = {x
More informationVector Spaces. Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms.
Vector Spaces Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms. For each two vectors a, b ν there exists a summation procedure: a +
More informationUniversity of Colorado Denver Department of Mathematical and Statistical Sciences Applied Linear Algebra Ph.D. Preliminary Exam June 8, 2012
University of Colorado Denver Department of Mathematical and Statistical Sciences Applied Linear Algebra Ph.D. Preliminary Exam June 8, 2012 Name: Exam Rules: This is a closed book exam. Once the exam
More informationMath 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Math 520 Exam 2 Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008 Exam 2 will be held on Tuesday, April 8, 7-8pm in 117 MacMillan What will be covered The exam will cover material from the lectures
More informationIterative Methods for Solving A x = b
Iterative Methods for Solving A x = b A good (free) online source for iterative methods for solving A x = b is given in the description of a set of iterative solvers called templates found at netlib: http
More informationMidterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015
Midterm for Introduction to Numerical Analysis I, AMSC/CMSC 466, on 10/29/2015 The test lasts 1 hour and 15 minutes. No documents are allowed. The use of a calculator, cell phone or other equivalent electronic
More informationJim Lambers MAT 610 Summer Session Lecture 1 Notes
Jim Lambers MAT 60 Summer Session 2009-0 Lecture Notes Introduction This course is about numerical linear algebra, which is the study of the approximate solution of fundamental problems from linear algebra
More information1 CHAPTER 4 THERMAL CONDUCTION
1 CHAPTER 4 THERMAL CONDUCTION 4. The Error Function Before we start this chapter, let s just make sure that we are familiar with the error function erf a. We may need it during this chapter. 1 Here is
More informationA Vector Space Justification of Householder Orthogonalization
A Vector Space Justification of Householder Orthogonalization Ronald Christensen Professor of Statistics Department of Mathematics and Statistics University of New Mexico August 28, 2015 Abstract We demonstrate
More informationIntroduction to Data Mining
Introduction to Data Mining Lecture #21: Dimensionality Reduction Seoul National University 1 In This Lecture Understand the motivation and applications of dimensionality reduction Learn the definition
More informationGuide for Ph.D. Area Examination in Applied Mathematics
Guide for Ph.D. Area Examination in Applied Mathematics (for graduate students in Purdue University s School of Mechanical Engineering) (revised Fall 2016) This is a 3 hour, closed book, written examination.
More informationMath 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination
Math 0, Winter 07 Final Exam Review Chapter. Matrices and Gaussian Elimination { x + x =,. Different forms of a system of linear equations. Example: The x + 4x = 4. [ ] [ ] [ ] vector form (or the column
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationNUMERICAL MATHEMATICS AND COMPUTING
NUMERICAL MATHEMATICS AND COMPUTING Fourth Edition Ward Cheney David Kincaid The University of Texas at Austin 9 Brooks/Cole Publishing Company I(T)P An International Thomson Publishing Company Pacific
More informationLaplace Transforms Chapter 3
Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important
More informationQuadratic forms. Here. Thus symmetric matrices are diagonalizable, and the diagonalization can be performed by means of an orthogonal matrix.
Quadratic forms 1. Symmetric matrices An n n matrix (a ij ) n ij=1 with entries on R is called symmetric if A T, that is, if a ij = a ji for all 1 i, j n. We denote by S n (R) the set of all n n symmetric
More informationPart IB - Easter Term 2003 Numerical Analysis I
Part IB - Easter Term 2003 Numerical Analysis I 1. Course description Here is an approximative content of the course 1. LU factorization Introduction. Gaussian elimination. LU factorization. Pivoting.
More informationNumerical Evaluation of Integrals with weight function x k Using Gauss Legendre Quadrature Rules
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-728, p-issn: 29-76X. Volume, Issue Ver. V (May - Jun. 2), PP 9-6 www.iosrjournals.org Numerical Evaluation of Integrals with weight function k Using Gauss
More informationNumerical Methods with MATLAB
Numerical Methods with MATLAB A Resource for Scientists and Engineers G. J. BÖRSE Lehigh University PWS Publishing Company I(T)P AN!NTERNATIONAL THOMSON PUBLISHING COMPANY Boston Albany Bonn Cincinnati
More informationLecture 1: Basic Concepts
ENGG 5781: Matrix Analysis and Computations Lecture 1: Basic Concepts 2018-19 First Term Instructor: Wing-Kin Ma This note is not a supplementary material for the main slides. I will write notes such as
More informationMath 413/513 Chapter 6 (from Friedberg, Insel, & Spence)
Math 413/513 Chapter 6 (from Friedberg, Insel, & Spence) David Glickenstein December 7, 2015 1 Inner product spaces In this chapter, we will only consider the elds R and C. De nition 1 Let V be a vector
More informationA Review of Linear Algebra
A Review of Linear Algebra Gerald Recktenwald Portland State University Mechanical Engineering Department gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab: Implementations
More informationA NOTE ON THE JORDAN CANONICAL FORM
A NOTE ON THE JORDAN CANONICAL FORM H. Azad Department of Mathematics and Statistics King Fahd University of Petroleum & Minerals Dhahran, Saudi Arabia hassanaz@kfupm.edu.sa Abstract A proof of the Jordan
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 7 Interpolation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More information1 angle.mcd Inner product and angle between two vectors. Note that a function is a special case of a vector. Instructor: Nam Sun Wang. x.
angle.mcd Inner product and angle between two vectors. Note that a function is a special case of a vector. Instructor: Nam Sun Wang Define angle between two vectors & y:. y. y. cos( ) (, y). y. y Projection
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
More information