MAT 167: Advanced Linear Algebra

Size: px
Start display at page:

Download "MAT 167: Advanced Linear Algebra"

Transcription

1 < Proem 1 (15 pts) MAT 167: Advanced Linear Agera Fina Exam Soutions (a) (5 pts) State the definition of a unitary matrix and expain the difference etween an orthogona matrix and an unitary matrix. Soution: A unitary matrix is a square matrix of size whose coumn vectors form an orthonorma asis for. In other words, a matrix is unitary if. In other words,. An orthogona matrix is its counterpart for. In other words, an orthogona matrix is a square matrix of size whose coumn vectors form an orthonorma asis for and satisfies!. () (5 pts) Prove if is unitary, then for any "#, we have % "&%()%"&*. Soution: e simpy use the definitions of the 2norm for vectors, the innerproduct, and a unitary matrix. % "! +, " "/./12 "4 2 "/ " "5" "6 +2""/./7%"& Since oth % "&% and %"&* are nonnegative, we can take nonnegative square roots of oth sizes of the aove equaity to get % ":;4)%"&*. (c) (5 pts) Show that the coumn vectors of any unitary matrix 9 are ineary independent. Soution: Let us consider a inear comination of coumns of : < CD is the G th coumn vector of. Using the notation, HIJ < where FE can e rewritten as: HC Mutipying from eft on the oth sides yieds: 9K < L, the aove Therefore, K9 HMHMC are ineary independent.

2 Œ Œ Score of this page: 2 Proem 2 (15 pts) Consider the foowing matrix: SR;TVUX ZU\[^]_ O&P Ù [a]_ R;TVUXc _dfeghid (a) ( pts) Show that O:P is an orthogona matrix. Soution: This can e easiy done y computing O!PO P () ( pts) Determine the Fourier expansion of " are the coumn vectors of O:P. Soution: Let the Fourier expansion of " e Hqwhere < O P "6 < R;TVUX U\[^]_ < O P O&P. k with respect to the asis whose vectors mu\[a]n RKToUp and < are the Fourier coefficients of ", and Hr s < O&P < \t. Since these coumn vectors form an orthonorma asis, the Fourier coefficients can e computed y mutipying on oth sides from eft: O P "5 Thus, the Fourier expansion in this case is: (c) ( pts) Show that O:P *( Soution: The 2norm of any matrix u R;TVUp RKToUp mu\[a]n R;TVUp Ù [a]_ O P O&P 5U\[a]n regardess of the vaues of. HMH ZU\[^]n R;TVUp is simpy the argest singuar vaue v need to do is to compute the argest eigenvaue of u&xu. In this case, O P O&P k {z wqx%y {z {z 1 z Thus the eigenvaues are: O&P we concude: * O&P (d) ( pts) Compute.. Thus the argest singuar vaue v independent of the vaues of. } of u. Thus a we, thus. Therefore Soution: Reca the 1norm of a matrix u is the maximum of the 1norm of the coumn vectors of u, i.e., *u~ ƒ E ˆ *um> Š G. Thus, in this case, the 1norm of the 1st coumn vector and that of the 2nd coumn vector are the same, i.e., O&P R;TVUp 7Œ U\[^] Œ

3 œ e e Œ e Œ Œ Œ O&P R;TVU U\[^]_ š7 )Œ Œ d First, notice that V R;TVU d ZÙ [a] using the asic trig. identities, such as V Ù [a] d RKToU, V. This means that d is periodic with period. Therefore, we R;TVUpœ ony need to Ù check [a]_œ this function over the interva e fe7d. In this interva, ceary and. Thus, R;TVUp U\[^]_ R;TVU fd Thus it is cear that d. The equaity hods if in this d interva. ow, it is aso cear that in this interva, and the equaity hods if. Using the periodicity of e, we concude that nd nd Ÿd and _d the first d equaity d hods if F, and the second equaity hods if. Score of this page: (e) ( pts) Show that O:P O&P %Ž7 for any vaue of. Soution: Reca the norm of a matrix u is the maximum of the 1norm of the row vectors of u, i.e., *u~%ž 1~ ƒ ˆ a *um \K B. Thus, in this case, the 1norm of the 1st row vector and that of the 2nd row vector are the same, i.e., O&P R;TVUp %Ž)Œ (f) Bonus proem (5 pts) Show that e hod O:P U\[a] O&P Œo7. For what vaues of do the equaities [Hint: R;TVUp Check Ù [a]n the periodicity R;TVU :6d of the 1norm in this case, and then use the trig. identity:.] Soution: Let }

4 L ( + DB. ( J F L () o (1J L O Score of this page: 4 Proem (1 pts) Compute the R factorization of the foowing matrix: u Soution: Use the GramSchmidt procedure here. Let Step 1: ª Step 2: ª Step : ª Finay: e get: 1J L 7; J * DB* e the coumn vectors of u. L J ˆ ( + 9./ ( + 9./ ( () o 1J L J L

5 A Score of this page: 5 Proem 4 (1 pts) For every innerproduct space ± over the scaar fied, show that if ² ³ ±, then ² µ is a suspace of ± even if ² itsef is not a suspace of ±. Soution: Take any " ² µ and any ¹. Take any memer º ². Then, y definition of» space we must have +," º}._ +2 4\º}.. ow, the iinearity of the innerproduct, we have: Thus " derived independenty from whether ² +," / º}. +2" º. +2 ( º./ +2¹¼" º./ ¹(+," º./ ¹ {² µ and ¹½"} {²#µ. This impies that ² µ is a suspace of ±, and this was is a suspace of ± or not. This competes the proof.

6 Score of this page: 6 Proem 5 (1 pts) Consider a inear system u&"#1¾, where u1 that satisfies u C. /, " Then show that this system of equations is consistent if and ony if +k¾ À. [Hint: Consider the orthogona decomposition: ÁMSu!DÂiÃÄ u:.], and ¾. for any Soution: 2@Å : Suppose u:"# ¾ is a consistent system. This means that ¾Ä 6ÁMSu!. ow, any vector # with u 6}C means that # rãæsu. But we know that ÁMSu!4» ÃÆSu in. Thus, ceary, ¾», i.e., +k¾q\ À./. 2Ç : Suppose ¾Z»È for any gãæsu!. Then, either ¾ZC or ¾6 gá u:. If ¾r C, then u&" C aways has a soution ")C. If ¾r Á u:, then this means that u&" ¾ is consistent. Thus, either way, u&" }¾ is consistent.

7 Ò s s z Score of this page: 7 Proem 6 (1 pts) Consider the foowing matrix: ué (a) (5 pts) Compute the Singuar Vaue Decomposition (SVD) of u. Soution: Foows the usua strategy. Thus, we have z wqx;y u u Êz {z Êz {z Êz {z z 2 Êz {z {z K` ;, Êz z~ê. Therefore, the singuar vaues are: V z ow, to determine Ë z, we sove the eigenvaueeigenvector proems. For k Thus we can set s Thus z we have For : Thus we can set s ength, z it ecomes For : Thus, we can set s t Ïs } Ð s t{ s ` t ( s t. ow, to compute Ñ, we use the formua: Ò E t, ut we need to normaize it to have a unit ength. t.. t, and after the normaization, to have a unit t. t, and again after the normaization, it ecomes ÓkÔ u E for v EMÕ Ð Ð Ð. Thus, :

8 v Ò A t Score of this page: Thus the SVD of u is: Ò4 Ð Ð Ð () (5 pts) Compute the rank approximation of u. Soution: The rank approximation of u is of course s Ð s Ð Ð Ð t

9 s z Score of this page: 9 Proem 7 (15 pts) Consider the matrices u Ö1 (a) (5 pts) hat is the orthogona proector Soution: t under ØcÉuM u u: u Therefore the image of s Ù () (1 pts) Same questions for Ö. onto Á u:, and what is the proection of the vector t under is: [Hint: Compute ÖMÚ, i.e., the pseudoinverse of Ö, and construct an orthogona proector using Ö Ú.] Soution: In this case, Ö is not of fu rank and Ö Ö Thus, we cannot do ÀÛ ÜÖÄ ÖÝÖ >Ö. Instead we need to compute the pseudoinverse ÖMÚ and construct the orthogona proector Û#7ÖMÖMÚ. To do so, we need to compute first the SVD of Ö. Êz wqx%y {z {z z So, and v, v z. As for the eigenvector for, k So, Þs t, and after the z normaization, we get: Ïs ow for the eigenvector for, t.

10 Ò s Score of this page: 1 Thus, with the normaization in mind, we get ( t. For Ñ, we use the formua Ò ÓKß u to get: k e need to find Ò4 and Ò( orthogona to Ò with *Ò4V~à*Ò4 o~. So, y ooking the components of Ò, we can easiy find, for exampe, Òm s t. Then, Ò can e Ò4 Ïs t. Thus, the SVD of Ö is: Ö Therefore, Ö Ú ËMá Ú Ñ > > k Therefore, Therefore the image of s ¼ÛmÖMÖ Ú Ù t under is:

11 ¹ â ã ã ¹ F â s ä Ð ä â Score of this page: 11 Proem (15 pts) Suppose we are given three data points in the sense of the east squares. (a) (5 pts) rite a system of equation in the form u:" ¾, where "5Ùs ¹ passes through a these three points. Soution: Since U\[^] e easiy written as: when _d d U\[^] * F d _d, F % ow, we want to find the est function to fit these points in the form of 1¹ #â in t as if this function, respectivey, the system of equation can () (5 pts) Sove the eastu\[^] squares proem using the norma equation, and write the soution in the. form of ¹ {â Soution: This is simpy can e soved y Su u: u ¾ Thus the east square soution is:! Ù [a] > o (c) (5 pts) Sove the east squares proem using the reduced R factorization and confirm that your hand computed resut agrees with that of (). Soution: Let the coumn vectors of u e this case is: ã ; _ and. Then the reduced R factorization of u ã s t in. ä + 9./ ã _ ã ( t Thus, u O æå å

12 ¹ â O å ¾ å O å Score of this page: 12 ow, we transform the inconsistent system u:"i ¾ using the orthogona proector å å a consistent system: u:"6 æå :D¾. Using the R factorization of u, this simpifies to: "5 M¾ å : to Therefore, "6 M¾ å Ð o Therefore, the east squares soution ased on the reduced R factorization is the same as that y the norma equation. (d) Bonus Proem (5 pts) After fitting the east squares soution you computed in () and (c) to the actua data, compute the residua (or error) vector and its ength in 2norm. Soution: The error vector is simpy, u:" o o o o Thus its 2norm is: o o Ð

Matrices and Determinants

Matrices and Determinants Matrices and Determinants Teaching-Learning Points A matri is an ordered rectanguar arra (arrangement) of numbers and encosed b capita bracket [ ]. These numbers are caed eements of the matri. Matri is

More information

221B Lecture Notes Notes on Spherical Bessel Functions

221B Lecture Notes Notes on Spherical Bessel Functions Definitions B Lecture Notes Notes on Spherica Besse Functions We woud ike to sove the free Schrödinger equation [ h d r R(r) = h k R(r). () m r dr r m R(r) is the radia wave function ψ( x) = R(r)Y m (θ,

More information

2M2. Fourier Series Prof Bill Lionheart

2M2. Fourier Series Prof Bill Lionheart M. Fourier Series Prof Bi Lionheart 1. The Fourier series of the periodic function f(x) with period has the form f(x) = a 0 + ( a n cos πnx + b n sin πnx ). Here the rea numbers a n, b n are caed the Fourier

More information

Math 124B January 17, 2012

Math 124B January 17, 2012 Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia

More information

Math 124B January 31, 2012

Math 124B January 31, 2012 Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat

More information

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form Exercises Fourier Anaysis MMG70, Autumn 007 The exercises are taken from: Version: Monday October, 007 DXY Section XY of H F Davis, Fourier Series and orthogona functions EÖ Some exercises from earier

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27

Physics 505 Fall 2007 Homework Assignment #5 Solutions. Textbook problems: Ch. 3: 3.13, 3.17, 3.26, 3.27 Physics 55 Fa 7 Homework Assignment #5 Soutions Textook proems: Ch. 3: 3.3, 3.7, 3.6, 3.7 3.3 Sove for the potentia in Proem 3., using the appropriate Green function otained in the text, and verify that

More information

Assignment 7 Due Tuessday, March 29, 2016

Assignment 7 Due Tuessday, March 29, 2016 Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges

More information

Math 220B - Summer 2003 Homework 1 Solutions

Math 220B - Summer 2003 Homework 1 Solutions Math 0B - Summer 003 Homework Soutions Consider the eigenvaue probem { X = λx 0 < x < X satisfies symmetric BCs x = 0, Suppose f(x)f (x) x=b x=a 0 for a rea-vaued functions f(x) which satisfy the boundary

More information

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS

LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL HARMONICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October 7, 202 Prof. Aan Guth LECTURE NOTES 9 TRACELESS SYMMETRIC TENSOR APPROACH TO LEGENDRE POLYNOMIALS AND SPHERICAL

More information

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Eectromagnetism II October, 202 Prof. Aan Guth LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

More information

Pricing Multiple Products with the Multinomial Logit and Nested Logit Models: Concavity and Implications

Pricing Multiple Products with the Multinomial Logit and Nested Logit Models: Concavity and Implications Pricing Mutipe Products with the Mutinomia Logit and Nested Logit Modes: Concavity and Impications Hongmin Li Woonghee Tim Huh WP Carey Schoo of Business Arizona State University Tempe Arizona 85287 USA

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION

ORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,

More information

Wave Equation Dirichlet Boundary Conditions

Wave Equation Dirichlet Boundary Conditions Wave Equation Dirichet Boundary Conditions u tt x, t = c u xx x, t, < x 1 u, t =, u, t = ux, = fx u t x, = gx Look for simpe soutions in the form ux, t = XxT t Substituting into 13 and dividing

More information

6 Wave Equation on an Interval: Separation of Variables

6 Wave Equation on an Interval: Separation of Variables 6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.

More information

FFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection

FFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection FFTs in Graphics and Vision Spherica Convoution and Axia Symmetry Detection Outine Math Review Symmetry Genera Convoution Spherica Convoution Axia Symmetry Detection Math Review Symmetry: Given a unitary

More information

Lecture Notes 4: Fourier Series and PDE s

Lecture Notes 4: Fourier Series and PDE s Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for

More information

BALANCING REGULAR MATRIX PENCILS

BALANCING REGULAR MATRIX PENCILS BALANCING REGULAR MATRIX PENCILS DAMIEN LEMONNIER AND PAUL VAN DOOREN Abstract. In this paper we present a new diagona baancing technique for reguar matrix pencis λb A, which aims at reducing the sensitivity

More information

Lecture Notes for Math 251: ODE and PDE. Lecture 32: 10.2 Fourier Series

Lecture Notes for Math 251: ODE and PDE. Lecture 32: 10.2 Fourier Series Lecture Notes for Math 251: ODE and PDE. Lecture 32: 1.2 Fourier Series Shawn D. Ryan Spring 212 Last Time: We studied the heat equation and the method of Separation of Variabes. We then used Separation

More information

Week 6 Lectures, Math 6451, Tanveer

Week 6 Lectures, Math 6451, Tanveer Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

MA 201: Partial Differential Equations Lecture - 10

MA 201: Partial Differential Equations Lecture - 10 MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary

More information

Framework for functional tree simulation applied to 'golden delicious' apple trees

Framework for functional tree simulation applied to 'golden delicious' apple trees Purdue University Purdue e-pubs Open Access Theses Theses and Dissertations Spring 2015 Framework for functional tree simulation applied to 'golden delicious' apple trees Marek Fiser Purdue University

More information

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased

PHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization

More information

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t. Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes

More information

More Scattering: the Partial Wave Expansion

More Scattering: the Partial Wave Expansion More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction

More information

Investigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l

Investigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l Investigation on spectrum of the adjacency matrix and Lapacian matrix of graph G SHUHUA YIN Computer Science and Information Technoogy Coege Zhejiang Wani University Ningbo 3500 PEOPLE S REPUBLIC OF CHINA

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

Partial permutation decoding for MacDonald codes

Partial permutation decoding for MacDonald codes Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics

More information

Vectors. Teaching Learning Point. Ç, where OP. l m n

Vectors. Teaching Learning Point. Ç, where OP. l m n Vectors 9 Teaching Learning Point l A quantity that has magnitude as well as direction is called is called a vector. l A directed line segment represents a vector and is denoted y AB Å or a Æ. l Position

More information

1 Heat Equation Dirichlet Boundary Conditions

1 Heat Equation Dirichlet Boundary Conditions Chapter 3 Heat Exampes in Rectanges Heat Equation Dirichet Boundary Conditions u t (x, t) = ku xx (x, t), < x (.) u(, t) =, u(, t) = u(x, ) = f(x). Separate Variabes Look for simpe soutions in the

More information

14 Separation of Variables Method

14 Separation of Variables Method 14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt

More information

Srednicki Chapter 51

Srednicki Chapter 51 Srednici Chapter 51 QFT Probems & Soutions A. George September 7, 13 Srednici 51.1. Derive the fermion-oop correction to the scaar proagator by woring through equation 5., and show that it has an extra

More information

Physics 566: Quantum Optics Quantization of the Electromagnetic Field

Physics 566: Quantum Optics Quantization of the Electromagnetic Field Physics 566: Quantum Optics Quantization of the Eectromagnetic Fied Maxwe's Equations and Gauge invariance In ecture we earned how to quantize a one dimensiona scaar fied corresponding to vibrations on

More information

Higher dimensional PDEs and multidimensional eigenvalue problems

Higher dimensional PDEs and multidimensional eigenvalue problems Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u

More information

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia

More information

LA PRISE DE CALAIS. çoys, çoys, har - dis. çoys, dis. tons, mantz, tons, Gas. c est. à ce. C est à ce. coup, c est à ce

LA PRISE DE CALAIS. çoys, çoys, har - dis. çoys, dis. tons, mantz, tons, Gas. c est. à ce. C est à ce. coup, c est à ce > ƒ? @ Z [ \ _ ' µ `. l 1 2 3 z Æ Ñ 6 = Ð l sl (~131 1606) rn % & +, l r s s, r 7 nr ss r r s s s, r s, r! " # $ s s ( ) r * s, / 0 s, r 4 r r 9;: < 10 r mnz, rz, r ns, 1 s ; j;k ns, q r s { } ~ l r mnz,

More information

Legendre Polynomials - Lecture 8

Legendre Polynomials - Lecture 8 Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha

More information

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

u(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0

u(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0 Bocconi University PhD in Economics - Microeconomics I Prof M Messner Probem Set 4 - Soution Probem : If an individua has an endowment instead of a monetary income his weath depends on price eves In particuar,

More information

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT 06269 3009 and Ronad J. Stern Department of Mathematics, Concordia

More information

Coupling of LWR and phase transition models at boundary

Coupling of LWR and phase transition models at boundary Couping of LW and phase transition modes at boundary Mauro Garaveo Dipartimento di Matematica e Appicazioni, Università di Miano Bicocca, via. Cozzi 53, 20125 Miano Itay. Benedetto Piccoi Department of

More information

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain CORRECTIONS TO CLASSICAL PROCEDURES FOR ESTIMATING THURSTONE S CASE V MODEL FOR RANKING DATA Aberto Maydeu Oivares Instituto de Empresa Marketing Dept. C/Maria de Moina -5 28006 Madrid Spain Aberto.Maydeu@ie.edu

More information

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations)

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations) Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =

More information

Multigrid Lecture Notes

Multigrid Lecture Notes Mutigrid Lecture Notes Wenqiang Feng, Steven Wise Abstract Mutigrid MG methods in numerica anaysis are a group of agorithms for soving differentia equations using a hierarchy of discretizations MG method

More information

9. EXERCISES ON THE FINITE-ELEMENT METHOD

9. EXERCISES ON THE FINITE-ELEMENT METHOD 9. EXERCISES O THE FIITE-ELEMET METHOD Exercise Thickness: t=; Pane strain proem (ν 0): Surface oad Voume oad; 4 p f ( x, ) ( x ) 0 E D 0 0 0 ( ) 4 p F( xy, ) Interna constrain: rigid rod etween D and

More information

Product Cosines of Angles between Subspaces

Product Cosines of Angles between Subspaces Product Cosines of Anges between Subspaces Jianming Miao and Adi Ben-Israe Juy, 993 Dedicated to Professor C.R. Rao on his 75th birthday Let Abstract cos{l,m} : i cos θ i, denote the product of the cosines

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

David Eigen. MA112 Final Paper. May 10, 2002

David Eigen. MA112 Final Paper. May 10, 2002 David Eigen MA112 Fina Paper May 1, 22 The Schrodinger equation describes the position of an eectron as a wave. The wave function Ψ(t, x is interpreted as a probabiity density for the position of the eectron.

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling

Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest

More information

Math-Net.Ru All Russian mathematical portal

Math-Net.Ru All Russian mathematical portal Math-Net.Ru A Russian mathematica porta D. Zaora, On properties of root eements in the probem on sma motions of viscous reaxing fuid, Zh. Mat. Fiz. Ana. Geom., 217, Voume 13, Number 4, 42 413 DOI: https://doi.org/1.1547/mag13.4.42

More information

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary

More information

An Introduction to Optimal Control Applied to Disease Models

An Introduction to Optimal Control Applied to Disease Models An Introduction to Optimal Control Applied to Disease Models Suzanne Lenhart University of Tennessee, Knoxville Departments of Mathematics Lecture1 p.1/37 Example Number of cancer cells at time (exponential

More information

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5].

(f) is called a nearly holomorphic modular form of weight k + 2r as in [5]. PRODUCTS OF NEARLY HOLOMORPHIC EIGENFORMS JEFFREY BEYERL, KEVIN JAMES, CATHERINE TRENTACOSTE, AND HUI XUE Abstract. We prove that the product of two neary hoomorphic Hece eigenforms is again a Hece eigenform

More information

An approximate method for solving the inverse scattering problem with fixed-energy data

An approximate method for solving the inverse scattering problem with fixed-energy data J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999

More information

An Example file... log.txt

An Example file... log.txt # ' ' Start of fie & %$ " 1 - : 5? ;., B - ( * * B - ( * * F I / 0. )- +, * ( ) 8 8 7 /. 6 )- +, 5 5 3 2( 7 7 +, 6 6 9( 3 5( ) 7-0 +, => - +< ( ) )- +, 7 / +, 5 9 (. 6 )- 0 * D>. C )- +, (A :, C 0 )- +,

More information

Problem set 6 The Perron Frobenius theorem.

Problem set 6 The Perron Frobenius theorem. Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator

More information

! " # $! % & '! , ) ( + - (. ) ( ) * + / 0 1 2 3 0 / 4 5 / 6 0 ; 8 7 < = 7 > 8 7 8 9 : Œ Š ž P P h ˆ Š ˆ Œ ˆ Š ˆ Ž Ž Ý Ü Ý Ü Ý Ž Ý ê ç è ± ¹ ¼ ¹ ä ± ¹ w ç ¹ è ¼ è Œ ¹ ± ¹ è ¹ è ä ç w ¹ ã ¼ ¹ ä ¹ ¼ ¹ ±

More information

An implicit Jacobi-like method for computing generalized hyperbolic SVD

An implicit Jacobi-like method for computing generalized hyperbolic SVD Linear Agebra and its Appications 358 (2003) 293 307 wwweseviercom/ocate/aa An impicit Jacobi-ike method for computing generaized hyperboic SVD Adam W Bojanczyk Schoo of Eectrica and Computer Engineering

More information

Volume 13, MAIN ARTICLES

Volume 13, MAIN ARTICLES Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi

More information

Chapter 5. Wave equation. 5.1 Physical derivation

Chapter 5. Wave equation. 5.1 Physical derivation Chapter 5 Wave equation In this chapter, we discuss the wave equation u tt a 2 u = f, (5.1) where a > is a constant. We wi discover that soutions of the wave equation behave in a different way comparing

More information

1D Heat Propagation Problems

1D Heat Propagation Problems Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2

More information

Introduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy

Introduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization

More information

ON THE REPRESENTATION OF OPERATORS IN BASES OF COMPACTLY SUPPORTED WAVELETS

ON THE REPRESENTATION OF OPERATORS IN BASES OF COMPACTLY SUPPORTED WAVELETS SIAM J. NUMER. ANAL. c 1992 Society for Industria Appied Mathematics Vo. 6, No. 6, pp. 1716-1740, December 1992 011 ON THE REPRESENTATION OF OPERATORS IN BASES OF COMPACTLY SUPPORTED WAVELETS G. BEYLKIN

More information

Scott Cohen. November 10, Abstract. The method of Block Cyclic Reduction (BCR) is described in the context of

Scott Cohen. November 10, Abstract. The method of Block Cyclic Reduction (BCR) is described in the context of ycic Reduction Scott ohen November, 99 bstract The method of ock ycic Reduction (R) is described in the context of soving Poisson's equation with Dirichet boundary conditions The numerica instabiityof

More information

AFormula for N-Row Macdonald Polynomials

AFormula for N-Row Macdonald Polynomials Journa of Agebraic Combinatorics, 21, 111 13, 25 c 25 Springer Science + Business Media, Inc. Manufactured in The Netherands. AFormua for N-Row Macdonad Poynomias ELLISON-ANNE WILLIAMS North Caroina State

More information

B œ c " " ã B œ c 8 8. such that substituting these values for the B 3 's will make all the equations true

B œ c   ã B œ c 8 8. such that substituting these values for the B 3 's will make all the equations true System of Linear Equations variables Ð unknowns Ñ B" ß B# ß ÞÞÞ ß B8 Æ Æ Æ + B + B ÞÞÞ + B œ, "" " "# # "8 8 " + B + B ÞÞÞ + B œ, #" " ## # #8 8 # ã + B + B ÞÞÞ + B œ, 3" " 3# # 38 8 3 ã + 7" B" + 7# B#

More information

Primal and dual active-set methods for convex quadratic programming

Primal and dual active-set methods for convex quadratic programming Math. Program., Ser. A 216) 159:469 58 DOI 1.17/s117-15-966-2 FULL LENGTH PAPER Prima and dua active-set methods for convex quadratic programming Anders Forsgren 1 Phiip E. Gi 2 Eizabeth Wong 2 Received:

More information

Wavelet shrinkage estimators of Hilbert transform

Wavelet shrinkage estimators of Hilbert transform Journa of Approximation Theory 163 (2011) 652 662 www.esevier.com/ocate/jat Fu ength artice Waveet shrinkage estimators of Hibert transform Di-Rong Chen, Yao Zhao Department of Mathematics, LMIB, Beijing

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

An explicit Jordan Decomposition of Companion matrices

An explicit Jordan Decomposition of Companion matrices An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057

More information

Work and energy method. Exercise 1 : Beam with a couple. Exercise 1 : Non-linear loaddisplacement. Exercise 2 : Horizontally loaded frame

Work and energy method. Exercise 1 : Beam with a couple. Exercise 1 : Non-linear loaddisplacement. Exercise 2 : Horizontally loaded frame Work and energy method EI EI T x-axis Exercise 1 : Beam with a coupe Determine the rotation at the right support of the construction dispayed on the right, caused by the coupe T using Castigiano s nd theorem.

More information

Do Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix

Do Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix VOL. NO. DO SCHOOLS MATTER FOR HIGH MATH ACHIEVEMENT? 43 Do Schoos Matter for High Math Achievement? Evidence from the American Mathematics Competitions Genn Eison and Ashey Swanson Onine Appendix Appendix

More information

ODE Homework 2. Since M y N x, the equation is not exact. 2. Determine whether the following equation is exact. If it is exact, M y N x 1 x.

ODE Homework 2. Since M y N x, the equation is not exact. 2. Determine whether the following equation is exact. If it is exact, M y N x 1 x. ODE Homework.6. Exact Equations and Integrating Factors 1. Determine whether the foowing equation is exact. If it is exact, find the soution pe x sin qdx p3x e x sin qd 0 [.6 #8] So. Let Mpx, q e x sin,

More information

Jackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Jackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.

More information

Data Search Algorithms based on Quantum Walk

Data Search Algorithms based on Quantum Walk Data Search Agorithms based on Quantum Wak Masataka Fujisaki, Hiromi Miyajima, oritaka Shigei Abstract For searching any item in an unsorted database with items, a cassica computer takes O() steps but

More information

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed

More information

Technische Universität Chemnitz

Technische Universität Chemnitz Technische Universität Chemnitz Sonderforschungsbereich 393 Numerische Simuation auf massiv paraeen Rechnern Zhanav T. Some choices of moments of renabe function and appications Preprint SFB393/03-11 Preprint-Reihe

More information

Wavelet Galerkin Solution for Boundary Value Problems

Wavelet Galerkin Solution for Boundary Value Problems Internationa Journa of Engineering Research and Deveopment e-issn: 2278-67X, p-issn: 2278-8X, www.ijerd.com Voume, Issue 5 (May 24), PP.2-3 Waveet Gaerkin Soution for Boundary Vaue Probems D. Pate, M.K.

More information

6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7

6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7 6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17 Soution 7 Probem 1: Generating Random Variabes Each part of this probem requires impementation in MATLAB. For the

More information

Distributed average consensus: Beyond the realm of linearity

Distributed average consensus: Beyond the realm of linearity Distributed average consensus: Beyond the ream of inearity Usman A. Khan, Soummya Kar, and José M. F. Moura Department of Eectrica and Computer Engineering Carnegie Meon University 5 Forbes Ave, Pittsburgh,

More information

Theory and implementation behind: Universal surface creation - smallest unitcell

Theory and implementation behind: Universal surface creation - smallest unitcell Teory and impementation beind: Universa surface creation - smaest unitce Bjare Brin Buus, Jaob Howat & Tomas Bigaard September 15, 218 1 Construction of surface sabs Te aim for tis part of te project is

More information

Maejo International Journal of Science and Technology

Maejo International Journal of Science and Technology Fu Paper Maejo Internationa Journa of Science and Technoogy ISSN 1905-7873 Avaiabe onine at www.mijst.mju.ac.th A study on Lucas difference sequence spaces (, ) (, ) and Murat Karakas * and Ayse Metin

More information

SCOTT PLUMMER ASHTON ANTONETTI

SCOTT PLUMMER ASHTON ANTONETTI 2742 A E E E, UE, AAAMA 3802 231 EE AEA 38,000-cre dversfed federl cmus 41,000 emloyees wth 72+ dfferent gences UE roosed 80-cre mster-lnned develoment 20 home stes 3,000F of vllgestyle retl 100,000 E

More information

arxiv: v1 [cs.it] 29 Nov 2010

arxiv: v1 [cs.it] 29 Nov 2010 SUBMITTED TO IEEE TRANS ON INFORMATION THEORY, 1 arxiv:10116121v1 [csit] 29 Nov 2010 On Beamformer Design for Mutiuser MIMO Interference Channes Juho Par, Student Member, IEEE, Youngchu Sung, Senior Member,

More information

Supplementary Appendix (not for publication) for: The Value of Network Information

Supplementary Appendix (not for publication) for: The Value of Network Information Suppementary Appendix not for pubication for: The Vaue of Network Information Itay P. Fainmesser and Andrea Gaeotti September 6, 03 This appendix incudes the proof of Proposition from the paper "The Vaue

More information

Fitting affine and orthogonal transformations between two sets of points

Fitting affine and orthogonal transformations between two sets of points Mathematica Communications 9(2004), 27-34 27 Fitting affine and orthogona transformations between two sets of points Hemuth Späth Abstract. Let two point sets P and Q be given in R n. We determine a transation

More information

Lecture 6: Moderately Large Deflection Theory of Beams

Lecture 6: Moderately Large Deflection Theory of Beams Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey

More information

arxiv: v3 [math.ac] 7 Dec 2009

arxiv: v3 [math.ac] 7 Dec 2009 SUPERTROPICAL MATRIX ALGEBRA II: SOLVING TROPICAL EQUATIONS ZUR IZHAKIAN AND LOUIS ROWEN arxiv:0902.2159v3 [math.ac] 7 Dec 2009 Abstract. We continue the study of matrices over a supertropica agebra, proving

More information

ζ-determinants of Laplacians with Neumann and Dirichlet boundary conditions

ζ-determinants of Laplacians with Neumann and Dirichlet boundary conditions INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 38 (25) 8967 8977 doi:1.188/35-447/38/41/9 ζ-determinants of Lapacians with Neumann and Dirichet boundary

More information

Inverse-Variance Weighting PCA-based VRE criterion to select the optimal number of PCs

Inverse-Variance Weighting PCA-based VRE criterion to select the optimal number of PCs Preprints of the 18th IFAC Word Congress Miano (Itay) August 28 - September 2, 211 Inverse-Variance Weighting PCA-based VRE criterion to seect the optima number of PCs Baigh Mnassri E Mostafa E Ade Mustapha

More information

Maximizing Sum Rate and Minimizing MSE on Multiuser Downlink: Optimality, Fast Algorithms and Equivalence via Max-min SIR

Maximizing Sum Rate and Minimizing MSE on Multiuser Downlink: Optimality, Fast Algorithms and Equivalence via Max-min SIR 1 Maximizing Sum Rate and Minimizing MSE on Mutiuser Downink: Optimaity, Fast Agorithms and Equivaence via Max-min SIR Chee Wei Tan 1,2, Mung Chiang 2 and R. Srikant 3 1 Caifornia Institute of Technoogy,

More information

Active Learning and Basis Selection for Kernel-Based Linear Models: A Bayesian Perspective

Active Learning and Basis Selection for Kernel-Based Linear Models: A Bayesian Perspective 1 Active Learning and Basis Seection for Kerne-Based Linear Modes: A Bayesian Perspective John Paisey, Xuejun Liao and Lawrence Carin Department of Eectrica and Computer Engineering Duke University, Durham,

More information

ETIKA V PROFESII PSYCHOLÓGA

ETIKA V PROFESII PSYCHOLÓGA P r a ž s k á v y s o k á š k o l a p s y c h o s o c i á l n í c h s t u d i í ETIKA V PROFESII PSYCHOLÓGA N a t á l i a S l o b o d n í k o v á v e d ú c i p r á c e : P h D r. M a r t i n S t r o u

More information

Solution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes...

Solution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes... : Soution Set Seven Northwestern University, Cassica Mechanics Cassica Mechanics, Third Ed.- Godstein November 8, 25 Contents Godstein 5.8. 2. Components of Torque Aong Principa Axes.......................

More information

11.1 One-dimensional Helmholtz Equation

11.1 One-dimensional Helmholtz Equation Chapter Green s Functions. One-dimensiona Hemhotz Equation Suppose we have a string driven by an externa force, periodic with frequency ω. The differentia equation here f is some prescribed function) 2

More information