numerical methods 源代码

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

function [L,n]=difflim(f,x,toler) %Input - f is the function % - x is the differentiation point % - t…

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

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

function U = finedif(f,g,a,b,c,n,m) %Input - f=u(x,0) % - g=ut(x,0) % - a and b right endpoints of [0,…

function A=abm(f,T,Y) %Input - f is the function % - T is the vector of abscissas % - Y is the vector of o…

function [R,quad,err,h]=romber(f,a,b,n,tol) %Input - f is the integrand % - a and b are upper and lower limits of …

function [C,L]=lagran(X,Y) %Input - X is a vector that contains a list of abscissas % - Y is a vector that contains…

function [C,D]=newpoly(X,Y) %Input - X is a vector that contains a list of abscissas % - Y is a vector that conta…

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