Вычислить и вывести массив значений функции

Uses Crt;
{$F+}
const n=61;
type matrix=array[1..2,1..n] of real;
type trig=function(x:real):real;
var
A:matrix;
k:integer;
function atan(x:real):real;
 begin
  atan:=sqr(arctan(x*exp(-(sqr(-x)/2))));
    end;
procedure mean(f:trig; var A:matrix);
  var
   k:integer;
    begin
     for k:=1 to n do begin
      a[1,k]:=-3+0.1*(k-1);
       a[2,k]:=f(a[1,k]);
         end;
           end;
procedure minmax(A:matrix);
  var
    min,max,c,d:real;
     k:integer;
       begin
        c:=a[2,1];
         d:=a[2,1];
           max:=a[1,k];
            min:=a[1,k];
       for k:=2 to n do begin
         if a[2,k]>c then begin
                       max:=a[1,k];
                      c:=a[2,k];
                     end;
         if a[2,k]<d then begin
                       min:=a[1,k];
                      d:=a[2,k];
                     end;
                    end;
          writeln(' ');
          writeln('Minima');
          writeln(min:0:4,' ',d:0:1);
          writeln('Maxima');
          writeln(max:0:4,' ',c:0:1);
            end;
procedure average(A:matrix);
  var
    k:integer;
     s1,s2,p,q:real;
       begin
      s1:=0;
      s2:=0;
    for k:=1 to n do begin
     s1:=s1+a[2,k];
      s2:=s2+sqr(a[2,k]);
        end;
     p:=s1/n;
     q:=sqrt(s2/n);
       writeln('Arithmetic mean',' ',p);
        writeln('Square mean',' ',q);
        end;
 
begin
Clrscr;
mean(atan,A);
for k:=1 to n do
write(a[1,k]:0:1,'  ',a[2,k]:0:3,'  ');
 minmax(A);
  average(A);
   readln;
    end.