Mathematical+Methods 源代码

function U=dirich(f1,f2,f3,f4,a,b,h,tol,max1) %Input - f1,f2,f3,f4 are boundary functions input as strings % - a and b right…

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

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…