For small means, the number of terms must be greater.  Taking no chances,
and using 50 times the mean,



limit=: 3 : 0
NB. threshold for mean y
(+/\(^-y)**/\1,y%}.i.50*y)I.1.0
)

limit 2

22

and the cumulative distribution function reaches 1 only when x=22.

plot(;limit"0)i.100

Adjusting the comparison tolerance may be a solution.

On Sat, Aug 9, 2008 at 4:44 PM, Raul Miller <[EMAIL PROTECTED]> wrote:

> On Sat, Aug 9, 2008 at 4:37 PM, Robert Cyr <[EMAIL PROTECTED]> wrote:
> > Sory about that
> >
> > try instead:
> >
> > plot(;limit"0)i.100
>
> That does succeed, mechanically.
>
> However:
>   limit 2
> 10
>   >./poissonrand 2 1e6
> 13
>
> --
> Raul
> ----------------------------------------------------------------------
> For information about J forums see http://www.jsoftware.com/forums.htm
>
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm

Reply via email to