Function CHEBEV(A, B:real; C:array of real; M:integer; X:real):real;
var
D,DD,Y,Y2,SV:real; J:integer;
begin
If (X - A) * (X - B) > 0 Then ShowMessage('X not in range.');
D:=0;
procedure CHEBPC(C:array of real;var D:array of real; N:integer);
var
DD:array[0..50] of real; J,K:integer;
SV:real;
begin
For J:=1 To N do
begin
D[J]:=0 ;
DD[
procedure PCSHFT(A, B:real;var D:array of real; N:integer);
var
FAC,CONST1:real; J,K:integer;
begin
CONST1:=2/ (B - A);
FAC:=CONST1;
For J:=2 To N do
begin
D[J]:=D
procedure CHDER(A,B:real;C:array of real;var CDER:array of real;N:integer);
var
J:integer; CON:real;
begin
CDER[N]:=0;
CDER[N - 1]:=2 * (N - 1) * C[N];
If N >= 3 Then
be
procedure CHEBFT(A, B:real; var C:array of real; N:integer);
var
K,J:integer; SUM,Y,BMA,BPA,FAC:real;
F:array[0..50] of real;
const
NMAX = 50; PI = 3.14159265358979;
begin
procedure EULSUM(var SUM:real;TERM:real;JTERM:integer);
var
J:Integer; TMP,DUM:real;
WKSP:array[0..41] of real;
begin
If JTERM = 1 Then
begin
NTERM^:=1;
WKSP[1
