numerical methods 源代码

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 …

function s=simprl(f,a,b,M) %Input - f is the integrand % - a and b are upper and lower limits of integration % …

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 m…

function R=rk4(f,a,b,ya,M) %Input - f is the function % - a and b are the left and right endpoints % - …

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

function s=traprl(f,a,b,M) %Input - f is the integrand % - a and b are upper and lower limits of integration %…

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

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 T=house (A) %Input - A is an nxn symmetric matrix %Output - T is a tridiagonal matrix % NUMERICAL METHODS: Matlab Prog…