Биноминальное распределение. Delphi

[boss89]

Нужна помощь в реализации формулы биноминального распределения:

Биномиальное распределение, распределение вероятностей числа появлений некоторого события при повторных независимых испытаниях. Если при каждом испытании вероятность появления события равна р, причём 0 £ p £ 1, то число m появлений этого события при n независимых испытаниях есть случайная величина, принимающая значения m = 1, 2,.., n с вероятностями

где q = 1 — p, a — биномиальные коэффициенты (отсюда название Б. р.). Приведённая формула иногда называется формулой Бернулли. Математическое ожидание и дисперсия величины m, имеющей Б. р., равны М (m) = np и D (m) = npq, соответственно. При больших n, в силу Лапласа теоремы, Б. р. близко к нормальному распределению, чем и пользуются на практике. При небольших n приходится пользоваться таблицами Б. р.

Кто знает как реализовать помогите пожалуйста..
Заранее спасибо!

[boss89]

Народ помогите пожалуйста! Очень надо..

[p51x]

Что вам конкретно надо?

[boss89]

Cобственно нужно вот эту формулу в делфи написать чтобы затем включить её в цикл и находить чему равно Pn..

[p51x]

Что сложного? Тип Integer взят для примера (возможно потребуется длинная арифметика).

function Factorial(n : Integer):Integer;
function P(n, m, p, q:Integer):Integer;
begin
  Result := Factorial(n)/(m*Factorial(n-m)) * Power(p, m) * Power(q, n-m);
end;

[boss89]

чёт не запускается..
А как кстати сделать следующее:
нужно чтобы он находил Pn и накапливал его в Е, допустим что Е:=первоначально равно 0, а потом Pn складывается с Е:


function Factorial(n : Integer):Integer;
function P(n, m, p, q:Integer):Integer;
begin
Pn := Factorial(n)/(m*Factorial(n-m)) * Power(p, m) * Power(q, n-m);
E:=pn+E;
end;
НУ и потом Е должна вставлятся в формулу ниже, до тех пор пока условия задачи не выполнятся.

function Zatrati ( t : Integer ) : Real;
begin
if t <> 0 then Result := N * ( C1*Pt*E) + C2 ) / T
else Result := 0;
end;

По такой схеме Результат будет расчитываться?

[p51x]

чёт не запускается..

А я и не говорил, что это точное решение... Я даже про тип сказал.

function Factorial(n : Integer):Integer;

Ее то написать надо.

нужно чтобы он находил Pn и накапливал его в Е, допустим что Е:=первоначально равно 0, а потом Pn складывается с Е:

Где-нибудь после описания функций.

e := 0;
for(...;...;...) e = e + P(...);