Triangularizar matriz

Programa que se utiliza como función en otros, como es el caso de resolución de sistemas mediante eliminación gaussiana, ortogonalización... Pienso que puede ser útil por si alguien lo quisiera utilizar como función en otro programa principal.

program triangulariza(input,output);
Const
  tanMax=20;
  precision=0.000001;

Type
  tpVector=array[1..tanMax] of real;
  tpMatriz=array[1..tanMax] of tpVector;
  tpDimension=1..tanMax;
  tpContador=1..MAXINT;

Var
  matriz:tpMatriz;
  min,pasoF,pasoC,filas,columnas:tpDimension;
  cont:tpContador;


Procedure LeerDatos(VAR f,c:tpDimension; VAR m:tpMatriz);
  Var
    i,j:tpContador;

  begin
    writeln('Introduzca el numero de filas: '); readln(f); writeln;
    writeln('Introduzca el numero de columnas: '); readln(c); writeln;
    for i:=1 to f do begin
      writeln('Termino de la fila ',i,': ');
      for j:=1 to c do begin
        write('Termino (',i,',',j,'): ');readln(m[i,j])
      end;
      writeln
    end;
    writeln
  end;


Procedure SeleccionarFila(f,c,pF,pC:tpDimension; VAR m:tpMatriz);
  Var
    i,j:tpContador;
    aux:real;

  Begin
    i:=pF;
    while (m[i,pC]<precision) and (i<f) and (m[i,pC]>-precision) do
      i:=i+1;
    If i<>f then
      begin
        for j:=pC to c do begin
          aux:=m[pF,j];
          m[pF,j]:=m[i,j];
          m[i,j]:=aux
        end
      end
  end;


Procedure restarFilas(f,c:tpDimension; VAR pF:tpDimension; VAR pC:tpDimension; VAR m:tpMatriz);
  var
    i,j:tpContador;
    aux:real;

  begin
    if m[pF,pC]<>0 then begin
      for i:=(pF+1) to f do begin
        aux:=m[i,Pc];
        for j:=pC to c do
          m[i,j]:=m[i,j]-aux*(m[pF,j]/m[pF,pC]);
      end;
      pF:=pF+1
    end;
    pC:=pC+1
  end;


procedure MinDim(f,c:tpDimension; VAR mini:tpDimension);
  begin
    if c<f then mini:=c else mini:=f
  end;


procedure EscribirMatriz(m:tpMatriz);
  var
    i,j:tpContador;

  begin
    for i:=1 to filas do begin
      for j:=1 to columnas do
        write(m[i,j]:7:2);
      writeln
    end
  end;

begin
  LeerDatos(filas,columnas,matriz);
  MinDim(filas,columnas,min);
  EscribirMatriz(matriz);writeln; writeln;
  pasoF:=1; pasoC:=1;
  for cont:=1 to min do
    begin
      SeleccionarFila(filas,columnas,pasoF,pasoC,matriz);
      restarFilas(filas,columnas,pasoF,pasoC,matriz)
    end;
  EscribirMatriz(matriz);
  readln
end.


Un saludo y que les sirva.