Rev 20 | Rev 3265 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed
| Rev | Author | Line No. | Line |
|---|---|---|---|
| 20 | reyssat | 1 | !if $wims_read_parm!=slib_header |
| 2 | !goto proc |
||
| 3 | !endif |
||
| 4 | slib_title=cumulative distribution function of Poisson law |
||
| 5 | slib_parms=2\ |
||
| 2073 | bpr | 6 | ,[data]\ |
| 7 | 1,p |
||
| 20 | reyssat | 8 | slib_author=Bernadette PERRIN-RIOU |
| 9 | slib_out=values on data of the cumulative distribution function of the Poisson law of parameter p |
||
| 10 | slib_comment=data can be a list of reals |
||
| 11 | slib_example=[2,4,7],6.1 |
||
| 12 | !exit |
||
| 13 | |||
| 14 | :proc |
||
| 15 | |||
| 16 | !distribute item $wims_read_parm into slib_data,slib_p |
||
| 17 | |||
| 18 | !default slib_N=1000 |
||
| 19 | !if $slib_p<=0 |
||
| 20 | slib_out=0 |
||
| 21 | !exit |
||
| 22 | !endif |
||
| 23 | |||
| 24 | slib_data=!declosing $slib_data |
||
| 25 | slib_cnt=!itemcnt $slib_data |
||
| 26 | slib_out=!exec pari slib_data=[$slib_data];\ |
||
| 27 | print(vector($slib_cnt, x,if(slib_data[x] >= 0, incgam(floor(slib_data[x])+1,$slib_p)/gamma(floor(slib_data[x])+1)))) |
||
| 28 | |||
| 29 | slib_out=!trim $slib_out |