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

[shmia]

Zitat:

Zitat von 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)

Beim Berechnen der Fakultät einer Zahl gibt es sehr schnell extrem hohe Werte.
Mit einem normalen Integer (MaxInt=2147483647)bist du da ganz schnell am Ende.
13! müsste einen 32bit Integer schon sprengen.

Der Trick an der Sache ist nun das geschickte Kürzen der Terme.
Beim Berechnen von N über K (bei dir in funcion ncr) muss alles versucht werden, damit Werte, die sich
sowieso wegkürzen, gar nicht erst berechnet werden.

[LiS]

ok... kürzen
werd ich mal versuchen
danke!

[BenBE]

Schau Dir mal in cvs.sourceforge.net/viewcvs.py/omorphia/omorphia/library/source/OMathGeneral.pas (Omorphia Projekt) die Funktion BinomialRow an.

Mit dieser kannst Du die NCR-Dinge auch sehr einfach Berechnen, ohne einen Overflow zu erhalten.

Die Extended-Implementierung hält bis Zeile 16391 mit, sollte also für deine Zwecke vollkommen ausreichen.