Перевод с++ на pascal. Метод Лобачевского-Греффе

Program TheBest;
 
Function St(Value,Degree:Real):Real;
Begin
    If Degree > 0 Then
        St:=Exp(Degree*Ln(Value))
    Else
        St:=1/(Exp(Degree*Ln(Value)));
End;
 
Var
    j,s,k,i,n,step,izo:Integer;
    summ,akms,akps,b:Real;
    Aa,A,x:Array[0..1000] Of Real;
Begin
    WriteLn('Введите степень уравнения');
    ReadLn(step);
    n:=step + 1;
    WriteLn('Введите коэффициенты при переменных');
    For i:=0 To Step Do
        ReadLn(Aa[i]);
    j:=2;
    While j <= 128 Do
        Begin
            For k:=0 To Step Do
                Summ:=0;
            For s:=1 To k Do
                Begin
                    If ((k-s)<0) Or ((k-s)>step) Then
                        akms:=0
                    Else
                        akms:=Aa[k-s];
                    If((k+s)<0) Or ((k+s)>step) Then
                        akms:=0
                    Else
                        akps:=Aa[k+s];
                    summ:=summ+St(-1,s)*akms*akps;
                End;
            j:=j*2;
        End;
    A[k]:=Aa[k]*Aa[k]+2*summ;
For i:=0 To Step Do
    Aa[i]:=A[i];
b:=1.0/128.0;
for i:=0 To step Do
x[i]:=St((Aa[i+1]/Aa[i]),b);
for i:=0 To step Do
Begin
izo:=i+1;
WriteLn('X',izo,'=',x[i]);
End;
WriteLn('Подставьте корни в исходное уравнение, меняя знаки корней на противоположные, если они не обращают его в тождество');
End.