File: binomial.sci

package info (click to toggle)
scilab 2.6-4
  • links: PTS
  • area: non-free
  • in suites: woody
  • size: 54,632 kB
  • ctags: 40,267
  • sloc: ansic: 267,851; fortran: 166,549; sh: 10,005; makefile: 4,119; tcl: 1,070; cpp: 233; csh: 143; asm: 135; perl: 130; java: 39
file content (13 lines) | stat: -rw-r--r-- 399 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
 
function pr=binomial(p,n)
//binomial probability distribution
//pr(k) = probability (X=k) , with X=B(n,p)
//Example:
//n=10;p=0.3;plot2d3("onn",(0:n)',binom(p,n)');
//n=100;p=0.3;
//moy=n*p;sigma=sqrt(n*p*(1-p));
//x=( (0:n)-moy )/sigma;
//plot2d(x,sigma*binom(p,n));
//deff('y=gauss(x)','y=1/sqrt(2*%pi)*exp(-(x.^2)/2)')
//plot2d(x,gauss(x));
z=poly(0,'z');
pr=coeff((1-p+z)^n).*horner(z^(0:n),p);