MATH661 Homework 2 - Multivariate approximation
|
|
- Marian Golden
- 5 years ago
- Views:
Transcription
1 MATH661 Homework - Multivariate approximation Posted: Sep 6 Due: 11:55PM, Sep 0 1 Problem statement Many physical phenomena are determined by the shape of surfaces separating one domain from another. Examples include: ow around a body; inclusions of one phase in another, e.g., air bubble in water; accumulation of defects in a crystal at boundaries, e.g., formation of metal grains. A rst step in modeling such phenomena is to approximate the bounding surface in an accurate and computationally ecient manner. The parametric form of the surface is f: R! R 3, with components (f 1 (u; v); f (u; v); f 3 (u; v)), and the goal is to seek an approximation (g 1 (u; v); g (u; v); g 3 (u; v)). The approximation will be constructed as a superposition of basis functions B jk (u; v) p g(u; v) = X l=0 Xq m=0 c lm B lm (u; v): (1) Formula (1) generalizes the univariate case, e.g., in Newton's form of the global interpolating polynomial n p n (x) = X i=0 i 1 c i B i (x); c i = [f 0 ; :::; f i ]; B i (x) = Y j=0 (x x j ): () A commonly encountered problem is time evolution of a bounding surface, which can be described through a separation of variables approach as p g(t; u; v) = X l=0 Xq m=0 c lm (t) B lm (u; v): (3) A simple example is a sphere with time varying radius r(t), f(t; ; ') = r(t)(cos sin '; sin sin '; cos '); for which p = q = 0, c 00 (t) = r(t), B 00 (; ') = (cos sin '; sin sin '; cos '). A more interesting example is a superposition of spherical harmonics as would occur in the oscillation of a liquid droplet, or the change in an equipotential surface in quantum mechanical descriptions of an atom p r(t; ; ') = X l=0 Xl m= l c lm (t) Y l m (; '); Y l m (; ') = e im' P l m (cos ) with P l m (cos ) denoting the associated Legendre polynomials dened in terms of the ordinary Legendre polynomials P l (x) by P l m (cos ) = ( 1) m (1 x ) m/ d m dx mp l(x): 1
2 This homework will investigate ways of approximating the spherical harmonics Y l m (; '), in conjunction with the univariate approximation of univariate functions c lm (t). The theoretical exercises build up the foundation of construction of orthogonal polynomial bases. The computational part uses these results to approximate the spherical harmonics, and produce physically relevant preditions. Theoretical exercises 1. K&C, 6.8.1, p.404 Solution. A series of plots for various values of is helpful: Mathematica In[8]:= lambda=table[lam,{lam,0,1,.}]; p=plot[{sin[x],lambda x},{x,0,pi/}]; SetDirectory["/home/student/courses/MATH661/homework"]; Export["HWFig1.png",p] HWFig1.png In[9]:= Figure 1. sin(x) (blue), x (yellow) for = 0; 0.; :::; 1 Introduce F (t; ) = sin t t, and apply Chebyshev alternation theorem to obtain system F (; ) = F ; ) sin = 1 ) sin (; ) = 0 ) cos = 0 cos 1 ) + cos sin = 1 (4) (Full credit for obtaining above system of equations). Solution of (4) can be determined numerically (Lesson 1) In[17]:= csi = csi /. FindRoot[(csi+Pi/) Cos[csi] - Sin[csi] == 1,{csi,Pi/ 3}][[1]] In[18]:= lambda = Cos[csi] In[4]:= p=plot[{sin[x],lambda x},{x,0,pi/},axeslabel->{"x","sin(x), lambda x"},gridlines->{{csi,pi/},none}]; SetDirectory["/home/student/courses/MATH661/homework"]; Export["HWFig.png",p]
3 HWFig.png In[5]:= Figure. Best approximant of sin x by linear function is x.. K&C, 6.8., p.404 (include a plot comparing exact function and approximations from exercises 1 and ) Solution. Using same notation as above, the quadratic norm best approximant is the solution of the problem Solve g 0 () = 0, ming; with Z g() = kf k = 0 / (sin x x) dx ; Z 0 / x(sin x x) dx = 0 ) = " Z 0 / x sin x dx #" Z 0 / x dx # 1 ) = 4/ 3 In[7]:= p=plot[{sin[x],lambda x, 4 x/pi^3},{x,0,pi/},axeslabel->{"x", "sin(x), lambda x"},gridlines->{{csi,pi/},none}]; SetDirectory["/home/student/courses/MATH661/homework"]; Export["HWFig3.png",p] HWFig3.png In[8]:= Figure 3. Comparison of k k, k k 1 best approximants of sin x. 3. K&C, 6.8.5, p.404 3
4 Solution. Consider the expansions of f ; g on the orthonormal set U = fu 1 ; :::; u n g f = a 1 u 1 + ::: + a n u n ; g = b 1 u 1 + :::; +b n u n ; a i ; b i R Since U is an orthonormal set obtain hf ; u j i = ha 1 u 1 + ::: + a n u n ; u j i = a j ; hg; u j i = hb 1 u 1 + :::; +b n u n ; u j i = b j ; * X n hf ; gi = a j u j ; X + n b k u k = Xn Xn a j b k hu j ; u k i j=1 k=1 j=1 k=1 n = X a j b j j=1 4. K&C, 6.8.1, p.405 Solution. Apply the Gram-Schmidt orthogonalization procedure on f1; x; x ; x 3 ; x 4 ; x 5 g with scalar product Z 1 hf ; gi = f(x) g(x) dx 1 Using code from Lesson 6 In[9]:= ScProd[f_,g_,w_,x_,a_,b_]:=Integrate[w f g,{x,a,b}]; GS[u_,w_,x_,a_,b_]:=Module[{q={},n=Length[u],v=u,r}, r=identitymatrix[n]; For [i=1, i<=n, i++, r[[i,i]] = ScProd[v[[i]],v[[i]],w,x,a,b]^(1/); AppendTo[q, v[[i]]/r[[i,i]] ]; For [j=i+1, j<=n, j++, r[[i,j]]=scprod[q[[i]],v[[j]],w,x,a,b]; v[[j]]=v[[j]]-r[[i, j]] q[[i]] ]]; Return[Simplify[q]] ]; Legendre = GS[{1,x,x^,x^3,x^4,x^5},1,x,-1,1] ( 1 p ; 70 x + 15) r 3 x; 1 r ) 5 (3 x 1); 1 r 7 x (5 x 3); 3 (35 x4 30 x + 3) p ; In[30]:= p=plot[legendre,{x,-1,1}]; SetDirectory["/home/student/courses/MATH661/homework"]; Export["HWFig4.png",p] HWFig4.png In[31]:= r 11 x (63 x 4 Figure 4. Orthonormal Legendre polynomials 4
5 The monic form of the Legendre polynomials from the text is In[33]:= TableForm[Expand[Legendre / Table[Coefficient[Legendre[[k+1]],x, k],{k,0,5}]]] 1 x x 1 3 x 3 3 x 5 x 4 6 x x 5 10 x3 9 In[34]:= x 1 5. K&C, 6.9.8, p.40. Use both elementary calculus and application of Chebyshev alternation theorem Solution. Since cosh(x) is an even function, the problem can be stated as min a;b max j"(x)j; "(x) = a + 06x61 bx cosh(x): Extrema of ": [0; 1]! R are either endpoint values "(0),"(1) (since [0; 1] is compact), or local maxima/minima "(t), with t solutions within (0,1) of " 0 (x) = bx sinh(x) = 0: The possible extrema are therefore y 1 (a) = "(0) = a 1; y (a; b) = "(1) = a + b cosh(1); y 3 (a; b) = "(t) = a + bt cosh t; bt sinh(t) = 0: From the Chebyshev alternating theorem (equioscillation theorem) y 1 (a) = y 3 (a; b) = y (a; b); and the best inf-norm approximant is b = cosh(1) 1, and a from solution of system 6. K&C, , p < a = 1 bt + cosh t bt = sinh(t) : b = cosh(1) 1 Solution. Introduce basis sets fa 1 (x); :::; a p (x)g, fb 1 (y); :::; b q (y)g of U ; V respectively. The dimension of a vector (linear) space is equal to the number of basis vectors (functions). A basis set for U V is B = fa 1 (x)b 1 (y); :::; a 1 (x)b q (y); a (x)b 1 (y); :::; a p (x)b q (yg, with pq elements, since elements of B are independent, and any w(x; y) = P i=1 : p Pq j=1 w ij a i (x)b j (y) 5
6 3 Computational application The following tasks will be fully formulated and partially solved in classroom work. 1. Generate data. Choose p; q f; 3; 4; 5; 6g, q 6 p, and coecients c lm (t) = a lm sin(! lm t), a lm random numbers uniformly distributed in [0,1],! lm = (l + m), l = 0; :::; p, m = 0; :::; q. Compute the positions of points on the surface at times t l = lt, t = 1/N, N = 3, l = 0; 1; :::; N, and angles j = j /P ; ' k = k /P, P = 18. Solution. Generate the random mode amplitudes from pylab import *; p=3; q=3; c=random((p+1,q+1)); print c [[ ] [ ] [ ] [ ]] omega=zeros((p+1,q+1)); for l in range(p+1): for m in range(q+1): omega[l,m]=*pi*(l+m) Load 3D plotting routines from mpl_toolkits.mplot3d import Axes3D from matplotlib import cm, colors Load the spherical harmonics library (1 bonus point awarded for coding of spherical harmonics instead of loading a library), and evaluate the spherical harmonics at requested ; ' values. from scipy.special import sph_harm N=3; t=linspace(0.,1.,n+1) P=18; theta=linspace(0.,pi,p+1); phi=linspace(0.,*pi,p+1) phi,theta=meshgrid(phi,theta) B=zeros((p+1,q+1,P+1,P+1)) for l in range(p+1): for m in range(l+1): B[l,m] = sph_harm(m,l,theta,phi).real any(isnan(b)) False Plot some of the spherical harmonics fig=plt.figure(figsize=plt.figaspect(1.)) xs=sin(phi)*cos(theta); ys=sin(phi)*sin(theta); zs=cos(phi); ax=fig.add_subplot(111,projection='3d') ax.plot_surface(xs,ys,zs,rstride=1,cstride=1) show(); 6
7 Bcolors=zeros(B.shape); for l in range(p+1): for m in range(l+1): mxb=b[l,m].max(); mnb=b[l,m].min(); Bcolors[l,m] = (B[l,m]-mnB)/max((mxB-mnB),B[l,m,0,0]) fig=plt.figure(figsize=plt.figaspect(1.)) ax=fig.add_subplot(111,projection='3d') ax.plot_surface(xs,ys,zs,rstride=1,cstride=1, facecolors=cm.rainbow(bcolors[1,1])) show(); fig=plt.figure(figsize=plt.figaspect(1.)) ax=fig.add_subplot(111,projection='3d') ax.plot_surface(xs,ys,zs,rstride=1,cstride=1, facecolors=cm.rainbow(bcolors[,1])) show(); The plots were save to disk and are show in Fig. 1. Figure 5. Some of the spherical harmonics Generate the surface at the requested times r=zeros((n+1,p+1,p+1)); for n in range(n+1): for l in range(p+1): for m in range(l+1): r[n] = r[n] + c[l,m]*b[l,m] Transform the spherical coordinates to Cartesian coordinates for plotting x=zeros(r.shape); y=zeros(r.shape); z=zeros(r.shape) for n in range(n+1): x[n] = r[n]*sin(phi)*cos(theta) y[n] = r[n]*sin(phi)*sin(theta) z[n] = r[n]*cos(phi) fig=figure(figsize=plt.figaspect(1.)) ax=fig.add_subplot(111,projection='3d') ax.plot_surface(x[0],y[0],z[0],rstride=1,cstride=1) show(); 7
8 Figure 6. Surface at t = 0.. Investigate recovery of time-varying coecients by projection onto the spherical harmonic basis. 3. Construct a piecewise linear approximation of the spherical harmonics. 4. Investigate approximation of the time-varying coecients by projection of the data onto the piecewise linear basis. Compare to exact coecients. 8
Approximation theory
Approximation theory Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 1 1.3 6 8.8 2 3.5 7 10.1 Least 3squares 4.2
More informationSection 7.5 Inner Product Spaces
Section 7.5 Inner Product Spaces With the dot product defined in Chapter 6, we were able to study the following properties of vectors in R n. ) Length or norm of a vector u. ( u = p u u ) 2) Distance of
More informationMa 530 Power Series II
Ma 530 Power Series II Please note that there is material on power series at Visual Calculus. Some of this material was used as part of the presentation of the topics that follow. Operations on Power Series
More informationChapter 8: Taylor s theorem and L Hospital s rule
Chapter 8: Taylor s theorem and L Hospital s rule Theorem: [Inverse Mapping Theorem] Suppose that a < b and f : [a, b] R. Given that f (x) > 0 for all x (a, b) then f 1 is differentiable on (f(a), f(b))
More informationTaylor series - Solutions
Taylor series - Solutions. f(x) sin(x) sin(0) + x cos(0) + x x ( sin(0)) +!! ( cos(0)) + + x4 x5 (sin(0)) + 4! 5! 0 + x + 0 x x! + x5 5! x! + 0 + x5 (cos(0)) + x6 6! ( sin(0)) + x 7 7! + x9 9! 5! + 0 +
More informationThis practice exam is intended to help you prepare for the final exam for MTH 142 Calculus II.
MTH 142 Practice Exam Chapters 9-11 Calculus II With Analytic Geometry Fall 2011 - University of Rhode Island This practice exam is intended to help you prepare for the final exam for MTH 142 Calculus
More information1 Ordinary points and singular points
Math 70 honors, Fall, 008 Notes 8 More on series solutions, and an introduction to \orthogonal polynomials" Homework at end Revised, /4. Some changes and additions starting on page 7. Ordinary points and
More informationMATH Max-min Theory Fall 2016
MATH 20550 Max-min Theory Fall 2016 1. Definitions and main theorems Max-min theory starts with a function f of a vector variable x and a subset D of the domain of f. So far when we have worked with functions
More informationMath 153 Calculus III Notes
Math 153 Calculus III Notes 10.1 Parametric Functions A parametric function is a where x and y are described by a function in terms of the parameter t: Example 1 (x, y) = {x(t), y(t)}, or x = f(t); y =
More informationTaylor and Maclaurin Series
Taylor and Maclaurin Series MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Background We have seen that some power series converge. When they do, we can think of them as
More informationAdvanced Calculus Math 127B, Winter 2005 Solutions: Final. nx2 1 + n 2 x, g n(x) = n2 x
. Define f n, g n : [, ] R by f n (x) = Advanced Calculus Math 27B, Winter 25 Solutions: Final nx2 + n 2 x, g n(x) = n2 x 2 + n 2 x. 2 Show that the sequences (f n ), (g n ) converge pointwise on [, ],
More informationLater in this chapter, we are going to use vector functions to describe the motion of planets and other objects through space.
10 VECTOR FUNCTIONS VECTOR FUNCTIONS Later in this chapter, we are going to use vector functions to describe the motion of planets and other objects through space. Here, we prepare the way by developing
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 2 Approximation Theory. Antony Jameson
Advanced Computational Fluid Dynamics AA5A Lecture Approximation Theory Antony Jameson Winter Quarter, 6, Stanford, CA Last revised on January 7, 6 Contents Approximation Theory. Least Squares Approximation
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
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 informationMultivariable Calculus
2 Multivariable Calculus 2.1 Limits and Continuity Problem 2.1.1 (Fa94) Let the function f : R n R n satisfy the following two conditions: (i) f (K ) is compact whenever K is a compact subset of R n. (ii)
More informationi x i y i
Department of Mathematics MTL107: Numerical Methods and Computations Exercise Set 8: Approximation-Linear Least Squares Polynomial approximation, Chebyshev Polynomial approximation. 1. Compute the linear
More informationAnalysis-3 lecture schemes
Analysis-3 lecture schemes (with Homeworks) 1 Csörgő István November, 2015 1 A jegyzet az ELTE Informatikai Kar 2015. évi Jegyzetpályázatának támogatásával készült Contents 1. Lesson 1 4 1.1. The Space
More informationAnnouncements. Topics: Homework: - sections , 6.1 (extreme values) * Read these sections and study solved examples in your textbook!
Announcements Topics: - sections 5.2 5.7, 6.1 (extreme values) * Read these sections and study solved examples in your textbook! Homework: - review lecture notes thoroughly - work on practice problems
More informationMath 221 Notes on Rolle s Theorem, The Mean Value Theorem, l Hôpital s rule, and the Taylor-Maclaurin formula. 1. Two theorems
Math 221 Notes on Rolle s Theorem, The Mean Value Theorem, l Hôpital s rule, and the Taylor-Maclaurin formula 1. Two theorems Rolle s Theorem. If a function y = f(x) is differentiable for a x b and if
More information(1) Recap of Differential Calculus and Integral Calculus (2) Preview of Calculus in three dimensional space (3) Tools for Calculus 3
Math 127 Introduction and Review (1) Recap of Differential Calculus and Integral Calculus (2) Preview of Calculus in three dimensional space (3) Tools for Calculus 3 MATH 127 Introduction to Calculus III
More informationMath Real Analysis II
Math 4 - Real Analysis II Solutions to Homework due May Recall that a function f is called even if f( x) = f(x) and called odd if f( x) = f(x) for all x. We saw that these classes of functions had a particularly
More informationThings you should have learned in Calculus II
Things you should have learned in Calculus II 1 Vectors Given vectors v = v 1, v 2, v 3, u = u 1, u 2, u 3 1.1 Common Operations Operations Notation How is it calculated Other Notation Dot Product v u
More informationAP CALCULUS (AB) Outline Chapter 4 Overview. 2) Recovering a function from its derivatives and a single point;
AP CALCULUS (AB) Outline Chapter 4 Overview NAME Date Objectives of Chapter 4 1) Using the derivative to determine extreme values of a function and the general shape of a function s graph (including where
More information3 Applications of partial differentiation
Advanced Calculus Chapter 3 Applications of partial differentiation 37 3 Applications of partial differentiation 3.1 Stationary points Higher derivatives Let U R 2 and f : U R. The partial derivatives
More informationProgram : M.A./M.Sc. (Mathematics) M.A./M.Sc. (Final) Paper Code:MT-08 Numerical Analysis Section A (Very Short Answers Questions)
Program : M../M.Sc. (Mathematics) M../M.Sc. (Final) Paper Code:MT-08 Numerical nalysis Section (Very Short nswers Questions) 1. Write two examples of transcendental equations. (i) x 3 + sin x = 0 (ii)
More information(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).
1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x
More informationMA 126 CALCULUS II Wednesday, December 10, 2014 FINAL EXAM. Closed book - Calculators and One Index Card are allowed! PART I
CALCULUS II, FINAL EXAM 1 MA 126 CALCULUS II Wednesday, December 10, 2014 Name (Print last name first):................................................ Student Signature:.........................................................
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part IB Thursday 7 June 2007 9 to 12 PAPER 3 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each question in Section
More informationANALYSIS 2: HOMEWORK SOLUTIONS, SPRING 2013
ANALYSIS 2: HOMEWORK SOLUTIONS, SPRING 2013 STEVEN HEILMAN Contents 1. Problem Set 1 1 2. Problem Set 2 9 3. Problem Set 3 16 4. Problem Set 4 30 5. Problem Set 5 44 6. Problem Set 6 62 7. Appendix: Notation
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 informationCompletion Date: Monday February 11, 2008
MATH 4 (R) Winter 8 Intermediate Calculus I Solutions to Problem Set #4 Completion Date: Monday February, 8 Department of Mathematical and Statistical Sciences University of Alberta Question. [Sec..9,
More information1 Solutions to selected problems
1 Solutions to selected problems 1. Let A B R n. Show that int A int B but in general bd A bd B. Solution. Let x int A. Then there is ɛ > 0 such that B ɛ (x) A B. This shows x int B. If A = [0, 1] and
More informationMA CALCULUS II Friday, December 09, 2011 FINAL EXAM. Closed Book - No calculators! PART I Each question is worth 4 points.
CALCULUS II, FINAL EXAM 1 MA 126 - CALCULUS II Friday, December 09, 2011 Name (Print last name first):...................................................... Signature:........................................................................
More informationMTH Calculus with Analytic Geom I TEST 1
MTH 229-105 Calculus with Analytic Geom I TEST 1 Name Please write your solutions in a clear and precise manner. SHOW your work entirely. (1) Find the equation of a straight line perpendicular to the line
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationWaves on 2 and 3 dimensional domains
Chapter 14 Waves on 2 and 3 dimensional domains We now turn to the studying the initial boundary value problem for the wave equation in two and three dimensions. In this chapter we focus on the situation
More informationMath 131 Final Exam Spring 2016
Math 3 Final Exam Spring 06 Name: ID: multiple choice questions worth 5 points each. Exam is only out of 00 (so there is the possibility of getting more than 00%) Exam covers sections. through 5.4 No graphing
More informationIntroduction to Algebraic and Geometric Topology Week 14
Introduction to Algebraic and Geometric Topology Week 14 Domingo Toledo University of Utah Fall 2016 Computations in coordinates I Recall smooth surface S = {f (x, y, z) =0} R 3, I rf 6= 0 on S, I Chart
More informationNumerical Analysis: Interpolation Part 1
Numerical Analysis: Interpolation Part 1 Computer Science, Ben-Gurion University (slides based mostly on Prof. Ben-Shahar s notes) 2018/2019, Fall Semester BGU CS Interpolation (ver. 1.00) AY 2018/2019,
More informationRecall 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 informationContinuity. MATH 161 Calculus I. J. Robert Buchanan. Fall Department of Mathematics
Continuity MATH 161 Calculus I J. Robert Buchanan Department of Mathematics Fall 2017 Intuitive Idea A process or an item can be described as continuous if it exists without interruption. The mathematical
More informationApproximation of Geometric Data
Supervised by: Philipp Grohs, ETH Zürich August 19, 2013 Outline 1 Motivation Outline 1 Motivation 2 Outline 1 Motivation 2 3 Goal: Solving PDE s an optimization problems where we seek a function with
More informationSolution: a) Let us consider the matrix form of the system, focusing on the augmented matrix: 0 k k + 1. k 2 1 = 0. k = 1
Exercise. Given the system of equations 8 < : x + y + z x + y + z k x + k y + z k where k R. a) Study the system in terms of the values of k. b) Find the solutions when k, using the reduced row echelon
More informationNumerical Methods I Orthogonal Polynomials
Numerical Methods I Orthogonal Polynomials Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Nov. 4th and 11th, 2010 A. Donev (Courant Institute)
More informationSYSTEMS OF ODES. mg sin ( (x)) dx 2 =
SYSTEMS OF ODES Consider the pendulum shown below. Assume the rod is of neglible mass, that the pendulum is of mass m, and that the rod is of length `. Assume the pendulum moves in the plane shown, and
More informationCalculus I Announcements
Slide 1 Calculus I Announcements Read sections 3.9-3.10 Do all the homework for section 3.9 and problems 1,3,5,7 from section 3.10. The exam is in Thursday, October 22nd. The exam will cover sections 3.2-3.10,
More informationInterpolation. 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter
Key References: Interpolation 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter 6. 2. Press, W. et. al. Numerical Recipes in C, Cambridge: Cambridge University Press. Chapter 3
More informationSeries Solutions. 8.1 Taylor Polynomials
8 Series Solutions 8.1 Taylor Polynomials Polynomial functions, as we have seen, are well behaved. They are continuous everywhere, and have continuous derivatives of all orders everywhere. It also turns
More informationDr. Allen Back. Oct. 6, 2014
Dr. Allen Back Oct. 6, 2014 Distribution Min=48 Max=100 Q1=74 Median=83 Q3=91.5 mean = 81.48 std. dev. = 12.5 Distribution Distribution: (The letters will not be directly used.) Scores Freq. 100 2 98-99
More information8.7 Taylor s Inequality Math 2300 Section 005 Calculus II. f(x) = ln(1 + x) f(0) = 0
8.7 Taylor s Inequality Math 00 Section 005 Calculus II Name: ANSWER KEY Taylor s Inequality: If f (n+) is continuous and f (n+) < M between the center a and some point x, then f(x) T n (x) M x a n+ (n
More informationFinal Exam Solutions
Final Exam Solutions Laurence Field Math, Section March, Name: Solutions Instructions: This exam has 8 questions for a total of points. The value of each part of each question is stated. The time allowed
More informationInner Product Spaces
Inner Product Spaces Introduction Recall in the lecture on vector spaces that geometric vectors (i.e. vectors in two and three-dimensional Cartesian space have the properties of addition, subtraction,
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationCMU CS 462/662 (INTRO TO COMPUTER GRAPHICS) HOMEWORK 0.0 MATH REVIEW/PREVIEW LINEAR ALGEBRA
CMU CS 462/662 (INTRO TO COMPUTER GRAPHICS) HOMEWORK 0.0 MATH REVIEW/PREVIEW LINEAR ALGEBRA Andrew ID: ljelenak August 25, 2018 This assignment reviews basic mathematical tools you will use throughout
More informationNo books, no notes, no calculators. You must show work, unless the question is a true/false, yes/no, or fill-in-the-blank question.
Math 304 Final Exam (May 8) Spring 206 No books, no notes, no calculators. You must show work, unless the question is a true/false, yes/no, or fill-in-the-blank question. Name: Section: Question Points
More informationMATH 202B - Problem Set 5
MATH 202B - Problem Set 5 Walid Krichene (23265217) March 6, 2013 (5.1) Show that there exists a continuous function F : [0, 1] R which is monotonic on no interval of positive length. proof We know there
More informationChapter 13: General Solutions to Homogeneous Linear Differential Equations
Worked Solutions 1 Chapter 13: General Solutions to Homogeneous Linear Differential Equations 13.2 a. Verifying that {y 1, y 2 } is a fundamental solution set: We have y 1 (x) = cos(2x) y 1 (x) = 2 sin(2x)
More informationStat751 / CSI771 Midterm October 15, 2015 Solutions, Comments. f(x) = 0 otherwise
Stat751 / CSI771 Midterm October 15, 2015 Solutions, Comments 1. 13pts Consider the beta distribution with PDF Γα+β ΓαΓβ xα 1 1 x β 1 0 x < 1 fx = 0 otherwise, for fixed constants 0 < α, β. Now, assume
More informationStaple or bind all pages together. DO NOT dog ear pages as a method to bind.
Math 3337 Homework Instructions: Staple or bind all pages together. DO NOT dog ear pages as a method to bind. Hand-drawn sketches should be neat, clear, of reasonable size, with axis and tick marks appropriately
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More informationInterpolation via weighted l 1 -minimization
Interpolation via weighted l 1 -minimization Holger Rauhut RWTH Aachen University Lehrstuhl C für Mathematik (Analysis) Matheon Workshop Compressive Sensing and Its Applications TU Berlin, December 11,
More informationTAYLOR AND MACLAURIN SERIES
TAYLOR AND MACLAURIN SERIES. Introduction Last time, we were able to represent a certain restricted class of functions as power series. This leads us to the question: can we represent more general functions
More informationx 2 + y 2 = 1. dx 2y and we conclude that, whichever function we chose, dx y2 = 2x + dy dx x2 + d dy dx = 2x = x y sinh(x) = ex e x 2
Implicit differentiation Suppose we know some relation between x and y, e.g. x + y =. Here, y isn t a function of x. But if we restrict attention to y, then y is a function of x; similarly for y. These
More informationSUMMARY PHYSICS 707 Electrostatics. E(x) = 4πρ(x) and E(x) = 0 (1)
SUMMARY PHYSICS 707 Electrostatics The basic differential equations of electrostatics are E(x) = 4πρ(x) and E(x) = 0 (1) where E(x) is the electric field and ρ(x) is the electric charge density. The field
More informationFriday 09/15/2017 Midterm I 50 minutes
Fa 17: MATH 2924 040 Differential and Integral Calculus II Noel Brady Friday 09/15/2017 Midterm I 50 minutes Name: Student ID: Instructions. 1. Attempt all questions. 2. Do not write on back of exam sheets.
More informationExtrema of Functions of Several Variables
Extrema of Functions of Several Variables MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Background (1 of 3) In single-variable calculus there are three important results
More informationScientific Computing
2301678 Scientific Computing Chapter 2 Interpolation and Approximation Paisan Nakmahachalasint Paisan.N@chula.ac.th Chapter 2 Interpolation and Approximation p. 1/66 Contents 1. Polynomial interpolation
More informationMath 113 (Calculus 2) Exam 4
Math 3 (Calculus ) Exam 4 November 0 November, 009 Sections 0, 3 7 Name Student ID Section Instructor In some cases a series may be seen to converge or diverge for more than one reason. For such problems
More informationMATH Final Review
MATH 1592 - Final Review 1 Chapter 7 1.1 Main Topics 1. Integration techniques: Fitting integrands to basic rules on page 485. Integration by parts, Theorem 7.1 on page 488. Guidelines for trigonometric
More informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationOutline. 1 Interpolation. 2 Polynomial Interpolation. 3 Piecewise Polynomial Interpolation
Outline Interpolation 1 Interpolation 2 3 Michael T. Heath Scientific Computing 2 / 56 Interpolation Motivation Choosing Interpolant Existence and Uniqueness Basic interpolation problem: for given data
More information+ 1 for x > 2 (B) (E) (B) 2. (C) 1 (D) 2 (E) Nonexistent
dx = (A) 3 sin(3x ) + C 1. cos ( 3x) 1 (B) sin(3x ) + C 3 1 (C) sin(3x ) + C 3 (D) sin( 3x ) + C (E) 3 sin(3x ) + C 6 3 2x + 6x 2. lim 5 3 x 0 4x + 3x (A) 0 1 (B) 2 (C) 1 (D) 2 (E) Nonexistent is 2 x 3x
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationVANDERBILT UNIVERSITY. MATH 2300 MULTIVARIABLE CALCULUS Practice Test 1 Solutions
VANDERBILT UNIVERSITY MATH 2300 MULTIVARIABLE CALCULUS Practice Test 1 Solutions Directions. This practice test should be used as a study guide, illustrating the concepts that will be emphasized in the
More informationReal Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi
Real Analysis Math 3AH Rudin, Chapter # Dominique Abdi.. If r is rational (r 0) and x is irrational, prove that r + x and rx are irrational. Solution. Assume the contrary, that r+x and rx are rational.
More informationThis exam contains 9 problems. CHECK THAT YOU HAVE A COMPLETE EXAM.
Math 126 Final Examination Autumn 2011 Your Name Your Signature Student ID # Quiz Section Professor s Name TA s Name This exam contains 9 problems. CHECK THAT YOU HAVE A COMPLETE EXAM. This exam is closed
More informationWe consider the problem of finding a polynomial that interpolates a given set of values:
Chapter 5 Interpolation 5. Polynomial Interpolation We consider the problem of finding a polynomial that interpolates a given set of values: x x 0 x... x n y y 0 y... y n where the x i are all distinct.
More informationLECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes).
CE 025 - Lecture 6 LECTURE 6 GAUSS QUADRATURE In general for ewton-cotes (equispaced interpolation points/ data points/ integration points/ nodes). x E x S fx dx hw' o f o + w' f + + w' f + E 84 f 0 f
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
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 informationCURRENT MATERIAL: Vector Calculus.
Math 275, section 002 (Ultman) Spring 2012 FINAL EXAM REVIEW The final exam will be held on Wednesday 9 May from 8:00 10:00am in our regular classroom. You will be allowed both sides of two 8.5 11 sheets
More informationExam 3 MATH Calculus I
Trinity College December 03, 2015 MATH 131-01 Calculus I By signing below, you attest that you have neither given nor received help of any kind on this exam. Signature: Printed Name: Instructions: Show
More informationHomework 27. Homework 28. Homework 29. Homework 30. Prof. Girardi, Math 703, Fall 2012 Homework: Define f : C C and u, v : R 2 R by
Homework 27 Define f : C C and u, v : R 2 R by f(z) := xy where x := Re z, y := Im z u(x, y) = Re f(x + iy) v(x, y) = Im f(x + iy). Show that 1. u and v satisfies the Cauchy Riemann equations at (x, y)
More informationMath 320: Real Analysis MWF 1pm, Campion Hall 302 Homework 2 Solutions Please write neatly, and in complete sentences when possible.
Math 320: Real Analysis MWF 1pm, Campion Hall 302 Homework 2 Solutions Please write neatly, and in complete sentences when possible. Do the following problems from the book: 1.4.2, 1.4.4, 1.4.9, 1.4.11,
More informationSection Taylor and Maclaurin Series
Section.0 Taylor and Maclaurin Series Ruipeng Shen Feb 5 Taylor and Maclaurin Series Main Goal: How to find a power series representation for a smooth function us assume that a smooth function has a power
More informationn=0 ( 1)n /(n + 1) converges, but not
Math 07H Topics for the third exam (and beyond) (Technically, everything covered on the first two exams plus...) Absolute convergence and alternating series A series a n converges absolutely if a n converges.
More informationFormalizing Basic Quaternionic Analysis
Formalizing Basic Quaternionic Analysis Andrea Gabrielli Marco Maggesi University of Florence, Italy Brasilia, ITP 2017 Andrea Gabrielli, Marco Maggesi (UniFI) Formalizing Basic Quaternionic Analysis Brasilia,
More informationMath 0230 Calculus 2 Lectures
Math 00 Calculus Lectures Chapter 8 Series Numeration of sections corresponds to the text James Stewart, Essential Calculus, Early Transcendentals, Second edition. Section 8. Sequences A sequence is a
More informationMath 190: Fall 2014 Homework 4 Solutions Due 5:00pm on Friday 11/7/2014
Math 90: Fall 04 Homework 4 Solutions Due 5:00pm on Friday /7/04 Problem : Recall that S n denotes the n-dimensional unit sphere: S n = {(x 0, x,..., x n ) R n+ : x 0 + x + + x n = }. Let N S n denote
More information10.2-3: Fourier Series.
10.2-3: Fourier Series. 10.2-3: Fourier Series. O. Costin: Fourier Series, 10.2-3 1 Fourier series are very useful in representing periodic functions. Examples of periodic functions. A function is periodic
More informationINFINITE SEQUENCES AND SERIES
11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES In section 11.9, we were able to find power series representations for a certain restricted class of functions. INFINITE SEQUENCES AND SERIES
More informationThis exam will be over material covered in class from Monday 14 February through Tuesday 8 March, corresponding to sections in the text.
Math 275, section 002 (Ultman) Spring 2011 MIDTERM 2 REVIEW The second midterm will be held in class (1:40 2:30pm) on Friday 11 March. You will be allowed one half of one side of an 8.5 11 sheet of paper
More informationHomework for Math , Spring 2013
Homework for Math 3220 2, Spring 2013 A. Treibergs, Instructor April 10, 2013 Our text is by Joseph L. Taylor, Foundations of Analysis, American Mathematical Society, Pure and Applied Undergraduate Texts
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 informationExercises for Multivariable Differential Calculus XM521
This document lists all the exercises for XM521. The Type I (True/False) exercises will be given, and should be answered, online immediately following each lecture. The Type III exercises are to be done
More informationCoordinate systems and vectors in three spatial dimensions
PHYS2796 Introduction to Modern Physics (Spring 2015) Notes on Mathematics Prerequisites Jim Napolitano, Department of Physics, Temple University January 7, 2015 This is a brief summary of material on
More informationFunctional Analysis Exercise Class
Functional Analysis Exercise Class Week: December 4 8 Deadline to hand in the homework: your exercise class on week January 5. Exercises with solutions ) Let H, K be Hilbert spaces, and A : H K be a linear
More informationPower Series Solutions to the Legendre Equation
Department of Mathematics IIT Guwahati The Legendre equation The equation (1 x 2 )y 2xy + α(α + 1)y = 0, (1) where α is any real constant, is called Legendre s equation. When α Z +, the equation has polynomial
More informationMATH 2083 FINAL EXAM REVIEW The final exam will be on Wednesday, May 4 from 10:00am-12:00pm.
MATH 2083 FINAL EXAM REVIEW The final exam will be on Wednesday, May 4 from 10:00am-12:00pm. Bring a calculator and something to write with. Also, you will be allowed to bring in one 8.5 11 sheet of paper
More information