Subversion Repositories wimsdev

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

slib_author=Bernadette, Perrin-Riou

slib_example=x^2+sin(x*y)+exp(y^2),[x,y],0,3,-3,2,50

!exit

:proc
!distribute items $wims_read_parm into slib_fn,slib_var,slib_x1,slib_x2,slib_y1,slib_y2,slib_steps
slib_fn=!rawmath $slib_fn
slib_var=!declosing $slib_var
slib_var1=!item 1 of $slib_var
slib_var2=!item 2 of $slib_var
!default slib_var1=x
!default slib_var2=y
slib_x1=$[$slib_x1]
!if $slib_x1=$empty or $slib_x1 isitemof NaN,Inf
  slib_x1=0
!endif
slib_x2=$[$slib_x2]
!if $slib_x2=$empty or $slib_x2 isitemof NaN,Inf
        slib_x2=1
!endif
slib_y1=$[$slib_y1]
!if $slib_y1=$empty or $slib_y1 isitemof NaN,Inf
  slib_y1=0
!endif
slib_y2=$[$slib_y2]
!if $slib_y2=$empty or $slib_y2 isitemof NaN,Inf
 slib_x2=1
!endif
!bound slib_steps between integer 1 and 1000 default 50
!bound slib_steps2 between integer 1 and 1000 default 50
slib_step=$[($slib_x2-($slib_x1))/$slib_steps]
slib_step2=$[($slib_y2-($slib_y1))/$slib_steps]
slib_min=$empty
slib_max=$empty

!for slib_i=0 to $slib_steps
        slib_fxn=!mathsubst $slib_var1= $[$slib_x1+$slib_i*$slib_step] in $slib_fn
        slib_fxn=!replace item NaN by $ in $slib_fxn
        slib_vals=!values $slib_fxn for $slib_var2=$slib_y1 to $slib_y2 step $slib_step
        slib_vals=!replace item NaN by $ in $slib_vals
        slib_vals=!nonempty items $slib_vals
        !if $slib_i>1
          slib_vals=!sort numeric items $slib_min,$slib_vals,$slib_max
        !else
                slib_vals=!sort numeric items $slib_vals
        !endif
        slib_vals=!trim $slib_vals
        slib_min=!item 1 of $slib_vals
        slib_max=!item -1 of $slib_vals
!next slib_i

slib_out=$slib_min,$slib_max