numerical methods 源代码

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

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

function U=crnich(f,c1,c2,a,b,c,n,m) %Input - f=u(x,0) % - c1=u(0,t) and c2=u(a,t) % - a and b right end…

function R = approot (f,X,epsilon) % Input - f is object function % - X is the vector of abscissas % …

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

function Z=srule(f,a0,b0,tol0) %Input - f is the integrand % - a0 and b0 are upper and lower limits of integrati…

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

function H=heun(f,a,b,ya,M) %Input - y'= f is the function % - a and b are the left and right endpoints % …

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

function U=forwdif(f,c1,c2,a,b,c,n,m) %Input - f=u(x,0) % - c1=u(0,t) and c2=u(a,t) % - a and b right endp…