Переписать программу с Pascal на С - C (СИ)

#include <stdio.h>
#include <math.h>
 
/*
 * Generated by TPTC - Translate Pascal to C
 *     Version 1.7 03/26/88   (C) 1988 S.H.Smith
 */
 
  #define eps  0.00001    /*  Точность определения минимума  */ 
  #define alfa 1.0    /*  Коэффициент отражения          */ 
  #define beta 0.5    /*  Коэффициент сжатия             */ 
  #define gamma 2.0    /*  Коэффициент растяжения         */ 
 
  #define t    0.2  
 
  #define n_max 6    /*  Макс. количество переменных  */ 
 
 
  double       x[n_max+1], y[n_max+1]; 
 
  double       fh, fl, f4, f5, f6; 
  int      h, l, it; 
 
 
  /*  *** Оптимизируемая функция ***  */ 
 
double       f(int      k)
{ double       r1, r2; 
 
  r1 = y[k] - x[k] * x[k]; r2 = 1 - x[k]; 
  return 100 * r1 * r1 + r2 * r2; 
} 
 
 
 
double       maxf(void)
{ double       f1, f2, r; 
 
  f1 = f(1); f2 = f(2); r = f(0); 
  h = 0; 
  if (r < f1) { 
    r = f1; h = 1; 
  } 
  if (r < f2) { 
    r = f2; h = 2; 
  } 
  return r; 
} 
 
 
double       minf(void)
{ double       f1, f2, r; 
 
  f1 = f(1); f2 = f(2); r = f(0); 
  l = 0; 
  if (f1 < r) { 
    r = f1; l = 1; 
  } 
  if (f2 < r) { 
    r = f2; l = 2; 
  } 
  return r; 
} 
 
 
 
 
 
void main(int   argc,
     char  *argv[])
{ 
  int      i, flag; 
  double       r, r1, x0, y0, x1, y1; 
    x0 = -1.2; y0 = 1; 
    x[0] = x0 - 0.5 * t; y[0] = y0 - t * sqrt(3) / 6; 
    x[1] = x0; y[1] = y0; r = f(1); y[1] = y0 + t * sqrt(3) / 3; 
    x[2] = x0 + 0.5 * t; y[2] = y[0]; 
 
    it = 0; 
    printf(" it=%5d x=%8.4f y=%8.4f f=%8.4f\n",it,x0,y0,r); 
    do { 
        fh = maxf(); fl = minf(); 
        x[3] = 0.5 * (x[0] + x[1] + x[2] - x[h]); 
        y[3] = 0.5 * (y[0] + y[1] + y[2] - y[h]); 
 
        x[4] = (1 + alfa) * x[3] - alfa * x[h]; 
        y[4] = (1 + alfa) * y[3] - alfa * y[h]; 
 
        f4 = f(4); 
        if ((f4 < fl)) { 
            x[5] = (1 - gamma) * x[3] + gamma * x[4]; 
            y[5] = (1 - gamma) * y[3] + gamma * y[4]; 
            f5 = f(5); 
            if ((f5 < fl)) { x[h] = x[5]; y[h] = y[5]; 
            } 
            else { x[h] = x[4]; y[h] = y[4]; 
            } 
            goto TheEnd;
        } 
 
        flag = 0; 
        for (i = 0; i <= 2; i++) { 
          if (((i != h) && (f4 > f(i)))) flag++; 
        } 
        if ((flag == 2)) { 
            x[6] = beta * x[h] + (1 - beta) * x[3]; y[6] = beta * y[h] + (1 - beta) * y[3]; 
            x[h] = x[6]; y[h] = y[6]; 
 
            goto TheEnd;
        } 
 
        if ((f(4) < fh)) { 
            for (i = 0; i <= 2; i++) { 
              x[i] = 0.5 * (x[i] + x[l]); y[i] = 0.5 * (y[i] + y[l]); 
            } 
        } 
        else { 
            x[h] = x[4]; y[h] = y[4]; 
        } 
 
        
TheEnd: r = 0; 
        for (i = 0; i <= 2; i++) { 
          r1 = f(i) - f(3); r = r + r1 * r1; 
        } 
        r = sqrt(r / 3); 
        it++; 
 
        r1 = (f(0) + f(1) + f(2)) / 3; 
        x0 = (x[0] + x[1] + x[2]) / 3; 
        y0 = (y[0] + y[1] + y[2]) / 3; 
 
        if (it % 20 == 0) 
          printf(" it=%5d x=%8.4f y=%8.4f f=%8.4f\n",it,x0,y0,r); 
    }  while (!((r < eps))); 
 
    printf(" it=%5d x=%8.4f y=%8.4f f=%8.4f\n",it,x0,y0,r); 
}