unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs,
  StdCtrls;
procedure ASUB(XIN:array of real;var XOUT:array of real);
procedure ATSUB(XIN:array of real;var XOUT:array of real);

type
  TForm1 = class(TForm)
    Button1: TButton;
    Memo1: TMemo;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;

implementation
//PROGRAM D1R9
//Driver for routine SPARSE
uses
  unit2;
  {$R *.DFM}
procedure ASUB(XIN:array of real;var XOUT:array of real);
var
   I,N:integer;
begin
    N:= 20;
    XOUT[1]:=XIN[1] + 2 * XIN[2];
    XOUT[N]:=-2 * XIN[N - 1] + XIN[N];
    For I:=2 To N - 1 do
        XOUT[I]:=-2 * XIN[I - 1] + XIN[I] + 2 * XIN[I + 1];
end;
procedure ATSUB(XIN:array of real;var XOUT:array of real);
var
   I,N:integer;
begin
    N:=20;
    XOUT[1]:=XIN[1] - 2 * XIN[2];
    XOUT[N]:=2 * XIN[N - 1] + XIN[N];
    For I:=2 To N - 1 do
        XOUT[I]:=2 * XIN[I - 1] + XIN[I] - 2 * XIN[I + 1];
end;

procedure TForm1.Button1Click(Sender: TObject);
var
    X,B,BCMP:Array[0..20] of real;
    I,J:Integer;   RSQ:real;
    F:TextFile;
const
    s1='%10.1f'; s2='%10.6f';  N = 20;
begin
  //輸入已知的方程組的系數矩陣
  For I:=1 To N do
  begin
      X[I]:= 0;
      B[I]:= 1;
  end;
  B[1]:= 3;
  B[N]:= -1;
  SPARSE(B,N,X,RSQ);
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p1\d1r9.dat');
  Rewrite(F);
  Writeln(F,'Sum-squared residual:  ',Format(s2,[RSQ]));
  Writeln(F, '計算出的方程組的解');
  For I:= 1 To 4 do
      Writeln(F,Format(s2,[X[1+(I-1)*5]]),Format(s2,[X[2+(I-1)*5]]),
                 Format(s2,[X[3+(I-1)*5]]),Format(s2,[X[4+(I-1)*5]]),
                 Format(s2,[X[5+(I-1)*5]]));
  Writeln(F, '將計算出的解B乘以系數矩陣，以驗證計算結果正確');
  ASUB(X,BCMP);
  Writeln(F, '解乘以系數矩陣','         方程組的右端向量');
  For I:=1 To N do
     Writeln(F,Format(s1,[BCMP[I]]),'           ',Format(s1,[B[I]]));
  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p1\d1r9.dat');
end;
end.