Конвертировать код Delphi под Pascal. Симплекс Метод

const mm = 100; nn = 100; shtraf = 1000;
 
var A : array[1..mm, 1..nn] of real;
    fun : array[1..nn] of integer; // Коэффициенты целевой функции
    m, n : integer; // m ограничений, n переменных
 
    basis : array[1..nn] of integer; // Номера базисных переменных
    i, j : integer;
    x : array[1..nn] of real; // Значения переменных при расшифровке плана
    Input,Output:text;
 
procedure solve;
var i, j, i0, j0 : integer;
    tmp : real;
    opt : boolean;
begin
    opt := false;
    repeat
        j0 := 1; i0 := 0;
        while (j0 < m+n+1) and (A[m+1, j0] >= 0) do inc(j0);
        if A[m+1, j0] >= 0 then opt := true;
 
        if not opt then begin
            tmp := 10000;
            for i := 1 to m do
                if (A[i, j0] > 0) and (A[i, m+n+1] / A[i, j0] < tmp) then
                begin
                    tmp := A[i, m+n+1] / A[i, j0];
                    i0 := i;
                end;
            // i0 - выводим, j0 - добавляем
            basis[i0] := j0; // Ввод нового элемента в базис
            // [i0, j0] - ведущий эл-т в Гауссе:
            for i := 1 to m + 1 do
                    if i <> i0 then
                    begin
                            tmp := A[i, j0];
                            for j := 1 to m + n + 1 do
                                   A[i,j] := A[i,j] - A[i0,j]*tmp/A[i0,j0];
                    end;
            tmp := A[i0, j0];
            for j := 1 to m + n + 1 do
                    A[i0, j] := A[i0, j] / tmp;
        end;
    until opt;
end;
 
begin
    Assign(Input, 'input.txt');
    Reset(Input);
    Assign(Output,'output.txt');
    Rewrite(Output);
    // Размерность матрицы
    read(n); read(m);
 
    for i := 1 to n do read(fun[i]); // Читаем коэффициенты целевой функции
 
    for i := 1 to m do
        for j := 1 to n do
            read(A[i, j]); // Читаем матрицу ограничений
 
    for i := 1 to m do
        read(A[i, n+m+1]); // Читаем правые части ограничений
 
    for i := 1 to m do // Вводим дополнительные переменные
        A[i, n+i] := 1;
    for i := 1 to m+1 do
    for j := 1 to m+1 do
        A[i,j] := 0;
   // fillchar(A[m+1], sizeof(A[m+1]), 0); // Заполняем нулями
 
    // Базис из доп. переменных
    for i := 1 to m do
        basis[i] := n + i;
    for j := 1 to n do
        A[m+1,j] := -fun[j]; // Оценки для переменных
 
    solve; // Запускаем вычисления
 
    // Вывод решения
    for i := 1 to m do
        if basis[i] <= n then
             x[basis[i]] := A[i, m+n+1];
 
    for i := 1 to n do writeLn('x[', i, '] = ', x[i]:0:3);
    writeLn('min f(x) = ', A[m+1, m+n+1]:0:3);
 
    Close(Input);
    Close(Output);
end.

for i := 1 to m+1 do
    for j := 1 to m+1 do
        A[i,j] := 0;