Subversion Repositories wimsdev

!if $wims_read_parm!=slib_header

 !goto proc

!endif

slib_title=Generation of a discrete law with positive coefficients

slib_parms=2\

  2,m\

  2*m,N\

  R, K

slib_author=S. Lemaire and Bernadette Perrin-Riou

slib_out=random  data of m>=2 positive  real numbers  p_1,p_2,..p_m of

sum 1 such that p_i*N are integers, (N>m). Outputs\

are real (K=R) or rational (K=Q)

slib_example= 5, 100\

7,10,Q

slib_require=pari

!exit

:proc

!distribute item $wims_read_parm into slib_mm,slib_NN,slib_K

!default slib_mm=2

!default slib_NN=$[2*$slib_mm]

!default slib_K= R

!if $slib_mm = 1

 slib_out =1

!else

 slib_cc = !values 1/$slib_mm for x = 1 to $[$slib_mm-1]

 !readproc slib/stat/multinomial 1, $[$slib_NN-$slib_mm], $slib_cc

 slib_vv= $slib_out

 !if $slib_K issametext Q

  slib_out=!exec pari slib_a=Vec([$slib_vv]);vector($slib_mm, x,(slib_a[x]+1)/$slib_NN) 

 !else

  slib_out=!values (x+1)/$slib_NN for x in $slib_vv

 !endif 

!endif

slib_out =!trim $slib_out