numerical methods 源代码

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 X=gseid(A,B,P,delta, max1) % Input - A is an N x N nonsingular matrix % - B is an N x 1 matrix % …

function E=euler(f,a,b,ya,M) %Input - y'=f is the function % - a and b are the left and right endpoints % …

function [lambda,V]= power1(A,X,epsilon,max1) %Input - A is an nxn matrix % - X is the nx1 starting vector % …

function T4=taylor(df,a,b,ya,M) %Input - df=[y' y'' y''' y''''] where y'=f(t,y) % - a and b are the left and right …

function quad=gauss(f,a,b,A,W) %Input - f is the integrand % - a and b upper and lower limits of integration % …