implementation
//PROGRAM D1R7
//Driver for routine TOEPLZ
uses
  unit2;
  {$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
var
    X,Y:Array[0..5] of real; R:Array[0..10] of double;
    I,J:Integer;   SUM1:real;
    F:TextFile;
const
    s1='%10.4f'; s2='.###00e+00';  N = 5;  N2 = 10;
begin
  //輸入已知的方程組的系數(shù)矩陣
  For I:=1 To N2-1 do
    R[I]:= 1 / I;
  //輸入已知的方程組的右端向量B
  For I:=1 To N do
    Y[I]:= 0.1 * I;
  //輸出計(jì)算結(jié)果到文件
  AssignFile(F, 'd:\delphi_shu\p1\d1r7.dat');
  Rewrite(F);
  Writeln(F, '已知的方程組的右端向量B');
  For I:= 1 To N do
    Writeln(F,Format(s1,[Y[I]]));
  TOEPLZ(R,X,Y,N);
  Writeln(F, '計(jì)算出的方程組的解');
  For I:= 1 To N do
    Writeln(F,'     ',FormatFloat(s2,X[I]));
  Writeln(F, '將計(jì)算出的解B乘以系數(shù)矩陣，以驗(yàn)證計(jì)算結(jié)果正確');
  Writeln(F, '解乘以系數(shù)矩陣','         方程組的右端向量');
  For I:=1 To N do
  begin
      SUM1:=0;
      For J:=1 To N do
      begin
        SUM1:=SUM1 + R[N+I-J]*X[J];
      end;
      Writeln(F,Format(s1,[SUM1]),'           ',Format(s1,[Y[I]]));
  end;
  CloseFile(F);
  //屏幕顯示計(jì)算結(jié)果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p1\d1r7.dat');
end;