# Register now to access 7 million high quality study materials (What's Course Hero?) Course Hero is the premier provider of high quality online educational resources. With millions of study documents, online tutors, digital flashcards and free courseware, Course Hero is helping students learn more efficiently and effectively. Whether you're interested in exploring new subjects or mastering key topics for your next exam, Course Hero has the tools you need to achieve your goals.

1 Page

### hwk6-5485

Course: MATH 5485, Fall 2009
School: UCF
Rating:

#### Document Preview

5485 Math October 23, 2006 Homework #6 Problems: 3.7 3.8 5b, 6. 3a, 9, 11 (for 9), 12 (for 9), 13. Due: Monday, November 6 Text: B. Bradie, A Friendly Introduction to Numerical Analysis. First Midterm: Wednesday, 1 November Will cover sections 1.24, 2.16, 3.13, 5, 6, 7 (diagonally dominant and tridiagonal matrices), and 8. You will be allowed to use one 8&quot; 11&quot; sheet of notes. No class:...

Register Now

#### Unformatted Document Excerpt

Coursehero >> Florida >> UCF >> MATH 5485

Course Hero has millions of student submitted documents similar to the one
below including study guides, practice problems, reference materials, practice exams, textbook help and tutor support.

Course Hero has millions of student submitted documents similar to the one below including study guides, practice problems, reference materials, practice exams, textbook help and tutor support.
5485 Math October 23, 2006 Homework #6 Problems: 3.7 3.8 5b, 6. 3a, 9, 11 (for 9), 12 (for 9), 13. Due: Monday, November 6 Text: B. Bradie, A Friendly Introduction to Numerical Analysis. First Midterm: Wednesday, 1 November Will cover sections 1.24, 2.16, 3.13, 5, 6, 7 (diagonally dominant and tridiagonal matrices), and 8. You will be allowed to use one 8" 11" sheet of notes. No class: Friday, November 3.
Find millions of documents on Course Hero - Study Guides, Lecture Notes, Reference Materials, Practice Exams and more. Course Hero has millions of course specific materials providing students with the best way to expand their education.

Below is a small sample set of documents:

UCF - MATH - 5485
Math 5485 November 15, 2006Homework #7Problems: 3.10 4.1 4.2 4.3 7 (just do Newton's Method), 11b. 2, 11, 14a, 15a. 2, 8, 10. 1b, 5.Due: Monday, November 27 Text: B. Bradie, A Friendly Introduction to Numerical Analysis.
UCF - MATH - 5485
Math 5485 November 27, 2006Homework #8Problems: 4.4 4.5 1ac, 4a, 5b, 8. 6 (ignore the Wilkinson shift), 12 (compare the convergence rate of the direct QR algorithm with that based on tridiagonalization).Due: Monday, December 4 Text: B. Bradie, A Friend
UCF - MATH - 5485
Math 5485 December 4, 2006Homework #9Problems: 5.1 5.2 5.3 5.4 2, 5, 8. 1b, 4, 9, 12. 4, 8, 11. 2, 3, 10 (only uniform and Chebyshev).Due: Wednesday, December 13 Text: B. Bradie, A Friendly Introduction to Numerical Analysis. Second Midterm: Friday, De
UCF - MATH - 5485
Chapter 17 Dynamics of Planar MediaIn this chapter, we continue our ascent of the dimensional ladder for linear systems. In Chapter 6, we embarked on our journey with equilibrium configurations of discrete systems - massspring chains, circuits, and struc
UCF - MATH - 5485
AIMS Lecture Notes 2006Peter J. Olver7. Iterative Methods for Linear SystemsLinear iteration coincides with multiplication by successive powers of a matrix; convergence of the iterates depends on the magnitude of its eigenvalues. We discuss in some det
UCF - MATH - 5485
AIMS Lecture Notes 2006Peter J. Olver5. Inner Products and NormsThe norm of a vector is a measure of its size. Besides the familiar Euclidean norm based on the dot product, there are a number of other important norms that are used in numerical analysis
UCF - MATH - 5485
Chapter 15 The Planar Laplace EquationThe fundamental partial differential equations that govern the equilibrium mechanics of multi-dimensional media are the Laplace equation and its inhomogeneous counterpart, the Poisson equation. The Laplace equation i
UCF - MATH - 5485
Very Basic MATLABPeter J. Olver January, 2009 Matrices: Type your matrix as follows: Use space or , to separate entries, and ; or return after each row. &gt; A = [4 5 6 -9;5 0 -3 6;7 8 5 0; -1 4 5 1] or &gt; A = [4,5,6,-9;5,0,-3,6;7,8,5,0;-1,4,5,1] or &gt; A = [
UCF - MATH - 5485
AIMS Lecture Notes 2006Peter J. Olver9. Numerical Solution of Algebraic SystemsIn this part, we discuss basic iterative methods for solving systems of algebraic equations. By far the most common is a vector-valued version of Newton's Method, which will
UCF - MATH - 5485
Chapter 19 Nonlinear SystemsNonlinearity is ubiquitous in physical phenomena. Fluid and plasma mechanics, gas dynamics, elasticity, relativity, chemical reactions, combustion, ecology, biomechanics, and many, many other phenomena are all governed by inhe
UCF - MATH - 5485
Chapter 22 Nonlinear Partial Differential EquationsThe ultimate topic to be touched on in this book is the vast and active field of nonlinear partial differential equations. Leaving aside quantum mechanics, which remains to date an inherently linear theo
UCF - MATH - 5485
AIMS Lecture Notes 2006Peter J. Olver11. Numerical Solution of the Heat and Wave EquationsIn this part, we study numerical solution methodss for the two most important equations of one-dimensional continuum dynamics. The heat equation models the diffus
UCF - MATH - 5485
AIMS Lecture Notes 2006Peter J. Olver10. Numerical Solution of Ordinary Differential EquationsThis part is concerned with the numerical solution of initial value problems for systems of ordinary differential equations. We will introduce the most basic
UCF - MATH - 5485
AIMS Lecture Notes 2006Peter J. Olver8. Numerical Computation of EigenvaluesIn this part, we discuss some practical methods for computing eigenvalues and eigenvectors of matrices. Needless to say, we completely avoid trying to solve (or even write down
UCF - MATH - 5485
AIMS Lecture Notes 2006Peter J. Olver13. Approximation and InterpolationWe will now apply our minimization results to the interpolation and least squares fitting of data and functions.13.1. Least Squares.Linear systems with more equations than unknow
UCF - MATH - 5485
AIMS Lecture Notes 2006Peter J. Olver2. Numerical Solution of Scalar EquationsMost numerical solution methods are based on some form of iteration. The basic idea is that repeated application of the algorithm will produce closer and closer approximation
UCF - MATH - 5485
Chapter 20 Nonlinear Ordinary Differential EquationsThis chapter is concerned with initial value problems for systems of ordinary differential equations. We have already dealt with the linear case in Chapter 9, and so here our emphasis will be on nonline
UCF - MATH - 5485
Chapter 18 Partial Differential Equations in ThreeDimensional SpaceAt last we have ascended the dimensional ladder to its ultimate rung (at least for those of us living in a three-dimensional universe): partial differential equations in physical space. A
UCF - MATH - 5485
Orthogonal Bases and the QR Algorithmby Peter J. Olver University of Minnesota1. Orthogonal Bases.Throughout, we work in the Euclidean vector space V = R n , the space of column vectors with n real entries. As inner product, we will only use the dot pr
UCF - MATH - 5485
AIMS Lecture Notes 2006Peter J. Olver14. Finite ElementsIn this part, we introduce the powerful finite element method for finding numerical approximations to the solutions to boundary value problems involving both ordinary and partial differential equa
UCF - MATH - 5485
AIMS Lecture Notes 2006Peter J. Olver12. MinimizationIn this part, we will introduce and solve the most basic mathematical optimization problem: minimize a quadratic function depending on several variables. This will require a short introduction to pos
UCF - MATH - 5587
Remark : On a connected domain R 2 , all harmonic conjugates to a given function u(x, y) only differ by a constant: v(x, y) = v(x, y) + c; see Exercise . Although most harmonic functions have harmonic conjugates, unfortunately this is not always the case.
UCF - MATH - 5587
Chapter 7 Complex Analysis and Conformal MappingThe term &quot;complex analysis&quot; refers to the calculus of complex-valued functions f (z) depending on a single complex variable z. To the novice, it may seem that this subject should merely be a simple reworkin
UCF - MATH - 5587
1 Re z Figure 7.1.1 Im z 1 Real and Imaginary Parts of f (z) = z .Therefore, if f (z) is any complex function, we can write it as a complex combination f (z) = f (x + i y) = u(x, y) + i v(x, y), of two inter-related real harmonic functions: u(x, y) = Re
UCF - MATH - 5587
Figure 7.4.Real and Imaginary Parts ofz.also have complex logarithms! On the other hand, if z = x &lt; 0 is real and negative, then log z = log | x | + (2 k + 1) i is complex no matter which value of ph z is chosen. (This explains why one avoids defining
UCF - MATH - 5587
The proof of the converse - that any function whose real and imaginary components satisfy the CauchyRiemann equations is differentiable - will be omitted, but can be found in any basic text on complex analysis, e.g., [3, 65, 118]. Remark : It is worth poi
UCF - MATH - 5587
is analytic everywhere except for singularities at the points z = 3 and z = -1, where its denominator vanishes. Since f (z) = h1 (z) , z-3 where h1 (z) = ez (z + 1)21 is analytic at z = 3 and h1 (3) = 16 e3 = 0, we conclude that z = 3 is a simple (order
UCF - MATH - 5587
if and only if it has vanishing divergence: v = u v + = 0. x y (7.36)Incompressibility means that the fluid volume does not change as it flows. Most liquids, including water, are, for all practical purposes, incompressible. On the other hand, the flow is
UCF - MATH - 5587
Using formula (7.19) for the complex derivative, d = -i = u - i v, dz x y so = u, x = v. yThus, = v, and hence the real part (x, y) of the complex function (z) defines a velocity potential for the fluid flow. For this reason, the anti-derivative (z) is k
UCF - MATH - 5587
gDFigure 7.14.Mapping to the Unit Disk.Remark : In this section, we have focused on the fluid mechanical roles of a harmonic function and its conjugate. An analogous interpretation applies when (x, y) represents an electromagnetic potential function;
UCF - MATH - 5587
Figure 7.16.The Effect of = z 2 on Various Domains.obtained by cutting the complex plane along the negative real axis. On the other hand, vertical lines Re z = a are mapped to circles | | = ea . Thus, a vertical strip a &lt; Re z &lt; b is mapped to an annulu
UCF - MATH - 5587
zph zFigure 7.18.Complex Curve and Tangent.notation x(t) = ( x(t), y(t) ) to complex notation z(t) = x(t)+ i y(t). All the usual vectorial curve terminology - closed, simple (non-self intersecting), piecewise smooth, etc. - is employed without modific
UCF - MATH - 5587
Center: .1 Radius: .5Center: .2 + i Radius: 1Center: 1 + i Radius: 1Center: -2 + 3 i Radius: 3 2 4.2426Center: .2 + i Radius: 1.2806 Figure 7.21.Center: .1 + .3 i Radius: .9487Center: .1 + .1 i Radius: 1.1045Center: -.2 + .1 i Radius: 1.2042Airfoi
UCF - MATH - 5587
Example 7.35. The goal of this example is to construct an conformal map that takes a half disk D+ = | z | &lt; 1, Im z &gt; 0 (7.73) to the full unit disk D = cfw_ | | &lt; 1 . The answer is not = z 2 because the image of D+ omits the positive real axis, resulting
UCF - MATH - 5587
7.5. Applications of Conformal Mapping.Let us now apply what we have learned about analytic/conformal maps. We begin with boundary value problems for the Laplace equation, and then present some applications in fluid mechanics. We conclude by discussing h
UCF - MATH - 5587
Figure 7.25.A NonCoaxial Cable.Example 7.39. A non-coaxial cable. The goal of this example is to determine the electrostatic potential inside a non-coaxial cylindrical cable, as illustrated in Figure 7.25, with prescribed constant potential values on th
UCF - MATH - 5587
0 Figure 7.29.15 Fluid Flow Past a Tilted Plate.30Note that = ( 1, 0 ), and hence this flow satisfies the Neumann boundary conditions (7.95) on the horizontal segment D = . The corresponding complex potential is (z) = z, with complex velocity f (z) = (
UCF - MATH - 5587
on the unit disk D for an impulse concentrated at the origin; see Section 6.3 for details. How do we obtain the corresponding solution when the unit impulse is concentrated at another point = + i D instead of the origin? According to Example 7.25, the lin
UCF - MATH - 5587
as long as n = -1. Therefore, we can use the Fundamental Theorem of Calculus (which works equally well for real integrals of complex-valued functions), to evaluate n+1 1 -1 = n = 2 k + 1 odd, 0, 2 t + i (t - 1) 2 z n dz = = , n = 2 k even. n+1 P t = -1 n+
UCF - MATH - 5587
Figure 7.32.Orientation of Domain Boundary.Theorem 7.48. If f (z) is analytic on a bounded domain C, then f (z) dz = 0.(7.118)Proof : If we apply Green's Theorem to the two real line integrals in (7.109), we find u dx - v dy = - u v - x y = 0,v dx +
UCF - MATH - 5587
Proof : Note that the integrand f (z) = 1/(z - a) is analytic everywhere except at z = a, where it has a simple pole. If a is outside C, then Cauchy's Theorem 7.48 applies, and the integral is zero. On the other hand, if a is inside C, then Proposition 7.
UCF - MATH - 5587
0 Figure 7.36.15 Kutta Flow Past a Tilted Airfoil.30which remains asymptotically 1 at large distances. By Cauchy's Theorem 7.48 coupled with formula (7.123), if C is a curve going once around the disk in a counter-clockwise direction, then i 1 dz = - 2
UCF - MATH - 5587
is analytic in the disk | z | 2 since its only singularity, at z = 3, lies outside the contour C. Therefore, by Cauchy's formula (7.135), we immediately obtain the integral ez dz = z2 - 2 z - 3 f (z) i dz = 2 i f (-1) = - . z+1 2eCCNote: Path independe
UCF - MATH - 5587
Chapter 12 Dynamics of Planar MediaIn previous chapters we studied the equilibrium configurations of planar media - plates and membranes - governed by the two-dimensional Laplace and Poisson equations. In this chapter, we analyze their dynamics, modeled
UCF - MATH - 5587
In this manner, we arrive at the basic conservation law relating the heat energy density and the heat flux vector w. As in our one-dimensional model, cf. (4.3), the heat energy density (t, x, y) is proportional to the temperature, so (t, x, y) = (x, y) u(
UCF - MATH - 5587
for the diffusion equation. See [35; p. 369] for a precise statement and proof of the general theorem. Qualitative Properties Before tackling examples in which we are able to construct explicit formulae for the eigenfunctions and eigenvalues, let us see w
UCF - MATH - 5587
Theorem 12.1. Suppose u(t, x, y) is a solution to the forced heat equation ut = u + F (t, x, y), for (x, y) , 0 &lt; t &lt; c,where is a bounded domain, and &gt; 0. Suppose F (t, x, y) 0 for all (x, y) and 0 t c. Then the global maximum of u on the set cfw_ (t, x
UCF - MATH - 5587
so there are no non-separable eigenfunctions . As a consequence, the general solution to the initial-boundary value problem can be expressed as a linear combination u(t, x, y) =m,n = 1cm,n um,n (t, x, y) =m,n = 1cm,n e- m,n t vm,n (x, y)(12.41)of
UCF - MATH - 5587
Let us start with the equation for q(). The second boundary condition in (12.50) requires that q() be 2 periodic. Therefore, the required solutions are the elementary trigonometric functions q() = cos m or sin m , where = m2 , (12.53)with m = 0, 1, 2, .
UCF - MATH - 5587
15 10 5 -4 -2 -5 -10 -15 2 4Figure 12.3.The Gamma Function.Thus, at integer values of x, the gamma function agrees with the elementary factorial. A few other values can be computed exactly. One important case is when x = 1 . Using 2 the substitution t
UCF - MATH - 5587
Remark : The definition of a singular point assumes that the other coefficients do not both vanish there, i.e., either q(x0 ) = 0 or r(x0 ) = 0. If all three functions happen to vanish at x0 , we can cancel any common factor (x - x0 )k , and hence, withou
UCF - MATH - 5587
we find that the only non-zero coefficients un are when n = 3 k +1. The recurrence relation u3 k+1 = u3 k-2 (3 k + 1)(3 k) yields u3 k+1 = 1 . (3 k + 1)(3 k)(3 k - 2)(3 k - 3) 7 6 4 3The resulting solution isx3 k+1 . (3 k + 1)(3 k)(3 k - 2)(3 k - 3) 7 6
UCF - MATH - 5587
two different Frobenius expansions. Usually, this expectation is valid, but there is an important exception, which occurs when the indices differ by an integer. The general result is summarized in the following list: (i ) If r2 - r1 is not an integer, the
UCF - MATH - 5587
We have thus found the series solution (-1)k xm+2k . 22k k(k - 1) 3 2 (r + k)(r + k - 1) (r + 2)(r + 1) k=0 k=0 (12.93) So far, we not paid attention to the precise values of the indices r = m. In order to continue the recurrence, we need to ensure that t
UCF - MATH - 5587
where h0 = 0, while = limkhk = 1 +1 1 1 + + + , 2 3 k (12.102)hk - log k .5772156649 . . .is known as Euler's constant. All Bessel functions of the second kind have a singularity at the origin x = 0; indeed, by inspection of (12.101), we find that th
UCF - MATH - 5587
of the Bessel boundary value problem (12.5455) are the squares of the roots of the Bessel function of order m. The corresponding eigenfunctions are wm,n (r) = Jm (m,n r) , n = 1, 2, 3, . . . , m = 0, 1, 2, . . . , (12.112)defined for 0 r 1. Combining (12
UCF - MATH - 5587
t=0t = .04t = .08t = .12 Figure 12.6.t = .16 Heat Diffusion in a Disk.t = .212.5. The Fundamental Solution of the Heat Equation.As we learned in Section 4.1, the fundamental solution to the heat equation measures the temperature distribution result
UCF - MATH - 5587
for the planar heat equation is given by the linear superposition formula u(t, x, y) = 1 4 t f (, ) e- [ (x-)2+(y-)2 ]/(4 t)d d.(12.125)We can interpret the solution formula (12.125) as a two-dimensional convolution u(t, x, y) = F (t, x, y) f (x, y)
UCF - MATH - 5587
Vibration of a Rectangular Drum Let us first consider the vibrations of a membrane in the shape of a rectangle R= 0 &lt; x &lt; a, 0 &lt; y &lt; b ,with side lengths a and b, whose edges are fixed to the (x, y)plane. Thus, we seek to solve the wave equation utt = c2
UCF - MATH - 5587
A table of their values (for the case c = 1) can be found in the preceding section. The Bessel roots do not follow any easily discernible pattern, and are not rational multiples of each other. This result, known as Bourget's hypothesis, [142; p. 484], was
UCF - MATH - 5587
following table, we display a list of all relative vibrational frequencies (12.158) that are &lt; 6. Once the lowest frequency 0,1 has been determined - either theoretically, numerically or experimentally - all the higher overtones m,n = m,n 0,1 are simply o