d1r6.txt

最近开发的,比较使用,曲线拟合,谢谢下载 ,没问题

implementation
//PROGRAM D1R6
//Driver for routine VANDER
uses
  unit2;
  {$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
var
    X,Q,W,TERM:Array[0..5] of real;
    I,J:Integer;   SUM1:real;
    F:TextFile;
const
    s1='%10.1f'; s2='%13.4f';  N = 5;
begin
  //输入已知的方程组的系数矩阵
  X[1]:=1;         X[2]:=1.5;
  X[3]:=2;         X[4]:=2.5;
  X[5]:=3;
  //输入已知的方程组的右端向量B
  Q[1]:=1;         Q[2]:=1.5;
  Q[3]:=2;         Q[4]:=2.5;
  Q[5]:=3;
  //输出计算结果到文件
  AssignFile(F, 'd:\delphi_shu\p1\d1r6.dat');
  Rewrite(F);
  Writeln(F, '已知的方程组的右端向量B');
  For I:= 1 To N do
    Writeln(F,Format(s1,[Q[I]]));
  VANDER(X,W,Q,N);
  Writeln(F, '计算出的方程组的解');
  For I:= 1 To N do
    Writeln(F,Format(s2,[W[I]]));
  Writeln(F, '将计算出的解B乘以系数矩阵，以验证计算结果正确');
  Writeln(F, '解乘以系数矩阵','         方程组的右端向量');
  SUM1:=0;
  For I:=1 To N do
  begin
      TERM[I]:=W[I];
      SUM1:=SUM1 + W[I];
  end;
  Writeln(F,Format(s1,[SUM1]),'           ',Format(s1,[Q[1]]));
  For I:=2 To N do
  begin
      SUM1:=0;
      For J:=1 To N do
      begin
        TERM[J]:=TERM[J]* x[J];
        SUM1:=SUM1 + TERM[J];
      end;
      Writeln(F,Format(s1,[SUM1]),'           ',Format(s1,[Q[I]]));
  end;
  CloseFile(F);
  //屏幕显示计算结果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p1\d1r6.dat');
end;