Subversion Repositories wimsdev

Rev

Rev 7692 | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
20 reyssat 1
!if $wims_read_parm!=slib_header
13586 bpr 2
  !goto proc
20 reyssat 3
!endif
4161 bpr 4
 
20 reyssat 5
slib_parms=2\
6
,[data]\
7
1,p
8
 
4330 bpr 9
slib_author=Sophie, Lemaire; Bernadette, Perrin-Riou
4161 bpr 10
 
20 reyssat 11
slib_example=[1,2,3],7
12
 
13
!exit
14
 
15
:proc
16
 
17
!distribute  item $wims_read_parm into slib_data,slib_p
18
 
19
!default slib_N=1000
4330 bpr 20
!if $slib_p<=0
13586 bpr 21
  slib_out=0
22
  !exit
20 reyssat 23
!endif
24
 
25
slib_data=!declosing $slib_data
26
slib_out=
27
!for slib_x in $slib_data
13586 bpr 28
  !if $slib_x<0
29
    slib_out=!append item 0 to $slib_out
30
  !else
31
    slib_out=!append item  $[$slib_p *exp(-$slib_p*$slib_x)] to $slib_out
32
  !endif
20 reyssat 33
!next slib_x
34
slib_out=!trim $slib_out