1、说明:
本书中所有的常用数值算法子过程按书中的章数分别放在以C开头的子目录中。
所有这些为验证上述子过程而编的验证过程按书中的章数分别放在以D开头的子目录中。
所有为验
源代码在线查看: unit2.pas
unit Unit2;
interface
Procedure FGAUSS(X:real; A:array of real;
var Y:real;var DYDA:array of real; NA:integer);
implementation
Procedure FGAUSS(X:real; A:array of real;
var Y:real;var DYDA:array of real; NA:integer);
var
I,II:integer; ARG,EX,FAC:real;
begin
Y:=0;
For II:=1 To (NA DIV 3) do
begin
I:=3*II-2;
ARG:=(X - A[I + 1]) / A[I + 2];
EX:=Exp(-Sqr(ARG));
FAC:=A[I] * EX * 2 * ARG;
Y:=Y + A[I] * EX;
DYDA[I]:=EX;
DYDA[I + 1]:=FAC / A[I + 2];
DYDA[I + 2]:=FAC * ARG / A[I + 2];
end;
end;
end.