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??數(shù)值方法（matlab版）\NumericalMethods7.1\abm.m

?????文件???????1937??2006-06-25?11:47??數(shù)值方法（matlab版）\NumericalMethods7.1\adapt.m

?????文件????????919??2006-06-29?08:57??數(shù)值方法（matlab版）\NumericalMethods7.1\approot.m

?????文件????????625??2006-06-24?06:28??數(shù)值方法（matlab版）\NumericalMethods7.1\backsub.m

?????文件???????1007??2006-06-29?08:32??數(shù)值方法（matlab版）\NumericalMethods7.1\bisect.m

?????文件???????1020??2006-06-29?12:41??數(shù)值方法（matlab版）\NumericalMethods7.1\cheby.m

?????文件???????1309??2006-06-28?11:26??數(shù)值方法（matlab版）\NumericalMethods7.1\crnich.m

?????文件???????1259??2004-01-19?05:04??數(shù)值方法（matlab版）\NumericalMethods7.1\csfit.m

?????文件???????1248??2006-06-25?07:30??數(shù)值方法（matlab版）\NumericalMethods7.1\diffext.m

?????文件???????1260??2006-06-25?07:28??數(shù)值方法（matlab版）\NumericalMethods7.1\difflim.m

?????文件????????774??2006-06-25?07:38??數(shù)值方法（matlab版）\NumericalMethods7.1\diffnew.m

?????文件???????1600??2006-06-28?11:58??數(shù)值方法（matlab版）\NumericalMethods7.1\dirich.m

?????文件????????891??2006-06-26?11:01??數(shù)值方法（matlab版）\NumericalMethods7.1\euler.m

?????文件????????145??2004-01-19?05:07??數(shù)值方法（matlab版）\NumericalMethods7.1\fib.m

?????文件???????1801??2006-06-25?15:02??數(shù)值方法（matlab版）\NumericalMethods7.1\fibonacci.m

?????文件???????1420??2006-06-30?16:20??數(shù)值方法（matlab版）\NumericalMethods7.1\findiff.m

?????文件???????1170??2006-06-30?16:43??數(shù)值方法（matlab版）\NumericalMethods7.1\finedif.m

?????文件???????1079??2006-06-29?08:35??數(shù)值方法（matlab版）\NumericalMethods7.1\fixpt.m

?????文件???????1061??2006-06-30?18:21??數(shù)值方法（matlab版）\NumericalMethods7.1\forwdif.m

?????文件????????843??2006-06-25?11:52??數(shù)值方法（matlab版）\NumericalMethods7.1\gauss.m

?????文件???????1584??2006-06-25?14:58??數(shù)值方法（matlab版）\NumericalMethods7.1\golden.m

?????文件???????2936??2006-06-26?06:49??數(shù)值方法（matlab版）\NumericalMethods7.1\grads.m

?????文件???????1123??2006-06-24?07:01??數(shù)值方法（matlab版）\NumericalMethods7.1\gseid.m

?????文件???????1214??2006-06-27?05:47??數(shù)值方法（matlab版）\NumericalMethods7.1\hamming.m

?????文件????????935??2006-06-30?11:35??數(shù)值方法（matlab版）\NumericalMethods7.1\heun.m

?????文件????????913??2004-01-19?05:08??數(shù)值方法（matlab版）\NumericalMethods7.1\house.m

?????文件???????1127??2006-06-28?15:48??數(shù)值方法（matlab版）\NumericalMethods7.1\invpow.m

?????文件????????907??2006-06-24?06:54??數(shù)值方法（matlab版）\NumericalMethods7.1\jacobi.m

?????文件???????1092??2006-06-28?07:24??數(shù)值方法（matlab版）\NumericalMethods7.1\jacobi1.m

?????文件????????937??2006-06-25?04:40??數(shù)值方法（matlab版）\NumericalMethods7.1\lagran.m

............此處省略34個文件信息