= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 
 
 L u ibU\R`^
 
 
 
 	 	 L u pe<P{ibU\R`^L u I M S L 3 2 . d l l 
 
 
 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 
 
 # L u I M S L O(uf{ 
 
 [ L u I M S L 3 2 . d l l ]   L u I M S L O(uf
 
 
 
         L u I M S L 3 2 . d l l /f N*NL u pe<P{ibU\R`^勓^[I M S L ^-NvRQpeۏLN\ň0
 
 
 
         (WL u I M S L -NvQpe/fǏN~Qpe}T
Tzz I M S L  Qv@b	gQpeGWwQ	g{|<O I M S L : : o d e ( . . . )  v<h_0O(u! u s i n g ( " I M S L " ) ; S{SL u I M S L -NvQpe0
 
 
 
         L u I M S L 3 2 . d l l  I M S L 3 2 . d l l b/fvQNS+TI M S L ^Qpev/ec^	v/ec0I M S L 3 2 . d l l (Wd"}vU_-Nb c[I M S L 3 2 . d l l vMOn0
 
 
 
         YL u I M S L 3 2 . d l l O(uvI M S L 3 2 . d l l b/fvQNS+TI M S L ^Qpev/ec^	^nNNBl(uF o r t r a n R^勨R`^QQpe
TN F _   _4YN/fS|^Qpe0OYL u I M S L /ecQpeF P S 3 H (WI M S L 3 2 . d l l -NQvQpe
T^/fF _ D F P S 3 H NO(uS|^QpeD F P S 3 H ۏL{0
 
 
 
         >f_R}I M S L 3 2 . d l l vel(WR}L u ibU\R`^Rh-NO(uYNellaQSTsSI M S L 3 2 . d l l v_SeN
T_NSNO(uvQNS+TI M S L ^Qpev/ec^0
 
 
 
                 " d l l \ L u I M S L 3 2 . d l l : d l l \ I M S L 3 2 . d l l " 
 
 
 
          NeL u I M S L 3 2 . d l l R}bR\(WS_MRvU_R^bSb _eN I m s l E r r . t x t  (uI S M L Qpev@b	gꁨRQ0R勇eN0L u I M S L 3 2 . d l l (W]\Oe
NNYnzzb RdeN I m s l E r r . t x t  0
 
 
 
         laI M S L ^QpeBl OF o r t r a n <h_vpe~0SNO(uQpes y s : : r e a r r a y (WL u pe~C pe~<h_	TF o r t r a n pe~KNۏLlbc0
 
 
 
         fJTI M S L 3 2 . d l l /f1uF o r t r a n ubvF o r t r a n -Nl	g_8^Uc:g6RelVwSXh NeSuN{/fvc~bkz^0
 
 
 
         L u I M S L 3 2 . d l l \ňvI M S L ^QpeN\Oo:yf[`NKN(uSI M S L ^vfechf}Y0WtُNQpev(ul0I M S L ^/fV i s u a l   N u m e r i c s lQSvNT(WO(uMR{_lQSvO(uS0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : C l e a r I m s l E r r { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : C l e a r I m s l E r r ( ) nzzI S M L QeN I m s l E r r . t x t  vQ[0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : E R S E T { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : E R S E T ( I E R S V R ,   I P A C T ,   I S A C T ) ͑nYtn0
 
 
 
         I E R S V R I M S L ~+R:N1 ^5 peW['Y%N͑0I E R S V R = 0 SbN@b	g&TRNS+TI E R S V R N
Nf%N͑v0
 
         I P A C T I P A C T = - 1 
N9eSS_MRnI P A C T = 0 
NQI P A C T = 1 QI P A C T = 2 ͑n:wn0
 
         I P A C T I P A C T = - 1 
N9eSS_MRnI P A C T = 0 
N\Pbkz^I P A C T = 1 _euSpeI P A C T = 2 ͑n:wn0
 
 
 
         ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : I E R C D { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : I E R C D ( ) ԏVOo`NxeeԏV0 0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : N 1 R T Y { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : N 1 R T Y ( i ) i = 0 bi = 1 ԏVOo`{|W0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : E I { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : E I ( x ) x ! = 0 . 0 v[pe&TR\ԏVЏL0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : S I { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : S I ( x ) x :N[pe0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : E R F { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : E R F ( x ) x :N[pe0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : B S J 0 { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : B S J 0 ( x ) ,{ N{|0 6B e s s e l Qpe0x :N[pe0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : B S J 1 { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : B S J 1 ( x ) ,{ N{|1 6B e s s e l Qpe0x :N[pe0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : B S Y 0 { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : B S Y 0 ( x ) ,{N{|0 6B e s s e l Qpe0x :N[pe0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : B S Y 1 { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : B S Y 1 ( x ) ,{N{|1 6B e s s e l Qpe0x :N[pe0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : B S I 0 { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : B S I 0 ( x ) Ockv,{ N{|0 6B e s s e l Qpe0x :N[pe0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : B S I 1 { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : B S I 1 ( x ) Ockv,{ N{|1 6B e s s e l Qpe0x :N[pe0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : B S K 0 { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : B S K 0 ( x ) Ockv,{	N{|0 6B e s s e l Qpe0x :N[pe0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : B S K 1 { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : B S K 1 ( x ) Ockv,{	N{|1 6B e s s e l Qpe0x :N[pe0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : C S I E Z { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : C S I E Z ( N D A T A ,   X D A T A ,   F D A T A ,   X V E C ,   V A L U E ) c u b i c   s p l i n e   	N!k7hagc<P0
 
 
 
         N D A T A eQvpencppev0
 
         X D A T A ^
N\NN D A T A vpe~0X[>epencpvX t<P0
 
         F D A T A ^
N\NN D A T A vpe~0X[>epencpvY t<P0
 
         X V E C   @bc<PvMOn0:wV A L U E RX V E C /f N*N[pecOSpeV A L U E RV A L U E TX V E C /f$N*Npe~0
 
         V A L U E 
N\NX V E C ^vpe~(uNԏVc<P_0Rv~g0SN:wSpe0
 
 
 
         ԏV<P	gSpeV A L U E RԏVV A L U E &TRԏV N*N[pec<P_0Rv~g	0
 
 
 
         ЏL1 Spe*YYb*Y\2 pencp*Y\bSpe
N9SM3 pe~
NX[(W4 pe~*Y\0
 
 
 
         [ OP[] 
 
 
 
 m a i n ( : x , y , a , b ) = 
 
 { 
 
         x = n e w [ r e a l s , 5 , d a t a :   1 . 6 1 5 ,   1 . 6 3 4 ,   1 . 7 0 2 ,   1 . 8 2 8 ,   1 . 9 2 1 ] , 
 
         y = n e w [ r e a l s , 5 , d a t a :   2 . 4 1 4 5 ,   2 . 4 6 4 5 9 ,   2 . 6 5 2 7 1 ,   3 . 0 3 0 3 5 ,   3 . 3 4 0 6 6 ] , 
 
         o { I M S L : : C S I E Z ( 5 , x , y , 1 . 6 8 2 ) , " \ r \ n " , I M S L : : C S I E Z ( 5 , x , y , 1 . 8 1 3 ) , " \ r \ n " } , 
 
         a = n e w [ r e a l s , 2 , d a t a :   1 . 6 8 2 , 1 . 8 1 3 ] ,   b = n e w [ r e a l s , 2 ] , 
 
         I M S L : : C S I E Z ( 5 , x , y , a , b ) , 
 
         o [ b ( 0 ) , " \ r \ n " , b ( 1 ) , " \ r \ n " ] 
 
 } ; 
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : Q D A G S { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : Q D A G S ( F , A , B , E R R A B S , E R R R E L , & E R R E S T ) { N*NQpevyRS(WzpX[(WGY_	
 
 
 
         F L u  NCQQpeSg0Qpe1u(u7b[IN0
 
         A , B yRNPSyR
NP0
 
         E R R A B S ggv~[|^0E R R A B S > = 0 0S:wSpe:w<P:N0 . 0 0
 
         E R R R E L ggvv[|^0E R R R E L > = 0 0S:wSpe:w<P:N1 e - 6 0
 
         E R R E S T ԏV{v~[]0S:wSpe0
 
         ԏV<PyR<P0
 
 
 
         ЏL1 
NAQL]WY(u2 Spe*YYb*Y\3 c[vh_
NX[(W4 
N/f NCQQpe5 Spe
N&{TBl6 Qpe^ԏV N*N[pe0
 
 
 
         [ OP[]   BlQpef ( x ) = l n ( x ) / s q r t ( x ) (W0 ^1 
NvyR<P
 
 
 
 f ( x ) = l n ( x ) / s q r t ( x ) ; 
 
 I M S L : : Q D A G S [ @ f , 0 . , 1 . ] ;       / / _0RyR<P
 
 ( : e r r ) = I M S L : : Q D A G S [ @ f , 0 . , 1 . , 0 . , 1 e - 6 , & e r r ] , e r r ;     / / _0RyR<Pv~[]
 
 } ; 
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : Q D A G { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : Q D A G ( F , A , B , I R U L E , E R R A B S , E R R R E L , & E R R E S T ) WNG a u s s - K r o n r o d lRO(u^{l{ N*NQpevyR
 
 
 
         F L u  NCQQpeSg0Qpe1u(u7b[IN0
 
         A , B yRNPSyR
NP0
 
         I R U L E yRlR0G a u s s - K r o n r o d lR(uNNbvppeI R U L E <Pppe	1 7 ^1 5 	02 1 1 ^2 1 	03 1 5 ^3 1 	04 2 0 ^4 1 	05 2 5 ^5 1 	06 3 0 ^6 1 	0'YYpeQpecPSI R U L E = 2 YgQpe	g N*NGY_v\<PSI R U L E = 1 YgQpe/fFdRvSI R U L E = 6 0S:wSpe:w<P:N2 0
 
         E R R A B S ggv~[|^0E R R A B S > = 0 0S:wSpe:w<P:N0 . 0 0
 
         E R R R E L ggvv[|^0E R R R E L > = 0 0S:wSpe:w<P:N1 e - 6 0
 
         E R R E S T ԏV{v~[]0S:wSpe0
 
         ԏV<PyR<P0
 
 
 
         ЏL1 
NAQL]WY(u2 Spe*YYb*Y\3 c[vh_
NX[(W4 
N/f NCQQpe5 Spe
N&{TBl6 Qpe^ԏV N*N[pe0
 
 
 
         [ OP[]   BlQpef ( x ) = l n ( x ) / s q r t ( x ) (W0 ^1 
NvyR<P
 
 
 
 f ( x ) = l n ( x ) / s q r t ( x ) ; 
 
 I M S L : : Q D A G [ @ f , 0 . , 1 . ] ;       / / _0RyR<P
 
 ( : e r r ) = I M S L : : Q D A G [ @ f , 0 . , 1 . , 2 , 0 . , 1 e - 6 , & e r r ] , e r r ;     / / _0RyR<Pv~[]
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : Q D A G P { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : Q D A G P ( F , A , B , P T S , E R R A B S , E R R R E L , & E R R E S T ) { N*N~[GY_pvQpevyR
 
 
 
         F L u  NCQQpeSg0Qpe1u(u7b[IN0
 
         A , B yRNPSyR
NP0
 
         P T S pe~X[>eyR:Svep\	g N*Nep08^ُNep/fyQpeQsGY_<Pv0We0
 
         E R R A B S ggv~[|^0E R R A B S > = 0 0S:wSpe:w<P:N0 . 0 0
 
         E R R R E L ggvv[|^0E R R R E L > = 0 0S:wSpe:w<P:N1 e - 6 0
 
         E R R E S T ԏV{v~[]0S:wSpe0
 
         ԏV<PyR<P0
 
 
 
         ЏL1 
NAQL]WY(u2 Spe*YYb*Y\3 c[vh_
NX[(W4 
N/f NCQQpe5 Spe
N&{TBl6 Qpe^ԏV N*N[pe0
 
 
 
         [ OP[]   BlQpef ( x ) = x ^ 3 * l n { a b s [ ( x ^ 2 - 1 ) * ( x ^ 2 - 2 ) ] } (W0 ^3 
NvyR<P
 
 
 
 f ( x ) = x ^ 3 * l n { a b s [ ( x ^ 2 - 1 ) * ( x ^ 2 - 2 ) ] } ; 
 
 I M S L : : Q D A G P [ @ f , 0 . , 3 . , m a t h : : r a 1 ( 1 , s q r t ( 2 . 0 ) ) ] ; 
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : Q D A G I { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : Q D A G I ( F , A , B , E R R A B S , E R R R E L , & E R R E S T ) { N*NQpe(Wewz:SbJSewz:SvyR
 
 
 
         F L u  NCQQpeSg0Qpe1u(u7b[IN0
 
         A , B yRNPSyR
NP0x /f N*NpeRSSv<P:N  A = - i n f , B = x - "x 	A = x , B = i n f x + "	A = - i n f , B = i n f - "+ "	0
 
         E R R A B S ggv~[|^0E R R A B S > = 0 0S:wSpe:w<P:N0 . 0 0
 
         E R R R E L ggvv[|^0E R R R E L > = 0 0S:wSpe:w<P:N1 e - 6 0
 
         E R R E S T ԏV{v~[]0S:wSpe0
 
         ԏV<PyR<P0
 
 
 
         ЏL1 
NAQL]WY(u2 Spe*YYb*Y\3 c[vh_
NX[(W4 
N/f NCQQpe5 Spe
N&{TBl6 Qpe^ԏV N*N[pe0
 
 
 
         [ OP[]   BlQpex * e x p [ - x ] (W0 + "	
NvyR<P0
 
 
 
 f ( x ) = x * e x p [ - x ] ; 
 
 I M S L : : Q D A G I [ @ f , 0 . 0 , i n f ] ; 
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : T W O D Q { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : T W O D Q ( F ,   A ,   B ,   G ,   H ,   I R U L E ,   E R R A B S ,   E R R R E L ,   & E R R E S T ) {N͑yR
 
 
 
         F L u NCQyQpeSg0Qpe1u(u7b[IN0
 
         A , B yRNPNyR
NP0
 
         G , H L u  NCQQpeSg1u(u7b[IN(ueg{QB\yRvNPN
NP0
 
         I R U L E 	byRvlR0G a u s s - K r o n r o d lR(uNNbvppeI R U L E <Pppe	1 7 ^1 5 	02 1 1 ^2 1 	03 1 5 ^3 1 	04 2 0 ^4 1 	05 2 5 ^5 1 	06 3 0 ^6 1 	0'YYpeQpecPSI R U L E = 2 YgQpe	g N*NGY_v\<PSI R U L E = 1 YgQpe/fFdRvSI R U L E = 6 0S:wSpe:w<P:N2 0
 
         E R R A B S ggv~[|^0E R R A B S > = 0 0S:wSpe:w<P:N0 . 0 0
 
         E R R R E L ggvv[|^0E R R R E L > = 0 0S:wSpe:w<P:N1 e - 6 0
 
         E R R E S T ԏV{v~[]0S:wSpe0
 
         ԏV<PyR<P0
 
 
 
         ЏL1 
NAQL]WY(u2 Spe*YYb*Y\3 c[vh_
NX[(W4 ꁚ[INQpeSpe*Npe
N&{TBl5 Spe
N&{TBl6 Qpe^ԏV N*N[pe0
 
 
 
         [ OP[]   {Qpef ( x , y ) = y * c o s ( x + y * y ) vя<OyR<P0YB\:SS[ 0 , 1 ] QB\:SS[ - 2 * x , 5 * x ] 0
 
 
 
 f ( x , y ) = y * c o s ( x + y * y ) ; 
 
 g ( x ) = - 2 * x ; 
 
 h ( x ) = 5 * x ; 
 
 I M S L : : T W O D Q [ @ f , 0 . , 1 . , @ g , @ h , 6 ] ; 
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : Q A N D { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : Q A N D ( F ,   m a t h : : r a 1 [ A 1 ,   A 2 ,   . . . ,   A n ,   B 1 ,   B 2 ,   . . . ,   B n ] ,   M A X F C N ,   E R R A B S ,   E R R R E L ,   & E R R E S T ) {Qpe(Wwb_NvyR
 
 
 
         F L u YCQQpeSg	gn *NSϑ0Qpe1u(u7b[IN0
 
         m a t h : : r a 1 [ A 1 ,   A 2 ,   . . . ,   A n ,   B 1 ,   B 2 ,   . . . ,   B n ] pe~^:Nn + n 0A 1 ,   A 2 ,   . . . ,   A n /fyRNP	gn *NNP<P0B 1 ,   B 2 ,   . . . ,   B n /fyR
NP	gn *N
NP<P0
 
         M A X F C N SAQvQpe0Ov g'Ypev0_{M A X F C N < = 2 5 6 ^ n 0S:wSpe:w<P:N6 4 ^ n 0
 
         E R R A B S ggv~[|^0E R R A B S > = 0 0S:wSpe:w<P:N0 . 0 0
 
         E R R R E L ggvv[|^0E R R R E L > = 0 0S:wSpe:w<P:N1 e - 6 0
 
         E R R E S T ԏV{v~[]0S:wSpe0
 
         ԏV<PyR<P0
 
 
 
         ЏL1 
NAQL]WY(u2 Spe*YYb*Y\3 c[vh_
NX[(W4 ꁚ[INQpe/feSQpe5 Spe
N&{TBl  6 Qpe^ԏV N*N[pe0
 
 
 
         [ OP[]   {Qpee x p [ - ( x 1 ^ 2 + x 2 ^ 2 + x 3 ^ 2 ) ] (WibU\vzeSOte*N	N~zz	
NvyRvя<O<P0yR:S:N[ - 2 , - 2 , - 2 ^2 , 2 , 2 ] 0
 
 
 
 f ( x 1 , x 2 , x 3 ) = e x p [ - ( x 1 ^ 2 + x 2 ^ 2 + x 3 ^ 2 ) ] ; 
 
 I M S L : : Q A N D [ @ f , m a t h : : r a 1 ( - 2 , - 2 , - 2 , 2 , 2 , 2 ) ] ; 
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : D E R I V { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : D E R I V ( F C N , K O R D E R , X , B G S T E P , T O L ) { NCQQpev N60N6b	N6[pe
 
 
 
         F C N L u  NCQQpeSg{Qpe(WX v[pe0Qpe1u(u7b[IN0
 
         K O R D E R [pev61 02 b3 	0
 
         X Bl[vp0
 
         B G S T E P (ueg{ekvwY<Pُ*Nek(uN{[pe0B G S T E P > 0 0S:wSpe:w<P:N0 . 1 0
 
         T O L AQvv[]0 ^1 	0S:wSpe:w<P:N1 e - 6 0
 
         ԏV<PQpev N60N6b	N6[pe0
 
 
 
         ЏL1 
NAQL]WY(u2 Spe*YYb*Y\3 c[vh_
NX[(W4 
N/f NCQQpe5 Spe
N&{TBl  6 Qpe^ԏV N*N[pe0
 
 
 
         [ OP[]   BlQpef ( x ) = 2 * x ^ 4 + 3 * x (W0 . 7 5 Yv	N6[pe
 
 
 
 f ( x ) = 2 * x ^ 4 + 3 * x ; 
 
 I M S L : : D E R I V [ @ f , 3 , 0 . 7 5 ] ; 
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : o d e { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : o d e ( F C N ,   m a t h : : r a 1 [ t 0 ,   t 1 ,   . . . ,   t k ] ,   m a t h : : r a 1 [ Y 1 , Y 2 , . . . , Y n ] ,   T O L ) Bl N!k8^_Rez[ y ' = f ( t , y )   , y ( t 0 ) = y 0 ] R u n g e - K u t t a el
 
 
 
         F C N (u7b[INvQpe(uN{_Rez~-NTezSzQpe<P1u(u7b0h_	g2 * n + 1 *NSpe,{ N*NSpe:NSϑTn *NSpe:NQpe<P gTn *NSpe:NSzQpe<PsS_Rezv<P	0n :N_Rez~-Nezv*Npe_N/f*gwQpev*Npe0b__YN
 
 
 
 f ( t , y 1 , y 2 , y 3 , d 1 , d 2 , d 3 ) = 
 
 {       d 1 = y 2 , 
 
         d 2 = - y 1 , 
 
         d 3 = - y 3 
 
 } ; 
 
 
 
         m a t h : : r a 1 [ t 0 ,   t 1 ,   . . . ,   t k ] pe~X[>erzSϑt 0eQRYe:Nt 0 (u7b^g_0Ret 1 ,   . . . ,   t k ev0
 
         m a t h : : r a 1 [ Y 1 , Y 2 , . . . , Y n ] X[>en *N*gwQpe(WwYpYvQpe<PY ( n ) = Y n ( t ) 0
 
         T O L AQ]0S:wSpe:w<P:N1 e - 6 0
 
 
 
         ԏV<P( 1 + k )  ( 1 + n ) N~pe~0X[>e( 1 + k ) ~~gS+TRY<P	k N~~gsSt , Y 1 , Y 2 , . . . , Y n 0
 
 
 
         ЏL1 
NAQL]WY(u2 Spe*YYb*Y\3 c[vh_
NX[(W4 QpeSpe*Npe
N&{TBl5 Spe
N&{TBl  6 QX[7 Qpe^ԏV N*N[pe0
 
 
 
         [ OP[]    N6_Rez~SR<P:N
 
                 r ' = 2 r - 2 r f ,   r ( 0 ) = 1 
 
                 f ' = - f + r f ,     f ( 0 ) = 3 
 
         {t = 1 , 2 , . . . , 1 0 evr 0f v<P0
 
 
 
         z^YN
 
 
 
 ! ! ! u s i n g [ " I M S L " , " m a t h " ] ; 
 
 f ( t , r , f , d r , d f ) = { d r = 2 * r - 2 * r * f ,   d f = - f + r * f } ;   / / Qpe[IN
 
 o d e [ @ f , r a 1 [ 0 , 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 , 1 0 ] , r a 1 [ 1 , 3 ] ] . o u t a [ ] ; 
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : M A C H { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : M A C H ( i ) ԏVS|^[W8^pe  0x :Ntepe0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 # I M S L : : I F N A N { 
 
 [ L u I M S L 3 2 . d l l ]   I M S L : : I F N A N ( x ) $Rex /f&T/fN A N N A N h:y
N/f N*Npe  0x :N[pe0ЏL1 Spe
N&{TBl0
 
 } # 
 
 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 
 
 