Вычислить произведение значений тех элементов массива, для которых справедливы неравенства - C (СИ)

var
  a: array[1..4, 1..4] of real;
  k, n: byte;
  f, g, p1, p2: real;
 
begin
  for k := 1 to 4 do
  begin
    for n := 1 to 4 do
    begin
      p1 := 1;
      p2 := 1;
      f := 2 * sin(k) * sin(2 * k - 1.5) * cos(2 * k + 1.5) - 6;
      g := abs(cos(sqr(n)) - 0.51) * sin(3 * n - 4) - 4.44;
      a[k, n] := n * f + sin(k) * g;
      write(a[k, n]:8:2);
    end;
    writeln;
  end;
  for k := 1 to 4 do
    for n := 1 to 4 do
      if a[k, n] < -1.0 then
        p1 := p1 * a[k, n];
  writeln('произведение элементов <-1: ', p1:8:2);
  if a[k, n] > 1.0 then begin
    p2 := p2 * a[k, n];
    writeln('произведение элементов >1: ', p2:8:2);
  end else
    writeln('нет элементов >1');
end.