Numerical Analysis Using M#math
|
|
- Frank McDowell
- 6 years ago
- Views:
Transcription
1 "msharpmath" revised on Numerical Analysis Using M#math Chapter 9 Eigenvalue Problems 9-1 Existence of Solution and Eigenvalues 9-2 Sturm-Liouville Problems 9-3 Eigenvalues on Stability and Convergence Rate 9-4 Exercises ============================== 9-5 Higher-order ODE and Eigenvalue Problems (optional) Eigenvalue problems occurring in electromagnetics, vibration, structure analysis and so on are related to the boundary value problems of ordinary differential equations. The eigenvalue-finding methods discussed in Chapter 4 provide the necessary platform to solve the eigenvalue problems of ordinary differential equations. Nevertheless, eigenvalues for relevant to an matrix consist of at most eigenvalues. This is in contrast to the fact that the eigenvalue problems are associated with an infinite number of eigenvalues. Therefore, there occurs an intrinsic limitation in obtaining eigenvalues for eigenvalue problems such as the well-known Sturm-Liouville problems. However, fortunately the fundamental eigenvalue that has the smallest absolute value is the most important in applications. And engineering importance diminishes as the absolute value of eigenvalues increases. In this regard, the methods treated in this chapter are of practical importance even though only a finite number of eigenvalues can be determined. First of all, discussed is the situation where an infinite number of solutions are existing. Next, the eigenvalue problems are properly defined. Then, numerical approaches to obtain eigenvalues are discussed. In doing this, the generalized nonlinear differential operator adopted earlier 1
2 is exploited again and the eigenvalue problems are expressed in the following compact form In the last section, the role of eigenvalues is explored in the initial value problems and in the iteration methods. //====================================================================== // syntax in M#math // Umbrella bvp is designed to solve BVPs the general forms of which are cz*{f''''} +cw*{f'''} +cfw*{ff'''} +cuw*{f'f'''} + cvw*{f''f'''} + cww*{f'''f'''} +cv*{f''} +cfv*{ff''} +cuv*{f'f''} + cvv*{f''f''} +cu*{f'} +cfu*{ff'} +cuu*{f'f'} +cf*{f} +cff*{ff} = cs bw*{f'''} + bv*{f''} + bu*{f'} + bf*{f} = bs Syntax bvp.x[n=101,g=1](a,b) ( DE, [ BCs ] ).Spoke.Spoke... x = (a,b).span(n,g=1); bvp [x] ( DE, [BCs] ).Spoke.Spoke... Spokes.peep.relax(relax=0.5).tol/abstol(abstol=1.e-5).reltol(reltol=1.e-5).maxiter(maxiter=500).shift(scalar).plot.plot(f,f',f'',...).return (f,f',f'',...) // peep the iteration procedure // under-relaxation factor // absolute tolerance // relative tolerance // maximum iteration count // eigenvalue shift // plot default data // plot user-defined data // return user-defined data 2
3 .togo (Y=y,Y1=y',...) // takeout each field Section 9-1 Existence of Solution and Eigenvalues In this section, the existence of an infinite number of solutions are discussed focusing on second-order ODEs, and the eigenvalue problems are defined. Existence of an infinite number of solutions. As was discussed in Chapter 8, solutions to boundary value problems are either 1 unique, 2 nonexisting, or 3 of an infinite number. Nevertheless, despite the nonexistence of mathematical solution, a meaningless numerical solution is always obtained due to the inherent limitation of computing with finite number of digits, as was confirmed in Example 8-5. As a beginning step, let us investigate the case where an infinite number of solutions are existing. For a particular second-order ODE, there are two types of boundary conditions with which an infinite number of solutions are obtained. infinite number of solutions with multiple integration constants infinite number of solutions with a single integration constant For the cases of equations (1) and (2), only a trivial solution is obtained numerically. This is because both the ODE and boundary conditions are homogeneous. Homogeneous conditions and eigenvalue problems. As was discussed above, the ODE with non-homogeneous boundary conditions cannot have a trivial solution. In this case, only one non-trivial solution is (1) (2) (3) 3
4 numerically obtained even though there exist an infinite number of solutions [the case of equation (3)]. The eigenvalue problems for the ODE over an interval treat only homogeneous boundary conditions and only nontrivial solutions. In the above, are not zero simultaneously. The homogeneous boundary condition guarantees that a constant multiple of an eigenfunction is also an eigenfunction. Of course, the ODE relevant to the eigenvalue problem is type of homogeneous differential equation. Eigenvalue problems, definition and property. Especially for the secondorder ODE, the eigenvalue problems are generalized as the Sturm-Liouville problem the detail of which will be discussed in the next section. Here, the characteristics of eigenvalue problems are addressed briefly. Let us find the mathematical solution to the following second-order ODE with homogeneous differential equation and boundary conditions (4) (5) First for the case of, the general solution is given as so that from the boundary condition. And from the boundary condition, we must have. However, since, it follows that. In other words, a trivial solution is the only possible solution. Next for the case of, from the general solution and from the boundary conditions, the only possible solution is found to be the trivial solution. Last for the case of, the general solution (6) (7) 4
5 is utilized to find from the boundary condition. This means that. And the remaining boundary condition enforces that (8) which can be satisfied by. At this point, the corresponding value is called the eigenvalue, and the corresponding solution is called the eigenfunction. Meanwhile, it is conventional that the integration constant is set to be unity, unless otherwise specified. From the above simple example, the eigenvalue problems are understood to exhibit the following characteristics. homogeneous ODE homogeneous boundary condition infinite number of eigenvalues eigenfunction corresponding to an eigenvalue, a constant multiple of, is also an eigenfunction the fundamental eigenvalue of smallest absolute value Next, let us investigate the trend of eigenfunctions in response to the increase of their absolute values. The graphical solution to the above-studied eigenvalue problem (9) is shown in Figure
6 Figure 9-1 Eigenfunctions corresponding to The eigenfunction corresponding to the smallest absolute eigenvalue is especially termed to be the fundamental eigenfunction. As shown in Figure 9-1, the fundamental eigenfunction does not have any root that satisfies inside the domain. Similarly, the second, third and fourth eigenfunctions have one, two and three roots of, respectively, inside the domain. This property is valid for other types of boundary conditions. For example, the profiles of eigenfunctions corresponding to the boundary conditions are displayed in Figure 9-2. It can be easily seen that the number of roots of is the same as that shown in Figure 9-1. Figure 9-2 Eigenfunctions corresponding to Section 9-2 Sturm-Liouville Problems The Sturm-Liouville problem, the most famous eigenvalue problem of second-order ODE, is studied in this section due to its practical importance in engineering and science. 6
7 Self-adjoint form. A standard second-order linear differential equation written as (10) can be recast into the so-called self-adjoint form by multiplying to both sides of the equation (11) The above self-adjoint form serves as the generic form producing an infinite number of orthogonal functions, and is used to define the Sturm-Liouville problem. Table 9-1 Sturm-Liouville Problem and eigenfunction expansion differential equation: boundary conditions: at each side or at both sides (periodic) for, is a real continuous function for an interval,, eigenvalue: all the eigenvalues are real and of infinite number eigenfunction: two distinct eigenfunctions are orthogonal, eigenfunction expansion: Sturm-Liouville problem. For the interval equation is defined by the following self-adjoint form, the Sturm-Liouville (12) 7
8 where a parameter is called the eigenvalue. In the above, are all real-valued and continuous over an interval. Especially over an interval. The self-adjoint form in equation (12) and the boundary conditions as given in Table 9-1 consist of the boundary value problem called the Sturm- Liouville problem. Frequently, the conditions 1, 2, 3 in Table 9-1 are combined together as (13) where and are not zeros simultaneously. The condition 4 in Table 9-1 is related to the singular Sturm-Liouville problem, and the condition 5 to the periodic Sturm-Liouville problem. In contrast, the regular Sturm-Liouville problem is defined by the conditions 1, 2 and 3. Discretization of Sturm-Liouville equation. The following Sturm- Liouville equation (14) can be discretized for the internal grids such that (15) where the central difference approximation over the interval shown in Figure 9-3 is employed. By rearranging, we have (16) which is further standardized to represent a typical eigenvalue problem (17) 8
9 Then, by utilizing the tridiagonal form of matrix found by the inverse power method., eigenvalues can be easily Figure 9-3 Grid configuration for discretization of Sturm-Liouville equation (Example 1) For a Sturm-Liouville problem,,, find the fundamental eigenvalue and the fundamental eigenfunction by numerical treatment. The given ODE belongs to the Sturm-Liouville problem and is characterized by. Re-writing equation (16) and inserting gives (18) where 5 uniform subintervals are used, i.e.. And the boundary condition at can be discretized from equation (8-47) such that or (19) 9
10 Also, the boundary condition at is simply written as (20) These two boundary conditions are inserted into equation (18), and the following linear equations are obtained. (21) By defining the matrix (22) and by using the Fadeev-Leverrier method, the characteristic equation is derived as (refer to the note ) the solution of which gives (23) (24) Indeed, for the given matrix, the same result can be obtained by the QR iteration method also. However, the above method is frequently used only for small matrices, and does not utilize the tridiagonal form. Therefore, it is desirable to use the inverse power method and the eigenvalue shifting. In order to exploit the inverse power method, by putting and solving, we have 10
11 (25) The above-obtained and is solved to give is again substituted into the right-hand side of (26) This procedure is repeated to obtain the following result iter 1/lambda f[1] f[2] f[3] f[4] from which the converged solution is found to be (27) The fundamental eigenvalue shows an error of 3.6% compared with the exact value (see a note in ). And the fundamental eigenfunction is shown as below [the boundary value is determined by using equation (19)]. 11
12 The exact eigenvalue can be determined from the eigenfunction is determined to be where. M#math to find the exact eigenvalue is as follows. and #> double ftn(x) = tan(x*pi)+x; #> w =.bisect(0.51,0.99, ftn); lam = w^2; // eigenvalue w = lam = #> bvp.x(0,pi) ( {y''}+{%y}=0, [ {y}-{y'}=0, {y}=0 ] ).togo(u = y); ans = #> x = (0,pi).span(6).tr; #> [ x, u.tr.rowskip(-20), sin(w*(pi-x)), sin(w*(pi-x))/ ]; ans = [ ] [ ] [ ] [ ] [ ] [ ] Use of generalized differential operator. Here, the following generalized differential operator (28) or the linearized differential operator as in equation (8-65) 12
13 (29) is employed to express the eigenvalue problem as (30) Meanwhile, the differential operator is thought to be discretized to be (31) as was discussed in Chapter 8. Therefore, the eigenvalue problem (32) can be discretized as above. At this point, it is observed that the term in the right-hand side of equation (32) contributes nothing to the discretization of the differential operator. The eigenvalues are obtained in sequence of smaller magnitudes by using the inverse power method and eigenvalue shifting. For example, the Sturm-Liouville equation can be expanded and modified as (33) where the coefficients of the differential operator are defined to be (34) (Example 2) The Bessel equation transformed into the Sturm-Liouville problem as below is For the case of and, find the three smallest eigenvalues. The boundary condition is given as. Use the form for numerical solution. 13
14 Re-writing the equation in the form with, we have (35) By the way, it is actually unnecessary to specify the boundary condition at due to the nature of singular Sturm-Liouville problem. Indeed, in order to exploit the compact form, it is required to specify the boundary condition at. Therefore, an additional condition is imposed considering the symmetry (in reality, corresponds to the center of the cylinder). M#math is written as below #> double eiglam(s) { Lam = bvp.x[201](0,1) ( {y''} +1/x*{y'}-1/x^2*{y}+{%y} = 0, [ {y'}=0, {y}=0 ] ).shift(s);; return sqrt(lam); } #> eiglam(0); eiglam.lam = ans = #> eiglam(50); eiglam.lam = ans = #> eiglam(100); eiglam.lam = ans = then the result is It can be seen that use of 200 subintervals gives very accurate eigenvalues. The exact eigenvalues are obtained from the following M#math. #> solve.x (.J_1(x) = 0 ).span(0.01,11); ans = [ e-008 ] [ e-009 ] [ e-008 ] 14
15 Section 9-3 Eigenvalues on Stability and Convergence Rate In contrast to the situations where the main objective is to find eigenvalues, it also happens that eigenvalues are not required to be sought even though they play an underlying serious effect. This can be illustrated with stiff differential equations and iteration methods on matrix equations. Eigenvalues associated with initial value problems. It was already discussed in section 7-8 that the stiffness of the differential equations is determined from the eigenvalues of the Jacobian matrix given in equation (7-91). In solving initial value problems, eigenvalues of the Jacobian matrix have the following nature. greatest eigenvalue determines step sizes due to stability smallest eigenvalue determines the number of steps for integration Therefore, there occurs a severe difficulty if a ratio of the greatest to smallest eigenvalues differ by several orders of magnitudes. An example is presented below to discuss this behavior. (Example 3) When a numerical solution to the following coupled linear equations is obtained by the Runge-Kutta method with, discuss the numerical stability. First, the coupled linear equations are written in the following matrix form (36) and M#math is written as #> double y1e(x) = 2*exp(-2*x)-exp(-1000*x); #> double y2e(x) = -exp(-2*x)+exp(-1000*x); #> matrix odesol(double h) { 15
16 } x = (0,h).march(6); return ode[x] ( y1' = 996*y1+1996*y2, y2' = -998*y1-1998*y2, y1 = 1, y2 = 0 ).return( x, y1/y1e(x), y2/y2e(x) ); #> odesol(0.01); ans = [ #IND ] [ ] [ ] [ e e+008 ] [ e e+011 ] [ e e+014 ] #> odesol(0.05); ans = [ #IND ] [ ] [ ] [ ] [ ] [ e+006 ] #> odesol(0.001); ans = [ #IND ] [ ] [ ] [ ] [ ] [ ] From the above results, it can be seen that solutions with diverge, whereas a considerably accurate solution is obtained with the smallest step size of. The reason why these results are appearing is that the eigenvalues associated with the corresponding matrix (37) differ in two orders of magnitude, i.e.. The greatest eigenvalue limits the step size through the Euler's stability condition 16
17 (38) The convergence of numerical solutions depends on whether the stability condition is satisfied or not. The exact solution to (Example 3) is given below. A proper use of M#math is utilizing Spoke '.stiff(tol=0.01)' #> ode.x[1000](0,1) ( y1' = 996*y1+1996*y2, y2' = -998*y1-1998*y2, y1=1, y2=0 ).stiff.return(x,y1,y2).skip(10).plot; Eigenvalues occurring in iteration methods. When the Gauss-Seidel iteration method is used to solve coupled linear equations, the diagonal dominance should be determined a priori as was discussed in section 3-6. The convergence characteristics of iteration methods is greatly affected by distribution of the eigenvalues of associated matrix. (Example 4) Solve the following coupled linear equation by the Gauss-Seidel iteration method until the error falls below. The diagonal elements of the given matrix are greater than the off-diagonal elements, and therefore the diagonal dominance is satisfied. This implies that solution can be obtained by iteration. To apply the Gauss-Seidel iteration method, the given equations are modified as (39) 17
18 By initializing and performing the first iteration, we have (40) Subsequent iterations can be performed by using the following M#math #> (a,b,c,d) = (3.2,3,5,5.2);; (x,y) = (0,0);; a=3.2; #> for.i(1,1000) { } r1 = 7-b*y;; // a*x + b*y = 7 r2 = -2-c*x;; // c*x + d*y = -2 x = r1/a;; y = r2/d;; err = 7-a*x-b*y + -2-c*x-d*y ;; [ i, x, y, err ]; if( err < 1.e-7) break; // 354 iterations are required ans = [ ] ans = [ ] ans = [ ]... ans = [ e-007 ] ans = [ e-008 ] Then, a final solution can be obtained by a total of 354 iterations. Even though the given matrix is just a simple matrix, it requires more than 300 iterations, the reason of which stems from the underlying effect of the eigenvalue of the given matrix. This will be discussed in detail below. Role of eigenvalues in iteration methods. In order to examine the iteration method from the standpoint of eigenvalues, matrices, are introduced as (41) then, the iteration given in equation (39) can be written as 18 (42)
19 where is given as an initial guess at the start of iteration. Meanwhile, the error at the th iteration (43) is defined, and equation (42) is recast into (44) or By utilizing the above expression, the error at the (45) th iteration becomes (46) or (47) If the matrix is selected to be, then. This means that the solution is obtained by the direct method since iteration is not required. Otherwise the error at the -th itration is calculated recursively from equation (47) such that (48) where the matrix is a key factor determining the convergence rate. In the case of (Example 4), the matrix is (49) 19
20 the eigenvalues of which are and their absolute values are very close to the unity (see the note in ). If the absolute value of eigenvalues of matrix exceeds the unity, the error diverges as iterations continue. In other words, by letting be the absolute eigenvalue close to the unity, the magnitude of error diminishes as (50) Therefore, convergence is very slow when, whereas convergence is considerably accelerated when. For example, for the following case (51) the matrix is calculated to be (52) the eigenvalues of which are. As a result, the error diminishes rapidly as iterations continue. Actual computation shows that iter x y residue where the rate of residue decreasing is about a factor multiplied by, i.e. a decreasing rate of one-tenth per iteration is found. M#math for equation (52) #> I =.I(2); A = [ 30,3; 5,50 ]; B = [ 30; 0,50 ]; #> (I-A*B.inv).eig; 20
21 ans = [ 0.1 ] [ -0.1 ] Section 9-4 Exercises =========================== 1 For an interval of, the Sturm-Liouville problem is given as. Find the fundamental eigenvalue for the following homeogenous boundary conditions. (a) (b) (c) 2 For an interval of, the Sturm-Liouville equation known as the Bessel function is given as Find the fundamental eigenvalue for the following boundary conditions (a) (b) (c) 3 When an axial load of is applied, the curvature of the slender rod can be determined from where. The final result becomes an eigenvalue problem written as For the two cases shown in the accompanying figure, find the axial road where values are given as. 21
22 4 An axisymmetric vibration of a circular membrane is governed by the following differential equation For the case of, find the fundamental eigenvalue. 5 The most well-known Sturm-Liouville problems are as follows (a) Bessel equation (b) Legendre equation (c) Tchebyshev equation (d) Hermite equation 22
23 (e) Laguerre equation // // end of Lectures // Section 9-5 (optional) Higher-order ODE and Eigenvalue Problems Eigenvalue problems associated with higher-order ODEs can be solved in a similar way. In this section, eigenvalue related to the bending and vibration of a beam is considered which is characterized by the fourth-order ODEs. Bending of a beam. When a distributed load is applied on a beam, the vertical displacement of the beam is governed by the following fourth-order differential equation (53) where is the modulus of elasticity, is the moment of inertia of the cross section. Especially, the product of is termed the flexural rigidity. The boundary conditions for the beam shown in Figure 9-4 are (54) (55) By converting the ODE in equation (53) into the form, we have (56) 23
24 where in general. Figure 9-4 Bending of a beam with a fixed end However, the boundary condition at takes a different form other than equation (55), depending on the state in which the right end of the beam is located. (Example 5) Find the displacement of a beam when a distributed load is applied to the beam. Here, boundary conditions are given as, and, (units are omitted for convenience). Find a numerical solution with 100 subintervals. Since is a constant, solution can be found easily by using equation (53) with the following program #> double ye(x) { L = 8; Lpi = L/pi; kei = 0.001; a2 = 0.5*L/pi^3; return kei*( Lpi^3*(Lpi*sin(x/Lpi)-x) + a2*x^2*(3*l-x) ); } #> kei = 0.001; L = 8; #> bvp.x(0,l) ( {y''''}-kei*sin(pi*x/l) = 0, [ {y}=0, {y'}=0, {y}=0, {y''}=0 ] ).return( x, ye(x), y/ye(x) ).rowskip(-10); A use of 100 subintervals produces a fairly exact solution as below. 24
25 kei = L = 8 ans = [ #INF ] [ ] [ ] [ ] [ ] [ ] [ ] [ ] [ ] [ ] [ #INF ] The exact solution is given as Vibration of beam. Now, the dynamic behavior related to the vibration of a beam is expressed by the following partial differential equation (57) where is the distributed load on the beam, is the mass per length. Consider the case where the displacement is given as from which equation (57) becomes ( ) (58) (59) If to be were constants, the exact solution to the above equation can be found 25
26 (60) from the boundary conditions [equations (54) and (55)]. (Example 6) For the case of, find the smallest three eigenvalues for the eigenvalue problem given in equation (59). Numerical solutions should be obtained with from 10 to 100 subintervals by adding 10 subintervals after each run. Find the effect of grid density on the eigenvalue obtained numerically. Exact solutions are found to be, and (see the note below). Since is a constant, the form of differential operator is very simple. M#math is written as #> L = 8; EIm = 1000; s = 0; // change s for other eigenvalues #> row = []; #> for.nx(10,100, 10) { Lam = bvp.x[nx+1](0,l) ( -EIm*{f''''}+{%f} = 0, [ {f} = 0, {f'} = 0, {f''}=0, {f'''}=0 ] ).shift(s); row _= [ nx, sqrt(lam)];; // vertical concatenation } #> row; row = [ ] [ ] [ ] [ ] [ ] [ ] [ ] [ ] [ ] [ ] Comparing with the exact solution, the result corresponding to 100 subintervals gives a consistency up to 4 significant digits. Finding other eigenvalues can be carried out by changing the shifintg value in the above M#math code. 26
27 // // end of file //
Chapter 8 Boundary Value Problems
"msharpmath" revised on 2014.09.01 Numerical Analysis Using M#math Chapter 8 Boundary Value Problems 8-1 Definition of Boundary Value Problems 8-2 Shooting Method 8-3 Second-order Finite Difference Method
More informationPARTIAL DIFFERENTIAL EQUATIONS and BOUNDARY VALUE PROBLEMS
PARTIAL DIFFERENTIAL EQUATIONS and BOUNDARY VALUE PROBLEMS NAKHLE H. ASMAR University of Missouri PRENTICE HALL, Upper Saddle River, New Jersey 07458 Contents Preface vii A Preview of Applications and
More informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
More informationCOPYRIGHTED MATERIAL. Index
Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,
More informationLecture IX. Definition 1 A non-singular Sturm 1 -Liouville 2 problem consists of a second order linear differential equation of the form.
Lecture IX Abstract When solving PDEs it is often necessary to represent the solution in terms of a series of orthogonal functions. One way to obtain an orthogonal family of functions is by solving a particular
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationPhysics 6303 Lecture 15 October 10, Reminder of general solution in 3-dimensional cylindrical coordinates. sinh. sin
Physics 6303 Lecture 15 October 10, 2018 LAST TIME: Spherical harmonics and Bessel functions Reminder of general solution in 3-dimensional cylindrical coordinates,, sin sinh cos cosh, sin sin cos cos,
More informationPreface. 2 Linear Equations and Eigenvalue Problem 22
Contents Preface xv 1 Errors in Computation 1 1.1 Introduction 1 1.2 Floating Point Representation of Number 1 1.3 Binary Numbers 2 1.3.1 Binary number representation in computer 3 1.4 Significant Digits
More informationMatrix Eigenvalue Problems
"msharpmath" revised 2014.09.01 Numerical Analysis Using m#math Chapter 4 Matrix Eigenvalue Problems www.msharpmath.com * Translator or contributor ** *[101] Msharpmath Inc. **[your ID] Your affiliation
More informationStructural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).
Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationChapter 2 Solution of Nonlinear Equations
"msharpmath" revised on 2014.09.01 Numerical Analysis Using M#math Chapter 2 Solution of Nonlinear Equations www.msharpmath.com * Translator or contributor ** *[101] Msharpmath Inc. **[your ID] Your affiliation
More informationMechanical Design in Optical Engineering
OPTI Buckling Buckling and Stability: As we learned in the previous lectures, structures may fail in a variety of ways, depending on the materials, load and support conditions. We had two primary concerns:
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationEigenvalues of Trusses and Beams Using the Accurate Element Method
Eigenvalues of russes and Beams Using the Accurate Element Method Maty Blumenfeld Department of Strength of Materials Universitatea Politehnica Bucharest, Romania Paul Cizmas Department of Aerospace Engineering
More informationAdvanced Vibrations. Distributed-Parameter Systems: Exact Solutions (Lecture 10) By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Exact Solutions (Lecture 10) By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Exact Solutions Relation between Discrete and Distributed
More informationAdvanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Approximate Methods Lecture 20 By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Approximate Methods Rayleigh's Principle The Rayleigh-Ritz
More informationAND NONLINEAR SCIENCE SERIES. Partial Differential. Equations with MATLAB. Matthew P. Coleman. CRC Press J Taylor & Francis Croup
CHAPMAN & HALL/CRC APPLIED MATHEMATICS AND NONLINEAR SCIENCE SERIES An Introduction to Partial Differential Equations with MATLAB Second Edition Matthew P Coleman Fairfield University Connecticut, USA»C)
More informationFOURIER SERIES, TRANSFORMS, AND BOUNDARY VALUE PROBLEMS
fc FOURIER SERIES, TRANSFORMS, AND BOUNDARY VALUE PROBLEMS Second Edition J. RAY HANNA Professor Emeritus University of Wyoming Laramie, Wyoming JOHN H. ROWLAND Department of Mathematics and Department
More informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
More informationPHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.
PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional
More informationBOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES
1 BOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES 1.1 Separable Partial Differential Equations 1. Classical PDEs and Boundary-Value Problems 1.3 Heat Equation 1.4 Wave Equation 1.5 Laplace s Equation
More informationModule 10: Finite Difference Methods for Boundary Value Problems Lecture 42: Special Boundary Value Problems. The Lecture Contains:
The Lecture Contains: Finite Difference Method Boundary conditions of the second and third kind Solution of the difference scheme: Linear case Solution of the difference scheme: nonlinear case Problems
More informationChapter 3 Higher Order Linear ODEs
Chapter 3 Higher Order Linear ODEs Advanced Engineering Mathematics Wei-Ta Chu National Chung Cheng University wtchu@cs.ccu.edu.tw 1 2 3.1 Homogeneous Linear ODEs 3 Homogeneous Linear ODEs An ODE is of
More informationNumerical Methods I Non-Square and Sparse Linear Systems
Numerical Methods I Non-Square and Sparse Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 25th, 2014 A. Donev (Courant
More informationLevel 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method
9210-203 Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method You should have the following for this examination one answer book No additional data is attached
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationNumerical Analysis. A Comprehensive Introduction. H. R. Schwarz University of Zürich Switzerland. with a contribution by
Numerical Analysis A Comprehensive Introduction H. R. Schwarz University of Zürich Switzerland with a contribution by J. Waldvogel Swiss Federal Institute of Technology, Zürich JOHN WILEY & SONS Chichester
More informationMaple in Differential Equations
Maple in Differential Equations and Boundary Value Problems by H. Pleym Maple Worksheets Supplementing Edwards and Penney Differential Equations and Boundary Value Problems - Computing and Modeling Preface
More informationACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY THE SUPERPOSITION METHOD
Journal of Sound and Vibration (1999) 219(2), 265 277 Article No. jsvi.1998.1874, available online at http://www.idealibrary.com.on ACCURATE FREE VIBRATION ANALYSIS OF POINT SUPPORTED MINDLIN PLATES BY
More informationTyn Myint-U Lokenath Debnath. Linear Partial Differential Equations for Scientists and Engineers. Fourth Edition. Birkhauser Boston Basel Berlin
Tyn Myint-U Lokenath Debnath Linear Partial Differential Equations for Scientists and Engineers Fourth Edition Birkhauser Boston Basel Berlin Preface to the Fourth Edition Preface to the Third Edition
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationLecture 10: Fourier Sine Series
Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. Not to be copied, used, or revised without eplicit written permission from the copyright owner. ecture : Fourier Sine Series
More informationReview for Exam 2 Ben Wang and Mark Styczynski
Review for Exam Ben Wang and Mark Styczynski This is a rough approximation of what we went over in the review session. This is actually more detailed in portions than what we went over. Also, please note
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 Lecture Notes 8 - PDEs Page 8.01 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationChapter 4 Analysis of a cantilever
Chapter 4 Analysis of a cantilever Before a complex structure is studied performing a seismic analysis, the behaviour of simpler ones should be fully understood. To achieve this knowledge we will start
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationAPPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems
APPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems Fourth Edition Richard Haberman Department of Mathematics Southern Methodist University PEARSON Prentice Hall PEARSON
More informationADVANCED ENGINEERING MATHEMATICS MATLAB
ADVANCED ENGINEERING MATHEMATICS WITH MATLAB THIRD EDITION Dean G. Duffy Contents Dedication Contents Acknowledgments Author Introduction List of Definitions Chapter 1: Complex Variables 1.1 Complex Numbers
More informationEFFECTIVE MODAL MASS & MODAL PARTICIPATION FACTORS Revision F
EFFECTIVE MODA MASS & MODA PARTICIPATION FACTORS Revision F By Tom Irvine Email: tomirvine@aol.com March 9, 1 Introduction The effective modal mass provides a method for judging the significance of a vibration
More informationADVANCED ENGINEERING MATHEMATICS
ADVANCED ENGINEERING MATHEMATICS DENNIS G. ZILL Loyola Marymount University MICHAEL R. CULLEN Loyola Marymount University PWS-KENT O I^7 3 PUBLISHING COMPANY E 9 U Boston CONTENTS Preface xiii Parti ORDINARY
More informationApplied Linear Algebra
Applied Linear Algebra Peter J. Olver School of Mathematics University of Minnesota Minneapolis, MN 55455 olver@math.umn.edu http://www.math.umn.edu/ olver Chehrzad Shakiban Department of Mathematics University
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationFIXED BEAMS IN BENDING
FIXED BEAMS IN BENDING INTRODUCTION Fixed or built-in beams are commonly used in building construction because they possess high rigidity in comparison to simply supported beams. When a simply supported
More informationENGINEERING MATHEMATICS I. CODE: 10 MAT 11 IA Marks: 25 Hrs/Week: 04 Exam Hrs: 03 PART-A
ENGINEERING MATHEMATICS I CODE: 10 MAT 11 IA Marks: 25 Hrs/Week: 04 Exam Hrs: 03 Total Hrs: 52 Exam Marks:100 PART-A Unit-I: DIFFERENTIAL CALCULUS - 1 Determination of n th derivative of standard functions-illustrative
More informationDeflection profile analysis of beams on two-parameter elastic subgrade
1(213) 263 282 Deflection profile analysis of beams on two-parameter elastic subgrade Abstract A procedure involving spectral Galerkin and integral transformation methods has been developed and applied
More informationUpon successful completion of MATH 220, the student will be able to:
MATH 220 Matrices Upon successful completion of MATH 220, the student will be able to: 1. Identify a system of linear equations (or linear system) and describe its solution set 2. Write down the coefficient
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 17
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil engineering Indian Institute of Technology, Madras Module - 01 Lecture - 17 In the last class, we were looking at this general one-dimensional
More informationAPPLIED PARTIAL DIFFERENTIAL EQUATIONS
APPLIED PARTIAL DIFFERENTIAL EQUATIONS AN I N T R O D U C T I O N ALAN JEFFREY University of Newcastle-upon-Tyne ACADEMIC PRESS An imprint of Elsevier Science Amsterdam Boston London New York Oxford Paris
More information9.1 Introduction to bifurcation of equilibrium and structural
Module 9 Stability and Buckling Readings: BC Ch 14 earning Objectives Understand the basic concept of structural instability and bifurcation of equilibrium. Derive the basic buckling load of beams subject
More informationNumerical Analysis: Solutions of System of. Linear Equation. Natasha S. Sharma, PhD
Mathematical Question we are interested in answering numerically How to solve the following linear system for x Ax = b? where A is an n n invertible matrix and b is vector of length n. Notation: x denote
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 11
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 11 Last class, what we did is, we looked at a method called superposition
More informationAdaptation of the Lanczos Algorithm for the Solution of Buckling Eigenvalue Problems
Original Article Adaptation of the Lanczos Algorithm for the Solution of Buckling Eigenvalue Problems Abstract An adaptation of the conventional Lanczos algorithm is proposed to solve the general symmetric
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 ecture Notes 8 - PDEs Page 8.0 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationPhysics 6303 Lecture 11 September 24, LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation
Physics 6303 Lecture September 24, 208 LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation, l l l l l l. Consider problems that are no axisymmetric; i.e., the potential depends
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationBoundary Value Problems - Solving 3-D Finite-Difference problems Jacob White
Introduction to Simulation - Lecture 2 Boundary Value Problems - Solving 3-D Finite-Difference problems Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Reminder about
More informationChapter Two Elements of Linear Algebra
Chapter Two Elements of Linear Algebra Previously, in chapter one, we have considered single first order differential equations involving a single unknown function. In the next chapter we will begin to
More informationChapter 7 Initial Value Problems
"msharpmath" revised on 2014.09.01 Numerical Analysis Using M#math Chapter 7 Initial Value Problems www.msharpmath.com * Translator or contributor ** *[101] Msharpmath Inc. **[your ID] Your affiliation
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationON NUMERICAL METHODS FOR ELLIPTIC TRANSMISSION/BOUNDARY VALUE PROBLEMS
ON NUMERICAL METHODS FOR ELLIPTIC TRANSMISSION/BOUNDARY VALUE PROBLEMS ANIRBAN ROY 349 LAFAYETTE DR, BETHLEHEM, PA, 182 Project Supervisors: Anna L. Mazzucato and Victor Nistor Department of Mathematics,
More informationSolution of Differential Equation by Finite Difference Method
NUMERICAL ANALYSIS University of Babylon/ College of Engineering/ Mechanical Engineering Dep. Lecturer : Dr. Rafel Hekmat Class : 3 rd B.Sc Solution of Differential Equation by Finite Difference Method
More informationIn this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation.
Eigen Function Expansion and Applications. In this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation. a/ The theory. b/ Example: Solving the Euler equation in two ways.
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationTopics for the Qualifying Examination
Topics for the Qualifying Examination Quantum Mechanics I and II 1. Quantum kinematics and dynamics 1.1 Postulates of Quantum Mechanics. 1.2 Configuration space vs. Hilbert space, wave function vs. state
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 9, February 8, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 9, February 8, 2006 The Harmonic Oscillator Consider a diatomic molecule. Such a molecule
More informationCHAPTER 8: Matrices and Determinants
(Exercises for Chapter 8: Matrices and Determinants) E.8.1 CHAPTER 8: Matrices and Determinants (A) means refer to Part A, (B) means refer to Part B, etc. Most of these exercises can be done without a
More information3. Overview of MSC/NASTRAN
3. Overview of MSC/NASTRAN MSC/NASTRAN is a general purpose finite element analysis program used in the field of static, dynamic, nonlinear, thermal, and optimization and is a FORTRAN program containing
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationNumerical simulation of the coil spring and investigation the impact of tension and compression to the spring natural frequencies
Numerical simulation of the coil spring and investigation the impact of tension and compression to the spring natural frequencies F. D. Sorokin 1, Zhou Su 2 Bauman Moscow State Technical University, Moscow,
More informationSystems of Algebraic Equations and Systems of Differential Equations
Systems of Algebraic Equations and Systems of Differential Equations Topics: 2 by 2 systems of linear equations Matrix expression; Ax = b Solving 2 by 2 homogeneous systems Functions defined on matrices
More informationNumerical Methods for Engineers. and Scientists. Applications using MATLAB. An Introduction with. Vish- Subramaniam. Third Edition. Amos Gilat.
Numerical Methods for Engineers An Introduction with and Scientists Applications using MATLAB Third Edition Amos Gilat Vish- Subramaniam Department of Mechanical Engineering The Ohio State University Wiley
More informationCh 10.1: Two Point Boundary Value Problems
Ch 10.1: Two Point Boundary Value Problems In many important physical problems there are two or more independent variables, so the corresponding mathematical models involve partial differential equations.
More informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationLecture 4.6: Some special orthogonal functions
Lecture 4.6: Some special orthogonal functions Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4340, Advanced Engineering Mathematics
More information7.5 Elastic Buckling Columns and Buckling
7.5 Elastic Buckling The initial theory of the buckling of columns was worked out by Euler in 1757, a nice example of a theory preceding the application, the application mainly being for the later invented
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationGyroscopic matrixes of the straight beams and the discs
Titre : Matrice gyroscopique des poutres droites et des di[...] Date : 29/05/2013 Page : 1/12 Gyroscopic matrixes of the straight beams and the discs Summarized: This document presents the formulation
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 informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationMODULE C: COMPRESSION MEMBERS
MODULE C: COMPRESSION MEMBERS This module of CIE 428 covers the following subjects Column theory Column design per AISC Effective length Torsional and flexural-torsional buckling Built-up members READING:
More information9. Iterative Methods for Large Linear Systems
EE507 - Computational Techniques for EE Jitkomut Songsiri 9. Iterative Methods for Large Linear Systems introduction splitting method Jacobi method Gauss-Seidel method successive overrelaxation (SOR) 9-1
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
More informationInternational Journal of Advanced Engineering Technology E-ISSN
Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,
More informationNotes for CS542G (Iterative Solvers for Linear Systems)
Notes for CS542G (Iterative Solvers for Linear Systems) Robert Bridson November 20, 2007 1 The Basics We re now looking at efficient ways to solve the linear system of equations Ax = b where in this course,
More informationIntroduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series
CHAPTER 5 Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series We start with some introductory examples. 5.. Cauchy s equation The homogeneous Euler-Cauchy equation (Leonhard
More informationProject: Vibration of Diatomic Molecules
Project: Vibration of Diatomic Molecules Objective: Find the vibration energies and the corresponding vibrational wavefunctions of diatomic molecules H 2 and I 2 using the Morse potential. equired Numerical
More informationNext topics: Solving systems of linear equations
Next topics: Solving systems of linear equations 1 Gaussian elimination (today) 2 Gaussian elimination with partial pivoting (Week 9) 3 The method of LU-decomposition (Week 10) 4 Iterative techniques:
More information