<?php
function nCr($n, $k){
   if ($k > $n)
     return NaN;
   if (($n - $k) < $k)
     return nCr($n, ($n - $k));
   $return = 1;
   for ($i=0; $i<$k; $i++){
     $return *= ($n - $i) / ($i + 1);
   }
   return $return;
}

function binompdf($n, $p, $k){
    $return = nCr($n, $k) * pow($p, $k) * pow((1 - $p), ($n - $k));
    return $return;
}
?>