-
Procedure MNBRAK(var AX, BX, CX, FA, FB, FC:real);
const
GOLD = 1.618034; GLIMIT = 100; TINY = 1E-20;
var
DUM,R,Q,U,ULIM,FU:real; DONE:boolean;
begin
FA:=FUNC(AX);
FB:=FUNC
http://www.codebf.com/read/168081/5448973
-
Procedure SIMP1(var A:matrx2; MP, NP, MM:integer; LL:array of integer;
NLL, IABF:integer;var KP:integer; var BMAX:real);
var
K:integer; TEST:real;
begin
KP:=LL
http://www.codebf.com/read/168081/5448975
-
procedure QSIMP(A, B:real; var S:real);
const
EPS = 0.000001; JMAX = 20;
var
J:integer; OST,OS,ST:real;
begin
OST:=-1E+30;
OS:=-1E+30;
For J:=1 To JMAX do
begin
http://www.codebf.com/read/168081/5448976
-
procedure QGAUS(A, B:real;var SS:real);
var
X, W:array[0..5] of real;
XM,XR,DX:real; J:integer;
begin
X[1]:=0.1488743389; X[2]:=0.4333953941;
X[3]:=0.6794095682; X[4]:=0.8
http://www.codebf.com/read/168081/5448977
-
procedure QROMB(A, B:real; var SS:real);
CONST
EPS = 0.000001; JMAX = 20;
VAR
S,H:ARRAY[0..21] of real;
JMAXP,J,K,KM:integer; DSS:real;
begin
JMAXP:=JMAX + 1;
K:=5;
http://www.codebf.com/read/168081/5448978
-
implementation
procedure TRAPZD(A, B:real;var S:real; N:integer);
var
J,IT:integer;
TNM,DEL,X,SUM:real;
begin
If N = 1 Then
begin
S:=0.5 * (B - A) * (FUNC(A) + FUNC(B))
http://www.codebf.com/read/168081/5448979
-
procedure MIDPNT(A, B:real;var S:real; N:integer);
var
TNM,DEL,DDEL,X,SUM:real;
J,IT:integer;
begin
If N = 1 Then
begin
S:=(B - A) * FUNC(0.5 * (A + B));
IT:=1;
http://www.codebf.com/read/168081/5448980
-
Function SQU(X, BB:real):real;
begin
SQU:=2 * X * FUNC(BB - Sqr(X));
end;
procedure MIDSQU(AA, BB:real; var S:real; N:integer);
var
A,B,TNM,X,DEL,DDEL,SUM:real; IT,J:integer;
begin
http://www.codebf.com/read/168081/5448981
-
Function INF(X:real):real;
begin
INF:=FUNC(1 / X) / Sqr(X);
End;
procedure MIDINF(AA, BB:real; var S:real; N:integer);
var
A,B,TNM,X,DEL,DDEL,SUM:real; J:integer;
begin
B:=1 / AA;
http://www.codebf.com/read/168081/5448982
-
Function SQL(X, AA:real):real;
begin
SQL:=2 * X * FUNC(AA + Sqr(X));
end;
procedure MIDSQL(AA, BB:real; var S:real; N:integer);
var
A,B,TNM,X,DEL,DDEL,SUM:real; IT,J:integer;
begin
http://www.codebf.com/read/168081/5448983