bernoulli programm

[LiS]

hi
ich hab mir mal auf die schnelle ein bernoulli programm geschrieben, allerdings gibt es bei geringen zahlen schon falsche ergebnisse. (bis 20 versuche funktioniert es noch)
hier der quelltext

Delphi-Quellcode:

function fak(k:integer):integer;
var i:integer;
begin
  result := 1;
  for i := 2 to k do result := result * i;
end;

procedure TForm1.Button1Click(Sender: TObject);
var
  p:extended;
function ncr(n,k:integer):integer;
var i:integer;
begin
  result := n-k+1;
  for i := 2 to k do result := result*(n-k+i);
  result := round(result / fak(k));
end;
begin
  p:=0;
  p:=power(strtofloat(edit1.Text),strtoint(edit3.Text));
  p:=p*power(1-strtofloat(edit1.Text),(strtoint(edit2.Text)-strtoint(edit3.Text)));
  p:=p*ncr(strtoint(edit2.Text),strtoint(edit3.Text));
  p:=100*p;
  label4.Caption := floattostr(p);
end;

edit1 = Wahrscheinlichkeit
edit2 = Versuche
edit3 = Erfolge