Subversion Repositories wimsdev

Rev

Rev 4328 | 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
 
4330 bpr 5
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
4161 bpr 6
 
4328 bpr 7
slib_parms=3\
8
,[data]\
4330 bpr 9
1,mu\
4328 bpr 10
1,sigma
11
 
12
 slib_example=[2,4,7] ,3,4\
13
 2,2,sqrt(2)\
14
 2,2,2\
15
 2,1,1
20 reyssat 16
 !exit
17
 
18
:proc
19
slib_out=
20
 
4328 bpr 21
!distribute item $wims_read_parm into slib_data,slib_mean,slib_var
20 reyssat 22
!default slib_mean=1
23
!default slib_var=1
4328 bpr 24
!! en fait c'est un écart type
20 reyssat 25
slib_data=!declosing $slib_data
26
 
27
slib_out=
28
!for slib_i in  $slib_data
29
 !if $slib_i > 0
30
   slib_out = !append item $[exp(-(log($slib_i)-$slib_mean)^2/(2*$slib_var^2))/($slib_i*$slib_var*(2*pi)^0.5)] to $slib_out
31
 !else
32
   slib_out = !append item 0 to $slib_out
33
 !endif
34
!next slib_i
35
 
36
 
37
slib_out=!trim $slib_out