Subversion Repositories wimsdev

Rev

Rev 4330 | 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
4161 bpr 4
 
20 reyssat 5
slib_parms=2\
6
  2,m\
7
  2*m,N\
8
  R, K
7692 bpr 9
 
4330 bpr 10
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
4161 bpr 11
 
20 reyssat 12
slib_example= 5, 100\
13
7,10,Q
3265 bpr 14
slib_require=pari
20 reyssat 15
!exit
16
 
17
:proc
18
!distribute item $wims_read_parm into slib_mm,slib_NN,slib_K
19
 
20
!default slib_mm=2
21
!default slib_NN=$[2*$slib_mm]
22
!default slib_K= R
23
 
24
!if $slib_mm = 1
25
 slib_out =1
26
!else
27
 slib_cc = !values 1/$slib_mm for x = 1 to $[$slib_mm-1]
28
 
29
 !readproc slib/stat/multinomial 1, $[$slib_NN-$slib_mm], $slib_cc
30
 slib_vv= $slib_out
31
 
32
 !if $slib_K issametext Q
7692 bpr 33
  slib_out=!exec pari slib_a=Vec([$slib_vv]);vector($slib_mm, x,(slib_a[x]+1)/$slib_NN)
20 reyssat 34
 !else
35
  slib_out=!values (x+1)/$slib_NN for x in $slib_vv
7692 bpr 36
 !endif
20 reyssat 37
!endif
7692 bpr 38
slib_out =!trim $slib_out