unit1.pas

用于开发税务票据管理的软件

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;

var
  Form1: TForm1;

implementation
//PROGRAM D10R7
//Driver for routine CORREL
uses
  unit2;
  {$R *.DFM}
procedure TForm1.Button1Click(Sender: TObject);
const
  s1='%14.6f'; s2 = '##0';
  N = 64;  N2 = 128;  PI = 3.1415926;
var
  DATA1,DATA2:array[0..64] of real;  ANS:array[0..128] of real;
  I,J:integer;   CMP:real;
  F:TextFile;
begin
  For I:=1 To N do
  begin
    DATA1[I]:=0;
    If (I > (N / 2 - N / 8)) And (I < (N / 2 + N / 8)) Then
      DATA1[I]:=1;
    DATA2[I]:=DATA1[I];
  end;
  CORREL(DATA1,DATA2, N, ANS);
  //Caculate directly
  //输出计算结果到文件
  AssignFile(F, 'd:\delphi_shu\p10\d10r7.dat');
  Rewrite(F);
  Writeln(F);
  Writeln(F, '   n       CORREL      Direct Calc');
  For I:=0 To 16 do
  begin
    CMP:=0;
    For J:=1 To N do
      CMP:=CMP + DATA1[((I + J - 1) Mod N) + 1] * DATA2[J];
    Writeln(F,'  ',FormatFloat(s2,I),Format(s1,[ANS[I+1]]),Format(s1,[CMP]));
  end;
  CloseFile(F);
  //屏幕显示计算结果
  memo1.Lines.LoadFromFile('d:\delphi_shu\p10\d10r7.dat');
end;

end.