Distribuţia Binomială: Modelare Statistică, Optimizare Numerică, cu ...

Distribuţia Binomială: Modelare Statistică, Optimizare Numerică, cu Aplicaţii în Bioinformatică şi Biochimie
end;
procedure get_nn(var nn : tip_n; n : integer; var f:text);
var
i : integer;
s : ansistring;
s1 : string[12];
vi : integer;
begin
readln(f,s);
delete(s,1,1);
for i := 0 to n-2 do begin
s1:=copy(s,1,pos(chr(9),s)-1);
delete(s,1,pos(chr(9),s));
val('$'+s1,nn[i],vi);
end;
val('$'+s,nn[n-1],vi);
end;
function bnyx(n, y, x :integer) : tip_prob;
var
ret, x_n : tip_prob;
ny , k : integer;
begin
if(x=0) then begin
if(y>0) then begin bnyx := 0.0; exit; end;
bnyx := 100.0; exit;
end;
x_n := 1.0 - x/n; ny := n - y;
if(n=x) then begin
if(ny>0) then begin bnyx := 0.0; exit; end;
bnyx := 100.0; exit;
end;
ret := 0;
if(y>0)then begin
ret := ret + ln(1.0-x_n)*y;
if(n>y)then begin
ret := ret + ln(x_n)*ny;
if(n-y>y) then y := n - y;
for k := y+1 to n do ret := ret + ln(k);
for k := 2 to n-y do ret := ret - ln(k);
end;
end else ret := ret + ln(x_n)*ny;
bnyx := 100.0*exp(ret);
end;
procedure prob_bino(n: integer; var bino : mat_prob);
var
n_2 , X , Y : integer;
begin
n_2 := n div 2;
for X := 0 to n_2+1 do //n_2+1
for Y := 0 to n do
bino[X][Y] := bnyx(n,Y,X);
for X := n downto n_2 + 2 do
for Y := 0 to n do
bino[X][Y] := bino[n-X][n-Y];
47(157)