Метод интервалов поиск минимума ф-ции (Delphi)

[Fly@ngel03]

Во время запуска (для необходимых функци1 (-x^3)-2x+x^4 и x^x) программа полностью зависает, а если ввести функцию -х^3 пишет * Floating point overflow*
пробовала ставить тип extended ничего не меняется
ну а в остальных функций работает нормально.Какие могут быть причины такого?

procedure TForm1.Button1Click(Sender: TObject);
function f(x:real):real;
begin
f:=power(-x,3)-2*x+ power(x,4);
end;
var a,b,h,T,P,P0,P1,P2,Pmin,min,n:real;
    begin
    a:=StrToFloat(Edit1.Text);
    b:=StrToFloat(Edit2.Text);
    h:=StrToFloat(Edit3.Text);
    T:=StrToFloat(Edit5.Text);
    Pmin:=a;
    n:=0;

  Repeat
  p0:=pmin;
  P:=(f(p0+0.01)-f(p0-0.01));
  If P>=0 then h:=-abs(h) else h:=abs(h);
        repeat
        P1:=P0+h;
        P2:=P0+2*h;
        if f(P1)>=f(P0) then h:=h/2;
        if (f(P2)<f(P1)) and (f(P1)<f(P0)) then  h:=h*2;
        until (f(P0)>f(P1)) and (f(P1)<f(P2));
  Pmin:= P0+((h*(f(P1)-3*f(P0)-f(P2)))/(4*f(P1)-2*f(P0)-2*f(P2)));
  n:=n+1;
  Until abs(f(Pmin)-f(p0))<=T;
min:=(Pmin+P0)/2;
Edit4.Text:=FloatToStr(min);
Edit6.Text:=FloatToStr(n);
end;
end.