Numerical+Methods 源代码

function R = approot (X,epsilon) % Input - f is object function saved as an M-file named f.m % - X is the vector of absciss…

function X=jacobi(A,B,P,delta, max1) % Input - A is an N x N nonsingular matrix % - B is an N x 1 matrix % - P is …

function Z=srule(f,a0,b0,tol0) %Input - f is the integrand input as a string 'f' % - a0 and b0 are upper and lower limit…

function [c,err,yc]=bisect(f,a,b,delta) %Input - f is the function input as a string 'f' % - a and b are the left and right e…

function [lambda,V]= power(A,X,epsilon,max1) % NUMERICAL METHODS: Matlab Programs % (c) 2004 by John H. Mathews and Kurtis D. Fi…

function H=heun(f,a,b,ya,M) %Input - f is the function entered as a string 'f' % - a and b are the left and right endpoi…

function C = lspoly(X,Y,M) %Input - X is the 1xn abscissa vector % - Y is the 1xn ordinate vector % - M is the de…

function D=qr2(A,epsilon) %Input - A is a symmetric tridiagonal nxn matrix % - epsilon is the tolerance %Output - D is the n…

function [L,n]=difflim(f,x,toler) %Input - f is the function input as a string 'f' % - x is the differentiation point % …

function [C,X,Y]=cheby(fun,n,a,b) %Input - fun is the string function to be approximated % - N is the degree of the Cheb…