procedure EVAL ( p:ZISU; n:integer; var VALREAL,VALIMAG:ZISU );
   type ZISU= array[0..99] of real;
    var a,b            : integer;
        ARG            : ZISU;
        OMEGAc,OMEGAs  : real;
    begin
       for a:=0 to n-1 do
          ARG[a]:=(2*3.14*a)/n;
       for a:=0 to n-1 do
          begin
             VALREAL[a]:=0;
             VALIMAG[a]:=0;
             for b:=0 to n-1 do
                begin
                   OMEGAc:=cos(ARG[a]*b);
                   OMEGAs:=sin(ARG[a]*b);
                   VALREAL[a]:=VALREAL[a]+p[b]*OMEGAc;
                   VALIMAG[a]:=VALIMAG[a]+p[b]*OMEGAs
                end
          end
    end;