Numerical+Methods 源代码

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

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 X=backsub(A,B) %Input - A is an n x n upper-triangular nonsingular matrix % - B is an n x 1 matrix %Output - X is …

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…

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

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

function X=gseid(A,B,P,delta, max1) % Input - A is an N x N nonsingular matrix % - B is an N x 1 matrix % - P is a…

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 [lambda,V]= power1(A,X,epsilon,max1) %Input - A is an nxn matrix % - X is the nx1 starting vector % - epsilon …

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