program matrix;

  const n = 50;

  var
    i,
    j,
    k: integer;
    s: real;
    a,
    b: array [1 .. n,1 .. n] of real;

  begin
    for i := 1 to n do
      for j := 1 to n do a[i, j] := 1.0;
    for i := 1 to n do
      for j := 1 to n do
        begin
          s := 0;
          for k := 1 to n do s := a[i, k] * a[k, j] + s;
          b[i, j] := s
        end
  end.