implementation
//PROGRAM D10R5
//Driver for routine COSFT
uses
  unit2;
  {$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
const
  s1='%13.6f'; s2 = '##';
  NP = 16;  WIDTH1 = 30;  PI = 3.1415926;  EPS = 0.001;
var
  DATA:array[0..16] of real;
  I,J,NLIM:integer;   PER,BIG,SCAL1,SMALL:real;
  STR1:String;  F:TextFile;
begin
  //輸出計算結果到文件
  AssignFile(F, 'd:\delphi_shu\p10\d10r5.dat');
  Rewrite(F);
  Writeln(F, 'Period of sinusoid in channels (2-',FormatFloat(S2,NP),
                                               ' ; or 0 to stop)');
  PER:= 16;
  Writeln(F,FormatFloat(S2,PER));
  If PER <= 0 Then Exit;
  For I:= 1 To NP do
    DATA[I]:= COS(2 * PI * (I - 1) / PER);
  COSFT(DATA, NP,+1);
  BIG:= -1.e10;
  SMALL:= 1e10;
  For I:= 1 To NP do
  begin
    If DATA[I] < SMALL Then SMALL:= DATA[I];
    If DATA[I] > BIG Then BIG:= DATA[I];
  end;
  SCAL1:= Width1 / (BIG - SMALL);
  For I:= 1 To NP do
  begin
    NLIM:= Trunc(SCAL1*(DATA[I]- SMALL) + EPS);
    STR1:='';
    For J:= 1 To NLIM + 1 do
      STR1:=STR1 + '*';
    Writeln(F,'  ',FormatFloat(s2,I),'  ', STR1);
  end;
  COSFT(DATA, NP,-1);
  BIG:= -1e10;
  SMALL:= 1e10;
  For I:= 1 To NP do
  begin
    If DATA[I] < SMALL Then SMALL:= DATA[I];
    If DATA[I] > BIG Then BIG:= DATA[I];
  end;
  SCAL1:= Width1 / (BIG - SMALL);
  For I:= 1 To NP do
  begin
    NLIM:= Trunc(SCAL1 * (DATA[I] - SMALL) + EPS);
    STR1:='';
    For J:= 1 To NLIM + 1 do
      STR1:=STR1 + '*';
    Writeln(F,'  ',FormatFloat(s2,I),'  ', STR1);
  end;
  CloseFile(F);
  //屏幕顯示計算結果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p10\d10r5.dat');
end;