Subversion Repositories wimsdev

Rev

Rev 4349 | Go to most recent revision | 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
4158 bpr 4
 
4349 bpr 5
slib_author=Gang, XIAO
4158 bpr 6
 
20 reyssat 7
slib_example=sin(x),x,0,5,numeric\
7116 bpr 8
sin(x),x,0,5\
9
sqrt(sin(x+4)^2 + cos(x+0.5)^2),x,0,5,numeric\
20 reyssat 10
sin(x),x
7116 bpr 11
 
3265 bpr 12
slib_require=maxima, pari
20 reyssat 13
!exit
14
 
15
:proc
16
!distribute items $wims_read_parm into slib_fn,slib_var,slib_x1,slib_x2,slib_opt
17
slib_fn=!rawmath $slib_fn
18
slib_var=!varlist $slib_var
19
slib_var=!item 1 of $slib_var
20
!default slib_var=x
21
slib_out=
22
!if $slib_x1=$slib_empty or $slib_x2=$slib_empty
23
 slib_out=!exec maxima integrate($slib_fn,$slib_var);
24
 !exit
25
!endif
26
slib_test=$[($slib_x1)+($slib_x2)]
27
!if numeric notin $slib_opt
28
 slib_out=!exec maxima integrate($slib_fn,$slib_var,$slib_x1,$slib_x2);
29
!endif
30
!if ($slib_out=$empty or integrate isin $slib_out or i isvarof $slib_out) \
31
  and NaN notin $slib_test
7116 bpr 32
   !!no more a true option (it was not documented), default is fast !
33
 !!!if fast isin $slib_opt
20 reyssat 34
  slib_out=!exec pari intnum($slib_var=$slib_x1,$slib_x2,$slib_fn)
7116 bpr 35
 !!!else
36
 !! slib_out=!exec pari intnum($slib_var=$slib_x1,$slib_x2,$slib_fn)
37
 !!!endif
20 reyssat 38
!endif
39