10.07.2015 Views

Numerical Methods Library for OCTAVE

Numerical Methods Library for OCTAVE

Numerical Methods Library for OCTAVE

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

N number of subdivisions.ReturnsRES result of integration.3.5.3 Newton-Cotes’ ruleRES = INTE_NEWTCOT(FUN,A,B,N) computes an approximation of the integral ofthe function FUN via the Newton-Cotes method (using N equispaced intervals).FUN accepts real scalar input x and returns a real scalar value. FUN canalso be an inline object.Parameters FUN integrated function.of subdivisions.A,B FUN is integrated on [A,B]. N numberReturns RES result of integration.3.6 Eigenvalue problemsThe standard algebraic eigenvalue problem is as follows: Given an n ×n matrix A, find a scalarλ and a nonzero vector x such thatAx = λx. (8)Such a scalar λ is called an eigenvalue, and x is a corresponding eigenvector. The set of all theeigenvalues of a matrix A, denoted by λ(A), is called the spectrum of A.An eigenvector of a matrix determines a direction in which the effect of the matrix is particularlysimple: The matrix expands or shrinks any vector lying in that direction by a scalarmultiple, and the expansion or contraction factor is given by the corresponding eigenvalue λ.Thus, eigenvalues and eigenvectors provide a means of understanding the complicated behaviorof a general linear trans<strong>for</strong>mation by decomposing it into simpler actions.Eigenvalue problems occur in many areas of science and engineering. For example, thenatural modes and frequencies of vibration of a structure are determined by the eigenvectorsand eigenvalues of an appropriate matrix. The stability of the structure is determined by thelocations of the eigenvalues, and thus their computation is of critical interest. We will also seelater in this book that eigenvalues can be very useful in analyzing numerical methods, such asthe convergence analysis of iterative methods <strong>for</strong> solving systems of algebraic equations, andthe stability analysis of methods <strong>for</strong> solving systems of differential equations.3.6.1 Power iteration[LAMBDA, V, NBIT] = EIG_POWER(A, X0, ITMAX, TOL) computes dominant eigenvalueand associated eigenvector of A with power iteration method. If EIG_POWERfails to converge after the maximum number of iterations or halts <strong>for</strong> anyreason, a message is displayed.Parameters15

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!