unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs,
  StdCtrls;

type
  TForm1 = class(TForm)
    Button1: TButton;
    Memo1: TMemo;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
type
  matrx2=array of array of real;
var
  Form1: TForm1;
  X:array[0..11] of real;     QCOL:matrx2;
implementation
//PROGRAM D12R8
//Driver for routine PZEXTR
uses
  unit2;
  {$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
const
  s1='%13.6f';  s2 = '%7.4f';  NV = 4;  NUSE = 5;
var
  F:TextFile;   YEST,YZ,DY:array[0..NV] of real;
  I,J,IEST:integer;    DUM,XEST:real;
begin
  //Feed values from a rational function
  //Fn(x)=(1-x+x^3)/(x+1)^n
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p12\d12r8.dat');
  Rewrite(F);
  For I:= 1 To 10 do
  begin
    IEST:= I;
    XEST:= 1 / I;
    DUM:= 1 - XEST + Exp(3*Ln(XEST));
    For J:= 1 To NV do
    begin
      DUM:= DUM / (XEST + 1);
      YEST[J]:= DUM;
    end;
    PZEXTR(IEST, XEST, YEST, YZ, DY, NV, NUSE);
    Writeln(F);
    Writeln(F,'IEST= ',FormatFloat('##', I),' XEST= ',Format(s2,[XEST]));
    Writeln(F,'Extrap. Function: ',Format(s1,[YZ[1]]),Format(s1,[YZ[2]]),
                                  Format(s1,[YZ[3]]),Format(s1,[YZ[4]]));
    Writeln(F,'Estimated Error:  ',Format(s1,[DY[1]]),Format(s1,[DY[2]]),
                                  Format(s1,[DY[3]]),Format(s1,[DY[4]]));
  end;
  Writeln(F);
  Writeln(F,'Actual Values:      ','1.0       1.0      1.0      1.0');
  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p12\d12r8.dat');
end;

end.