Перевести на Паскаль

[Camelot_2012]

Всем привет! Помогите перевести код в Паскаль! Заранее благодарен!!!

double ab(double x1,double x2,double z2,double &a,double &b,double &N0)
{  int i=0;
   double y1,y2,d_alfa=1;
   y1=f(x1,x2);
   N0++;
   do{  
y2=f(x1,x2-d_alfa*(i+1)*z2);
N0++;
i++;
y1=y2;
      }while(y1<y2);
     b=(i+1)*d_alfa;
     if(i==0) a=0;
     else a=(i-1)*d_alfa;
     return a,b;
}
double minimiz_f_ZS(double x1,double x2,double z1,double z2,double a,double b,double &N,double &alfa)
{
  double y1,y2,t,alfa1,alfa2;
  t=(1+sqrt(5))/2;
  alfa1=b-(b-a)/t;
  alfa2=a+(b-a)/t;
  y1=f(x1-alfa1*z1,x2-alfa1*z2);
  y2=f(x1-alfa2*z1,x2-alfa2*z2);
  do{
     if(y1<y2){b=alfa2;alfa2=alfa1;y2=y1;
       alfa1=a+b-alfa2;
       y1=f(x1-alfa1*z1,x2-alfa1*z2);
               N++;
       }
    else{a=alfa1;alfa1=alfa2;y1=y2;
  alfa2=a+b-alfa1;
  y2=f(x1-alfa2*z1,x2-alfa2*z2);
  N++;
  }
     }while(b-a>2*E);
     if(y1<y2) b=alfa2;
     else b=alfa1;
     alfa=(a+b)/2;
     return alfa;
}
void main()
{
 double N1,N0,N,x1[100],x2[100],z1,z2,alfa,a,b,y1,y2;
 int i,k;
 k=0;
 N=0;
 x1[0]=-1;
 x2[0]=0;
 z1=dx1(x1[k],x2[k]);
 z2=dx2(x1[k],x2[k]);
 N1=2;
 clrscr();
 cout<<setw(5)<<"x1"<<setw(15)<<"x2"<<endl;
 do{
     ab(x1[k],x2[k],z2,a,b,N0);
     minimiz_f_ZS(x1[k],x2[k],z1,z2,a,b,N,alfa); //Z.S.
     x1[k+1]=x1[k]-alfa*z1;
     x2[k+1]=x2[k]-alfa*z2;
     z1=dx1(x1[k+1],x2[k+1]);
     z2=dx2(x1[k+1],x2[k+1]);
     N1=N1+2;
     k++;
     cout<<x1[k]<<setw(15)<<x2[k]<<endl;
   }while(sqrt(z1*z1+z2*z2)>E);
   cout<<endl;
   cout<<"x1~="<<x1[k]<<endl<<"x2~="<<x2[k]<<endl
   <<"y~="<<f(x1[k],x2[k]);
   cout<<endl<<"N="<<N+N0+N1<<endl
   <<"k="<<k;
   getch();
}

[grandmixer]

В коде не было переменных dx1,dx2,E и еще какой-то функции f()
что смог перевел xD

f

unction ab(var x1,x2,z2,a,b,N0:double):double;
var i:integer;
y1,y2,d_alfa:double;
begin
i:=1;
y1:=1;
y2:=1;
d_alfa:=1;
  y1:=f(x1,x2);
  N0:=N0+1;

  while y1<y2 do
  begin
        y2:=f(x1,x2-d_alfa*(i+1)*z2);
        N0:=N0+1;
        inc(i);
        y1:=y2;
  end;

  b:=(i+1)*d_alfa;

  if i=0 then a:=0
  else a:=(i-1)*d_alfa;
end;

function minimiz_f_ZS(var x1,x2,z1,z2,a,b,N,alfa:double):double;
var
y1,y2,t,alfa1,alfa2:double;
begin
  t:=(1+sqrt(5))/2;
  alfa1:=b-(b-a)/t;
  alfa2:=a+(b-a)/t;
  y1:=f(x1-alfa1*z1,x2-alfa1*z2);
  y2:=f(x1-alfa2*z1,x2-alfa2*z2);

  while b-a>2*E do
  begin
        if y1<y2 then
        begin
              b:=alfa2;
              alfa1:=a+b-alfa2;
              y1:=f(x1-alfa1*z1,x2-alfa1*z2);
              N:=N+1;
        end
        else begin
              a:=alfa1;
              alfa2:=a+b-alfa1;
              y2:=f(x1-alfa2*z1,x2-alfa2*z2);
              N:=N+1;
        end;
  end;
  if y1<y2 then b:=alfa2
  else b:=alfa1;

  alfa:=(a+b)/2;
  minimiz_f_ZS:=alfa;
end;

var
  N1,N0,N,z1,z2,alfa,a,b,y1,y2:double;
  x1,x2:array [0..100] of double;
  i,k:integer;
begin
k:=0;
N:=0;
x1[0]:=-1;
x2[0]:=0;
z1:=dx1(x1[k],x2[k]);
z2:=dx2(x1[k],x2[k]);
N1:=2;
writeln('     ','x1','               ','x2');

while sqrt(z1*z1+z2*z2)>E do
begin
     ab(x1[k],x2[k],z2,a,b,N0);
     minimiz_f_ZS(x1[k],x2[k],z1,z2,a,b,N,alfa);
     x1[k+1]:=x1[k]-alfa*z1;
     x2[k+1]:=x2[k]-alfa*z2;
     z1:=dx1(x1[k+1],x2[k+1]);
     z2:=dx2(x1[k+1],x2[k+1]);
     N1:=N1+2;
     inc(k);
     writeln(x1[k],'               ',x2[k]);
end;

writeln('x1~=',x1[k]);
writeln('x2~=',x2[k]);
writeln('y~=',f(x1[k],x2[k]));
writeln('N=',N+N0+n1);
writeln('k=',k);

end.