из Pascal в С++

program pminf;
{квадратичной интерполяции-экстраполяции}
uses crt;
var xm, a, b, e, h, x: real; v: integer;

function fun(x: real): real;
{ вычисление значения функции в точке }
begin
  fun:=0.1*x*x*x-2*x*x+10*sqrt(x)
end;

function mkv(h, x1, e: real; var xm: real): real;
{ функция реализует метод квадратичной интерполяции-экстраполяции }
var f0, f1, f2, x0, x2, c, b: real;
begin
  xm:=x1;
  repeat
      x1:=xm;     x0:=x1-h; x2:=x1+h;
      f0:=fun(x0); f1:=fun(x1); f2:=fun(x2);    c:=1/2/h/h*(f0-2*f1+f2);
      b:=((-1)*f0*(2*x1+h)+4*f1*x1-f2*(2*x1-h))/2/h/h;
      xm:=(-1)*b/2/c
  until ((xm-x1)<e);
  mkv:=fun(xm);
end;
[COLOR="Red"]//++++++++++++++все что дальше не обязательно++++++++++++++++++[/COLOR]
begin  clrscr;
  writeln ('Нахождение минимума функции.');
  writeln ('Выберите метод вычисления минимума:');
  writeln ('1. Метод золотого сечения.');
  writeln ('2. Метод квадратичной интерполяции-экстраполяции.');
  writeln ('- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -');
  write ('Введите номер варианта: '); read (v);
  case v of
    1:
      begin
        write ('Координата начала отрезка : '); read (a);
        write ('Координата конца отрезка  : '); read (b);
        write ('Точность                  : '); read (e);
        writeln ('- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -');
        h:=mzols(a,b,e,xm);
        writeln ('F(x)=',h:12:5, ' x=',xm:12:5);
      end;
    2:
      begin
        write ('Начальный шаг поиска : '); read (h);
        write ('Точность             : '); read (e);
        write ('Начальное значение x : '); read (x);
        writeln ('- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -');
        h:=mkv(h,x,e,xm);
        write ('F(x)= ', h:12:5, ' x=',xm:12:5)
      end
  else { иначе case}
    writeln ('Неверный ввод!')
  end; {конец case}
end.