home *** CD-ROM | disk | FTP | other *** search
- PROGRAM d12r8(input,output,dfile);
- (* driver for routine SPCTRM *)
- CONST
- m=16;
- m4=64; (* m4=4*m *)
- TYPE
- glmarray = ARRAY [1..m] OF real;
- gl4marray = ARRAY [1..m4] OF real;
- gldarray = gl4marray;
- VAR
- j,k : integer;
- ovrlap : boolean;
- p,q,w2 : glmarray;
- w1 : gl4marray;
- dfile : text;
-
- (*$I MODFILE.PAS *)
- (*$I FOUR1.PAS *)
-
- (*$I SPCTRM.PAS *)
-
- BEGIN
- glopen(dfile,'spctrl.dat');
- k := 8;
- ovrlap := true;
- spctrm(p,m,k,ovrlap,w1,w2);
- close(dfile);
- glopen(dfile,'spctrl.dat');
- k := 16;
- ovrlap := false;
- spctrm(q,m,k,ovrlap,w1,w2);
- close(dfile);
- writeln('Spectrum of data in file SPCTRL.DAT');
- writeln(' ':13,'overlapped ',' ':5,'non-overlapped');
- FOR j := 1 to m DO BEGIN
- writeln (j:6,' ':5,p[j]:13,' ':5,q[j]:13)
- END
- END.
-