Вот, скажем, мы прошли второй ряд. В a[2,1] будет ноль.
for j:=1 to 3 do
a[2,j]:=a[2,j]-a[1,j]*a[2,1]/a[1,1];
Аналогично нужно пройти третий ряд, но 2 раза, поскольку необходимо "обнулить" 2 элемента.