Numerical+Methods 源代码

function [lambda,V]= power1(A,X,epsilon,max1) %Input - A is an nxn matrix % - X is the nx1 starting vector % - epsilon …

function D=qr1(A,epsilon) %Input - A is a symmetric tridiagona nxn matrix % - epsilon is the tolerance %Output - D is the ve…

function T4=taylor(df,a,b,ya,M) %Input - df=[y' y'' y''' y'''']entered as a string 'df' % where y'=f(t,y) % - a…

function quad=gauss(f,a,b,A,W) %Input - f is the integrand input as a string 'f' % - a and b upper and lower limits of i…

function M=milne(f,T,Y) %Input - f is the function entered as a string 'f' % - T is the vector of abscissas % - Y…

function L=ls(F1,F2,a,b,alpha,beta,M) %Input - F1 and F2 are the systems of first-order equations % representing the …

function L=linsht(F1,F2,a,b,alpha,beta,M) %Input - F1 and F2 are the systems of first-order equations % representing …

function X=trisys (A,D,C,B) %Input - A is the sub diagonal of the coefficient matrix % - D is the main diagonal of the coeff…

function [A,B]=tpcoeff(X,Y,M) %Input - X is a vector of equally spaced abscisssas in [-pi, pi] % - Y is a vector of or…

NUMERICAL METHODS: Matlab Programs (c) 2004 by John H. Mathews and Kurtis D. Fink Complementary Software to accompany the textbo…