数值方法matlab版源代码

John H.Mathews及Kurtis D.Fink著<>(第四版）

function A=abm（fTY）

%Input   - f is the function 
%           - T is the vector of abscissas
%           - Y is the vector of ordinates
%Remark.  The first four coordinates of T and Y must
%               have starting values obtained with RK4
%Output - A=[T‘ Y‘] where T is the vector of abscissas and
%              Y is the vector of ordinates

% If f is an M-file function call A=abm（@fTY）.
% If f is an anonymous function call A=abm（fTY）.

%  NUMERICAL METHODS: Matlab Programs
% （c） 2004 by John H. Mathews and Kurtis D. Fink
%  Complementary Software to accompany the textbook:
%  NUMERICAL METHODS: Using Matlab Fourth Edition
%  ISBN: 0-13-065248-2
%  Prentice-Hall Pub. Inc.
%  One Lake Street
%  Upper Saddle River NJ 07458

n=length（T）;
if n<5return end;

F=zeros（14）;
F=f（T（1:4）Y（1:4））;
h=T（2）-T（1）;

for k=4:n-1
   %Predictor
   p=Y（k）+（h/24）*（F*[-9 37 -59 55]‘）;
   T（k+1）=T（1）+h*k;
   F=[F（2） F（3） F（4） f（T（k+1）p）];
   %Corrector
   Y（k+1）=Y（k）+（h/24）*（F*[1 -5 19 9]‘）;
   F（4）=f（T（k+1）Y（k+1））;
end

A=[T‘ Y‘];

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1044  2006-06-30 14:57  数值方法（matlab版）\NumericalMethods7.1\abm.m

     文件       1937  2006-06-25 11:47  数值方法（matlab版）\NumericalMethods7.1\adapt.m

     文件        919  2006-06-29 08:57  数值方法（matlab版）\NumericalMethods7.1\approot.m

     文件        625  2006-06-24 06:28  数值方法（matlab版）\NumericalMethods7.1\backsub.m

     文件       1007  2006-06-29 08:32  数值方法（matlab版）\NumericalMethods7.1\bisect.m

     文件       1020  2006-06-29 12:41  数值方法（matlab版）\NumericalMethods7.1\cheby.m

     文件       1309  2006-06-28 11:26  数值方法（matlab版）\NumericalMethods7.1\crnich.m

     文件       1259  2004-01-19 05:04  数值方法（matlab版）\NumericalMethods7.1\csfit.m

     文件       1248  2006-06-25 07:30  数值方法（matlab版）\NumericalMethods7.1\diffext.m

     文件       1260  2006-06-25 07:28  数值方法（matlab版）\NumericalMethods7.1\difflim.m

     文件        774  2006-06-25 07:38  数值方法（matlab版）\NumericalMethods7.1\diffnew.m

     文件       1600  2006-06-28 11:58  数值方法（matlab版）\NumericalMethods7.1\dirich.m

     文件        891  2006-06-26 11:01  数值方法（matlab版）\NumericalMethods7.1\euler.m

     文件        145  2004-01-19 05:07  数值方法（matlab版）\NumericalMethods7.1\fib.m

     文件       1801  2006-06-25 15:02  数值方法（matlab版）\NumericalMethods7.1\fibonacci.m

     文件       1420  2006-06-30 16:20  数值方法（matlab版）\NumericalMethods7.1\findiff.m

     文件       1170  2006-06-30 16:43  数值方法（matlab版）\NumericalMethods7.1\finedif.m

     文件       1079  2006-06-29 08:35  数值方法（matlab版）\NumericalMethods7.1\fixpt.m

     文件       1061  2006-06-30 18:21  数值方法（matlab版）\NumericalMethods7.1\forwdif.m

     文件        843  2006-06-25 11:52  数值方法（matlab版）\NumericalMethods7.1\gauss.m

     文件       1584  2006-06-25 14:58  数值方法（matlab版）\NumericalMethods7.1\golden.m

     文件       2936  2006-06-26 06:49  数值方法（matlab版）\NumericalMethods7.1\grads.m

     文件       1123  2006-06-24 07:01  数值方法（matlab版）\NumericalMethods7.1\gseid.m

     文件       1214  2006-06-27 05:47  数值方法（matlab版）\NumericalMethods7.1\hamming.m

     文件        935  2006-06-30 11:35  数值方法（matlab版）\NumericalMethods7.1\heun.m

     文件        913  2004-01-19 05:08  数值方法（matlab版）\NumericalMethods7.1\house.m

     文件       1127  2006-06-28 15:48  数值方法（matlab版）\NumericalMethods7.1\invpow.m

     文件        907  2006-06-24 06:54  数值方法（matlab版）\NumericalMethods7.1\jacobi.m

     文件       1092  2006-06-28 07:24  数值方法（matlab版）\NumericalMethods7.1\jacobi1.m

     文件        937  2006-06-25 04:40  数值方法（matlab版）\NumericalMethods7.1\lagran.m

............此处省略34个文件信息