Re: Low hanging fruit - Accelerated random distribution functions

[Paul Kienzle]

On Feb 22, 2007, at 7:29 PM, Paul Kienzle wrote:

> On Feb 22, 2007, at 11:46 AM, David Bateman wrote:
>
> > pascal_rnd,
>
> Needs work.

Looking at the heart of pascal, it is a poisson
variate which a lambda chosen at random from a
beta distribution:

  function r = nbinrnd(n,p,varargin)
    r = randp( (1-p)./p .* randg(n,varargin{:}) );
  endfunction

I renamed it nbinrnd for compatibility.

- Paul