Subversion Repositories wimsdev

Rev

Rev 4328 | Rev 4330 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 4328 Rev 4329
Line 2... Line 2...
2 
!if $wims_read_parm!=slib_header
 2 
!if $wims_read_parm!=slib_header
3 
!goto proc
 3 
!goto proc
4 
!endif
 4 
!endif
5 
 
 5 
 
6 
slib_author=Bernadette PERRIN-RIOU
 6 
slib_author=Bernadette PERRIN-RIOU
- 
 
 7 
 
- 
 
 8 
slib_parms=5\
- 
 
 9 
,[data]\
- 
 
 10 
0,mu\
- 
 
 11 
1,sigma
7 
 
 12 
 
8 
slib_example=[0.2,0.4,0.5],0,1
 13 
slib_example=[0.2,0.4,0.5],0,1
9 
slib_require=octave
 14 
slib_require=octave
10 
!exit
 15 
!exit
11 
 
 16 
 
12 
:proc
 17 
:proc
13 
 
 18 
 
14 
!distribute  item $wims_read_parm into slib_data,slib_mean,slib_var,slib_x1,slib_x2
 19 
!distribute  item $wims_read_parm into slib_data,slib_mean,slib_var
15 
slib_out=
 20 
slib_out=
16 
!default slib_mean=0
 21 
!default slib_mean=0
17 
!default slib_var=1
 22 
!default slib_var=1
18 
!default slib_x1=$[($slib_mean-5)*($slib_var)]
 - 
 
19 
!default slib_x2=$[($slib_mean+5)*($slib_var)]
 - 
 
- 
 
 23 
 
20 
slib_data=!declosing $slib_data
 24 
slib_data=!declosing $slib_data
21 
 
 25 
 
22 
slib_out=!exec octave normal_inv([$slib_data], $slib_mean,$slib_var^2)
 26 
slib_out=!exec octave normal_inv([$slib_data], $slib_mean,$slib_var^2)
23 
!!slib_out=!exec octave norminv([$slib_data], $slib_mean,$slib_var^2)
 27 
!!slib_out=!exec octave norminv([$slib_data], $slib_mean,$slib_var)
24 
 
 28 
 
25 
 
 29 
 
26 
slib_out=!words2items $slib_out
 30 
slib_out=!words2items $slib_out
27 
 
 31 
 
28 
slib_out=!trim $slib_out
 32 
slib_out=!trim $slib_out