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 |