Subversion Repositories wimsdev

Rev

Rev 4330 | Blame | Compare with Previous | Last modification | View Log | RSS feed

!if $wims_read_parm!=slib_header
  !goto proc
!endif

slib_parms=4\
 ,n\
 ,a\
 ,b\
 ,m
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou

slib_example=10, 3,7,11
slib_require=pari
!exit

:proc
slib_out=

!distribute item $wims_read_parm into slib_n,slib_a,slib_b,slib_m
slib_A=!randint $slib_m

slib_out=!exec pari {slib_L=vector($slib_n);slib_A=$slib_A;\
for(slib_k=1,$slib_n,\
slib_A=lift(Mod($slib_a*slib_A+$slib_b,$slib_m));\
slib_L[slib_k]=slib_A);print(slib_L);}

slib_out =!trim $slib_out