Вычислить произведение значений тех элементов массива, для которых справедливы указанные неравенства - Free Pascal

Формулировка задачи:

Создать массив с элементами a_kn=n*f(k)+sin(k)*g(n), где k,n=1,2,3,4; f(k)=2sin*ksin(2k-1,5)*cos(2k+1,5)-6 g(n)=|cosn²-0,51|*sin(3n-4)-4,44 Вывести его на экран в виде таблицы. Вычислить произведение значений тех элементов для которых справедливы неравенства a_kn<-1 или a_kn>1

Решение задачи: «Вычислить произведение значений тех элементов массива, для которых справедливы указанные неравенства»

Листинг программы

function F(k: byte): real;
begin
  f := 2 * sin(k) * sin(2 * k - 1.5) * cos(2 * k + 1.5) - 6
end;
 
function G(n: byte): real;
begin
  g := abs(cos(n * n) - 0.51) * sin(3 * n - 4) - 4.44
end;
 
const
  k = 4;
  n = 4;
  
var
  A: array[1..k, 1..n] of real;
  D: real;
  i, j: byte;
 
begin
  D := 1;
  for i := 1 to k do begin
    for j := 1 to n do begin
      a[i, j] := j * F(i) + sin(i) * G(j);
      write(a[i, j]:9:3);
      if (a[i, j] < -1) or (a[i, j] > 1) then D := D * a[i, j]
    end;
    writeln
  end;
  writeln('Произведение: ', D);
  readln
end.

Объяснение кода листинга программы

В этом коде реализованы две функции:

F(k: byte): real; Вычисляет значение функции для каждого k, удовлетворяющего условию.
G(n: byte): real; Вычисляет значение функции для каждого n, удовлетворяющего условию. В коде определен массив A, размером k на n, инициализированный случайными значениями. В цикле по i вычисляются значения элементов массива A по формуле a[i, j] = j F(i) + sin(i) G(j). Затем выводится значение каждого элемента массива A с шагом 3 знака после запятой. Далее проверяется условие и, если оно выполняется, то произведение D умножается на текущее значение элемента массива A. После завершения циклов выводится произведение всех элементов массива A, для которых выполняется условие. Для ввода/вывода используется функция Writeln, которая выводит строку или число, а также функцию Readln, которая считывает строку из стандартного ввода.