producto de matrices pascal

procedure ProductoAxB_C(a : MatrizMxN; b : MatrizNxP; var c : MatrizMxP);
var
    i, j, k, aux : Integer;
begin
    aux := 0;
    for i := 1 to M do
        for j := 1 to P do
        begin
            for k := 1 to N do
            begin
                c[i,j] := a[i,k] * b[k,j];
                aux := aux + c[i,j];
            end;
            c[i,j] := aux;
            aux := 0;
        end;
end;