1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
|
!if $wims_read_parm!=slib_header
!goto proc
!endif
slib_parms=2\
1,m\
1,lambda
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
slib_example= 5,2.5
slib_require=octave
!exit
:proc
!distribute item $wims_read_parm into slib_data, slib_lambda
slib_out=$empty
!default slib_m=1
!default slib_lambda=1
!!slib_ex=$[exp(-$slib_lambda)]
!!slib_q=$slib_ex
!!
!!
!!slib_U=!random 0,1 repeat $slib_m
!!
!!!for slib_i=1 to $slib_m
!! slib_V=!item $slib_i of $slib_U
!! slib_q=$slib_ex
!! slib_p=$slib_ex
!! slib_j=0
!!
!!!while $slib_V >= $slib_q
!! !advance slib_j
!! slib_p=$[$slib_p*$slib_a/$slib_j]
!! slib_q=$[$slib_q+$slib_p]
!!!endwhile
!!
!!slib_out=!append item $slib_j to $slib_out
!!!next slib_i
slib_data=!declosing $slib_data
slib_cnt=!itemcnt $slib_data
!if $slib_cnt=1
slib_data=1,$slib_data
!endif
!!deprecated syntax
!!slib_out=!exec octave poisson_rnd($slib_lambda,$slib_data)
slib_out=!exec octave poissrnd($slib_lambda,$slib_data)
slib_out=!trim $slib_out
slib_out=!replace internal $\
$ by ; in $slib_out
slib_out=!words2items $slib_out
|
