Вычислить функцию Р = cosA-tgB+ sin^2C

unit functin;
interface 
  type mas=array [1..100] of real;
  var n:integer;
  procedure vvod(var x,y,z:mas; var n:integer);
  procedure vivod(var x,y,z:mas; n:integer);
  function fun (var x,y,z:mas; var a,b,c,n:integer; p:real):real;
implementation 
  procedure vvod (var x,y,z:mas; var  n:integer);
var i,j,k: integer;
begin
  writeln('Размерночть массива n');
  readln(n);
  writeln('Введите элементы массива x ');
    for i:=1 to n do
      begin
      write('x[',i,']= ');
      readln(x[i]);
      end;
      
      writeln('Введите элементы массива y ');
    for j:=1 to n do
      begin
      write('y[',j,']= ');
      readln(y[j]);
      end;
      
      writeln('Введите элементы массива z ');
    for k:=1 to n do
      begin
      write('z[',k,']= ');
      readln(z[k]);
      end;
end;
 
procedure vivod (var  x,y,z: mas; n:integer);
var i,j,k: integer;
begin
  writeln('Массив x');
  for i:=1 to n do
   begin
      write(x[i]:3);   
    end;
     writeln;
    writeln('Массив y');
  for j:=1 to n do
   begin
      write(y[j]:3);    
    end;
    writeln;
    writeln('Массив z');
  for k:=1 to n do
   begin
      write(z[k]:3);
    end;
      writeln;    
end;
 
function fun (var x,y,z:mas; var a,b,c,n:integer; p:real):double;
var 
   i,j,k:integer;
begin
  for i:=1 to n do 
  if (x[i]>1) and (x[i]<10)
      then 
       inc(a);
       
       for j:=1 to n do 
  if (y[j]>1) and (y[j]<10)
      then inc(b);
      
     for k:=1 to n do
  if (z[k]>1) and (z[k]<10)
      then inc(c);
   writeln;
  P:= cos(a)-(sin(b)/cos(b))+sqr(sin(c));
  writeln('P=', p:2:3);
end;
begin
end.

uses functin;
var x,y,z:mas;
 a,b,c:integer ;
 n:integer;
 p:real;
 
begin
  vvod(x,y,z,n);
  vivod(x,y,z,n);
  fun(x,y,z,a,b,c,n,p);
end.