d5r6.txt

Numeric Programs

implementation
//PROGRAM D5R6
//Driver for routine CHINT
uses
  unit2;
  {$R *.DFM}
Function FUNC(X:real):real;
begin
    FUNC:= Sqr(X) * (Sqr(X) - 2) * Sin(X);
end;

Function FINT(X:real):real;
var
    AAA:real;
begin
    //Integral of FUNC
    AAA:= 4 * X * (Sqr(X) - 7) * Sin(X);
    FINT:= AAA - (Sqr(Sqr(X)) - 14 * Sqr(x) + 28) * Cos(X);
end;

procedure TForm1.Button1Click(Sender: TObject);
var
  F:TextFile;
  C,CINT1:array[0..40] of real;
  X,AAA:real;  I,MVAL:integer;
const
  s1='%13.6f'; s2='##0'; NVAL = 40;
  PIO2 = 1.5707963;
  A = -PIO2;   B = PIO2;
begin
  CHEBFT(A, B, C, NVAL);
  //Test integral
  //输出计算结果到文件
  AssignFile(F, 'd:\delphi_shu\p5\d5r6.dat');
  Rewrite(F);
  Writeln(F, 'How many terms in Chebyshev evaluation?');
  //Input MVAL , between 6 and 40, MVAL=0 TO END
  MVAL:= 20;
  Writeln(F,FormatFloat(s2,MVAL));
  If (MVAL <= 0) Or (MVAL > NVAL) Then Exit;
  CHINT(A,B,C,CINT1,MVAL);
  Writeln(F, '         X         Actual    Cheby. Integ.');
  For I:= -8 To 8 do
  begin
    X:= I * PIO2 / 10;
    AAA:=CHEBEV(A,B,CINT1,MVAL,X);
    Writeln(F,Format(s1,[X]),Format(s1,[FINT(X)-FINT(-PIO2)]),
                       Format(s1,[AAA]));
  end;
  CloseFile(F);
  //屏幕显示计算结果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p5\d5r6.dat');
end;