Numerical+Methods 源代码

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 …

function X = uptrbk(A,B) %Input - A is an N x N nonsingular matrix % - B is an N x 1 matrix %Output - X is an N x 1 matrix …

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