Numerical+Methods 源代码

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…

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 U = finedif(f,g,a,b,c,n,m) %Input - f=u(x,0) as a string 'f' % - g=ut(x,0) as a string 'g' % - a and b right e…

function A=abm(f,T,Y) %Input - f is the function entered as a string 'f' % - T is the vector of abscissas % - Y i…

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 [C,L]=lagran(X,Y) %Input - X is a vector that contains a list of abscissas % - Y is a vector that contains a list…

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

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