Numerical+Methods 源代码

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

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 U=forwdif(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 e…

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

function [T,Z]=rks4(F,a,b,Za, M) %Input - F is the system input as a string 'F' % - a and b are the endpoints of the i…

function T=rctrap(f,a,b,n) %Input - f is the integrand input as a string 'f' % - a and b are upper and lower limits of int…

function [R,quad,err,h]=romber(f,a,b,n,tol) %Input - f is the integrand input as a string 'f' % - a and b are upper and lo…

function T=house (A) %Input - A is an nxn symmetric matrix %Output - T is a tridiagonal matrix % NUMERICAL METHODS: MATLAB Progr…

function [A,df]=diffnew(X,Y) %Input - X is the 1xn abscissa vector % - Y is the 1xn ordinate vector %Output - A is the …