Numerical+Methods 源代码

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…

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 containing th…

function [P,iter]= seidel(G,P,delta, max1) %Input - G is the nonlinear system saved in the M-file G.m % - P is the initi…

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 D=qr1(A,epsilon) %Input - A is a symmetric tridiagona nxn matrix % - epsilon is the tolerance %Output - D is the ve…

function T4=taylor(df,a,b,ya,M) %Input - df=[y' y'' y''' y'''']entered as a string 'df' % where y'=f(t,y) % - a…

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

function L=ls(F1,F2,a,b,alpha,beta,M) %Input - F1 and F2 are the systems of first-order equations % representing the …