В матрице найти произведение элементов главной диагонали минус сумма элементов побочной диагонали

program lab_5;
 
uses
  crt;
 
var
  
  a: array[1..3, 1..3] of real;
  i, j: integer;
  b: real;
 
begin
  clrscr;
  for i := 1 to 3 do
  for j := 1 to 3 do
  a[1, 1] := 0.01;
  a[1, 2] := 0.04;
  a[1, 3] := -0.03;
  a[2, 1] := 0.02;
  a[2, 2] := 0.08;
  a[2, 3] := -0.07;
  a[3, 1] := 0.1;
  a[3, 2] := -0.15;
  a[3, 3] := 0.04;
  
  for i := 1 to 3 do
    for j := 1 to 3 do
  if i = j then 
  begin;
  b:=a[i, j]*a[i,j]*a[i,j];
  end;
  
  for i := 1 to 3 do
  begin
    for j := 1 to 3 do 
      write(a[i, j], ' ');
    writeln;
  end;
  write (b);
  readkey;
end.

uses
  crt;
const
  n = 3;
  a: array[1..n, 1..n] of real = ((0.01,  0.04, -0.03),
                                  (0.02,  0.08, -0.07),
                                  ( 0.1, -0.15,  0.04));
var
  i, j: integer;
  b: real;
 
begin
  clrscr;
  for i:=1 to n do
  begin
    for j:=1 to n do
      write(a[i,j]:6:2);
    writeln;
  end;
 
  {считаем произведение элементов главной диагонали }
  b:=1;
  
  for i:=1 to n do
    b:=b*a[i,i];
    
 { a[i,n-i+1] - элемент побочной }
  for i:=1 to n do
    b:=b-a[i,n-i+1];
    
  writeln('b = ',b:0:7);
  readkey;
end.