+
    i                        ^ RI HtHtHtHtHtHt ^ RIHtH	t	 ^ RI
Ht ^ RIHtHtHt ^ RIHt ^ RIHtHt ^ RIHtHt ^ RIHtHt ^ R	IHt ^ R
IHtHt ^ RI H!t! ^ RI"H#t# ^ RI$H%t%H&t&H't'H(t( ^ RI)H*t* ^ RI+H,t,H-t- ^ RI.H/t/ ^ RI0H1t1 ^ RI2H3t3 ^ RI4H5t5 ^ RI6H7t7 ^ RI8H9t9 R t:R t;R t<R t=R t>R t?R t@R tAR tBR tCR  tDR! tER" tFR# tGR$ tHR% tIR& tJR' tKR( tLR) tMR* tNR+ tOR, tPR- tQR. tRR/ tSR0 tTR1 tUR2 tVR3# )4    )DifferentialOperatorHolonomicFunctionDifferentialOperators
from_hyperfrom_meijergexpr_to_holonomic)RecurrenceOperatorsHolonomicSequence)
EulerGamma)IRationalpi)S)Symbolsymbols)explog)asinhcosh)sqrt)cossin)besselj)beta)CiSierferfc)gamma)hypermeijerg)sstr)O)hyperexpand)ZZ)QQ)RRc            	         \        R 4      p \        \        P                  ! V 4      R4      w  rW!P                  8X  g   Q hV\        VP                  P                  VP                  P                  .V4      8X  g   Q hW,          V ^,          V^,          ,          ,           \        ^ W ^,          .V4      8X  g   Q hV ^,          ^,           V,           V V^,          ,          ,           \        V ^,          ^,           ^^ ^ ^ V .V4      8X  g   Q hW,          V ^,          ,           ^,           W ^,          V ,           ,          ,
          ^,          RV ^,          ,          RV ^,          ,          V,          ,           RV ^,          ,          V^,          ,          ,           V ^	,          ) V^,          ,          ,           8X  g   Q hW^,          ,          V ^,          ^,           V^,          ,          ,           W ^,          ,           ,          p^V ,          ^V ^,          ,          V,          ,           V ^,          V^,          ,          ,           ^V ^,          ,          V ,           ^<,           V^,          ,          ,           ^
V ^,          ,          ^V ,          ,           V^,          ,          ,           V ^,          ^V ^,          ,          ,           V^,          ,          ,           V ^,          ^,           V^,          ,          ,           pW48X  g   Q hR# )xDxNiii)	r   r   r&   old_poly_ringderivative_operatorr   basezerooner)   Rr*   pqs        b/var/www/html/photoedit/myenv/lib/python3.14/site-packages/sympy/holonomic/tests/test_holonomic.pytest_DifferentialOperatorr5      sM   A!""2"21"5t<EA&&&&&%qvv{{AFFJJ&?CCCC6AqD2q5L $8!Q1q$IIIIqD1H?Q
A %q!taxAq!Q&?CD D DFQTMAdQh/!3ad
	q!tr8 1a4Z2q508156TERU?8C C C C	
UadQh"a%'	'BAI6A	
Q1q!t8r/!QTRUN2	adQ	r1u$	%(*QT	BF(:b!e'C	D	
AAqDBE!	"%&TAXQ$6	7A 6M6    c                     \        R 4      p \        \        P                  ! V 4      R4      w  r\	        V^,          V ,          V 4      p\	        ^V,          W^,          ,          ,           V 4      pW48X  g   Q h\	        W,          ^,           V 4      p\	        V^,           V 4      p\	        V ^,
          V ^,          ^,
          V,          ,           V ^,          V ,
          V^,          ,          ,           V 4      pW4,           V8X  g   Q h\	        W,          V^,          V ^,          ^,           ,          ,           V 4      p\	        V^,
          V 4      p\	        RV ^,          ,          ^~V ,          ,
          ^,
          RV ^,          ,          ^V ^,          ,          ,
          RV ,          ,
          ^,           V,          ,           R	V ^,          ,          ^V ^,          ,          ,
          ^V ,          ,           ^,
          V^,          ,          ,           ^	V ^,          ,          ^V ^,          ,          ,           ^&V ^,          ,          ,           ^V ,          ,           ^(,           V^,          ,          ,           V 4      pW4,           V8X  g   Q h\	        V^,          ^,
          V 4      p\	        V ^,          V,           V 4      p\	        V ^,          ) ^-V ^,          ,          ,           RV ^
,          ,          ,
          RV ^,          ,          ,           V ^,          ,
          RV ^,          ,          ,
          V ^,          ) ^V ^,          ,          ,           ^V ^,          ,          ,
          ^V ^,          ,          ,           ^,
          V,          ,           V ^,          ^-V ^,          ,          ,
          RV ^
,          ,          ,           RV ^,          ,          ,
          V ^,          ,           RV ^,          ,          ,           V^,          ,          ,           V ^,          ^V ^,          ,          ,
          ^V ^,          ,          ,           ^V ^,          ,          ,
          ^,           V^,          ,          ,           V 4      pW4,           V8X  g   Q hV ^,          ^V ,          ,           ^,           pV ^,          ^V ,          ,
          ^,           pW2,          VP                  4       ,
          pWB,          VP                  4       ,
          p\	        W04      \	        W@4      ,           p\	        ^V ^,          ,          ^V ,          ,           ^,           R
V ^,          ,          ^V ^,          ,          ,
          ^>V ,          ,
          ^
,           V,          ,           V ^,          ^V ^,          ,          ,           ^V ^,          ,          ,           ^
V ,          ,
          ^G,
          V^,          ,          ,           V 4      pWV8X  g   Q hR# )r)   r*   i  i  i'  iv  Niiyi)r   r   r%   r+   r   diff)r)   r1   r*   r2   r3   rss          r4   test_HolonomicFunction_additionr<   +   sG   A!""2"21"5t<EA"a%!)Q'A1(aq5[0!4A6M6!&1*a(A"q&!$A1q5QTAXO3q!tax2q56HH!LA5A::!&2q5AqD1H#55q9A"q&!$A3A:a/#5$A+aQRd
:RUX[\U\:\_b:bfh9hh1*rAqDy(261C72q5@Aad(R!Q$Y&ad2R!V;R?2q5HIJK	MA 5A::"a%!)Q'A!Q$)Q'AArE6Bq"uH,s1b5y8419Dq!tKcRSUVRVhV
R%"QU(	SAX	%AqD	01	4b89<=rEBq"uH<LsSTVXSXy<XQT	=qD=q!t8=$%'U<++./ebBh.>QT.ICPQSTPTH.T	/
q5. 	A 3!8O8	1qs
QA	1qs
QA	qvvxA	qvvxA!"3A"99A1QT6BqD=2-"QT'Bq!tG2Cbd2JR2OQS1SS	
A!Q$AqD	 2a4	'"	,b!e3456	8A6M6r6   c                  &   \        R 4      p \        \        P                  ! V 4      R4      w  r\	        W ,           W^,          ,          ,           V 4      p\	        W,          W ,          ,           V^,          ,           V 4      p\	        ^V ^,          ,          ^V ^,          ,          ,           ^V ^,          ,          ,           ^,           ^V ^,          ,          ^V ^,          ,          ,
          ^V ,          ,           V,          ,           ^V ^,          ,          ^V ^,          ,          ,           ^V ^,          ,          ,           ^,           V^,          ,          ,           ^V ^,          ,          ^V ^,          ,          ,           ^V ,          ,           V^,          ,          ,           ^V ^,          ,          V ^,          ,           V^,          ,          ,           V 4      pW4,          V8X  g   Q h\	        V^,          ^,           V 4      p\	        V^,
          V 4      p\	        ^RV,          ,           ^V^,          ,          ,           V 4      pW4,          V8X  g   Q h\	        V^,          ^,           V ,           V,           V 4      p\	        W ,          ^,
          ^,          V 4      p\	        ^V ^,          ,          ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,
          ^V ^,          ,          ,
          ^V ,          ,
          ^,
          ^V ^,          ,          ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,
          ^V ^,          ,          ,
          ^V ,          ,
          ^,
          V,          ,           ^V ^,          ,          ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,
          ^V ^,          ,          ,
          ^V ,          ,
          ^,
          V^,          ,          ,           ^V ^,          ,          ^
V ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,
          ^V ,          ,
          V^,          ,          ,           ^V ^,          ,          ^V ^,          ,          ,           V ^,          ,
          V^,          ,          ,           V 4      pW4,          V8X  g   Q h\	        W^,          ,          ^,
          V 4      p\	        W ,          V ,
          V 4      p\	        V ^,
          RV ,          ^,           V,          ,           W^,          ,          ,           V 4      pW4,          V8X  g   Q hR# )r)   r*   N)r   r   r%   r+   r   r)   r1   r*   r2   r3   r:   s         r4   %test_HolonomicFunction_multiplicationr@   M   s   A!""2"21"5t<EA"$qQw,*A!$rt)BE/1-A1QT6AadF?Qq!tV3a7Bq!tGa1f<LrRSt<SUW;WW	
1a4"QT'	Bq!tG	#a	'Q./23AqD&1QT6/AaC2GQ1NO	
1a4!Q$A 	"A 3!8O8"a%'1%A"Q$"A1R-1b!e)3Q7A3!8O8"a%'!)B,*A246A+q)A1QT6Bq!tG+bAg5!Q$>1a4G!AqD&PSTUVSVVYZZ	
1a4"QT'	Bq!tG	#a1f	,r!Q$w	61	<q	@"DE	
1a4"QT'	Bq!tG	#a1f	,qAv	5!	;a	?QFGJKAqD&q!tGK1fK AvK&()!K,-/UJ33 781fqAvo16LbRSe5ST VW	XA 3!8O8!E'!)Q'A"$q&!$A1q5RTAXrM1QAI=qAA3!8O8r6   c                     \        R 4      p \        \        P                  ! V 4      R4      w  r\	        W ,           W^,          ,          ,           V 4      p\	        W 4      p\        ^
4       F  pWCV,          8X  g   Q hWC,          pK  	  R# )r)   r*   N)r   r   r%   r+   r   range)r)   r1   r*   r2   ans         r4   test_HolonomicFunction_powerrE   g   si    A!""2"21"5t<EA"$qQw,*A" A2YqDyy	 r6   c            
         \        R 4      p \        \        P                  ! V 4      R4      w  r\	        V^,
          V ^ ^.4      p\	        V^,          ^,           V ^ ^^ .4      p\	        RV,           V^,          ,
          V^,          ,           V ^ . RO4      pW4,           V8X  g   Q h\	        W ,
          V^,          ,           V ^ ^^.4      p\	        V^,          V ,           V ^ ^^ .4      p\	        V ^,          ) V ^,          ^,          ,
          V ^,          ,
          \        ^^4      ,           V ^,          V ^,          ^,          ,           V \        ^^4      ,          ,           ^,           V,          ,           V \        R^4      ,          \        ^^4      ,           V^,          ,          ,           V ^,          V \        ^^4      ,          ,
          \        ^^4      ,           V^,          ,          ,           V ^,          V ^,          ,           \        P                  ,           V^,          ,          ,           V ^ . RO4      pW4,           V8X  g   Q h\	        V^,          ^V ,          V,          ,           V ^,          ,           V ^ ^^.4      p\	        V^,          ^,           V ^ ^^.4      p\	        V ^,          ^V ^,          ,          ,           ^V ^,          ,          ,
          ^,
          ^V ^,          ,          ^$V ^,          ,          ,           ^ V ,          ,
          V,          ,           V ^,          ^V ^,          ,          ,           ^V ^,          ,          ,           ^	,
          V^,          ,          ,           ^V ^,          ,          ^$V ^,          ,          ,           ^ V ,          ,
          V^,          ,          ,           V ^,          ^
V ^,          ,          ,           ^,
          V^,          ,          ,           V ^ . R	O4      pW4,           V8X  g   Q h\	        V^,          V ,           V ^. R
O4      p\	        V^,
          V ^^.4      p\	        V ^,          ) V ,
          ^,           V ^,          V ,           V,          ,           V ) ^,
          V^,          ,          ,           V ^,           V^,          ,          ,           V ^. RO4      pW4,           V8X  g   Q h\        \        V 4      4      p\        ^V ,          ^R7      p\	        V ^,          ^,           V ^,          ^V ,          ,           V,          ,           V ^,          ^,           V^,          ,          ,           V ^,          ^V ,          ,           V^,          ,          ,           V ^\        ^4      ^,           R\        ^4      ,           \        ^4      ) ^,           .4      pW4,           V8X  g   Q h\        R4      p\        \        V 4      4      p\        \        V ^,          V ,
          4      4      pW4,           P                  4       P                  V\        ) ^,          4      P                  4       pV\        \        V 4      ,          \        V ) ^,           4      ,          \        V 4      ,           8X  g   Q hR# )r)   r*   x0C_1N)         )rM   rM   r>   rM   )rK      rJ   i)rL   r      )rK   rP   rM   )r   r   r&   r+   r   r   r   Halfr   r   r   r   to_exprsubsr   expand)r)   r1   r*   r2   r3   r:   rI   s          r4   test_addition_initial_conditionrV   q   s   A!""2"21"5t<EA"Q$1qc*A"a%'1a!Q0A"r'BE/BE11aCA5A::"&2q5.!QA7A"a%!)QAq62AAqD51a46>AqD08Aq>AadQPQTRSVmVWX`abdeXfVfFfijFjlnEnn	
8B?	Xa^	+RU2367dQx1~=M6MPXYZ\]P^6^`bde`e5efijlmimpqrspsisvwv|v|i|  A  CD  D  iDD FG  IJ  LY	ZA5A::"a%!A#b&.1a4/A1v>A"a%!)QAq62A1a4!AqD&=1QT61A5!AqD&2ad7:JRPQT:QSU9UU
Q$1a4-!AqD&
 1
$b!e	+,/0Av1a4/?"Q$/FA.MNQRTUQUq!tGRRQQ !1o	7A 5A::"a%!)Q95A"q&!Q,AAqD519q=QTAXrM9aR!VRUNJ	
QA!^	-A5A::#a&!A!A#!$A1a4!81qs
B6!Q$(BE9IIQPQTTUVWTWZY[]^Y^L^^	1s1vz2A;Q!4	6A5A::
%.C$q'"A$q!tAv,'A	
sQBqD)002A$q'	$rAv,&a0000r6   c            
         \        R 4      p \        \        P                  ! V 4      R4      w  r\	        V^,          W,          ,           ^,
          V ^ ^^.4      p\	        V^,          ^,           V ^ ^^.4      p\	        V ^,          ^V ^,          ,          ,           ^<,           ^V ,          V,          ,           V ^,          ^	V ^,          ,          ,           ^,           V^,          ,          ,           ^V ^,          ,          ^V ,          ,           V^,          ,          ,           V ^,          ^
,           V^,          ,          ,           V ^ . RO4      pW4,          V8X  g   Q h\	        V^,          V ,           V ^ ^^ .4      p\	        V^,          V ^,          ,
          V ^ . RO4      p\	        V ^,          ^%V ^,          ,          ^,          ,
          ^
V ^,          ,          ^,          ,
          ^V ^,          ,          ^	,          ,
          ^V ^,          ,          ^	,          ,
          ^V ^,          ,          ^,          ,           RV ^,          ,          ^	,          ,           ^V ,          ,           \        ^(^4      ,           ^V ^,          ,          ^V ^,          ,          ^	,          ,
          ^bV ^,          ,          ^	,          ,
          ^V ^,          ,          ^	,          ,
          ^V ^,          ,          ^	,          ,           ^V ^,          ,          ,           V \        ^(^	4      ,          ,           ^(,
          V,          ,           ^V ^,          ,          ^RV ^,          ,          ^	,          ,
          ^LV ^,          ,          ^	,          ,           ^V ^,          ,          ^,          ,           ^V ^,          ,          ^	,          ,           V \        ^P^4      ,          ,
          V^,          ,          ,           RV ^,          ,          ^V ^,          ,          ^,          ,           ^V ^,          ,          ^,          ,
          ^PV ^,          ,          ^	,          ,
          \        ^^	4      ,           V^,          ,          ,           ^V ^,          ,          ^@V ^,          ,          ^	,          ,
          ^V ^,          ,          ^	,          ,
          ^V ^,          ,          ,           V \        ^^	4      ,          ,
          \        ^^4      ,
          V^,          ,          ,           RV ^,          ,          ^@V ^,          ,          ^	,          ,           V \        ^^4      ,          ,           V^,          ,          ,           V ^,          ^@V ^,          ,          ^,          ,
          ^V ^,          ,          ^,          ,
          \        ^^	4      ,           V^,          ,          ,           V ^ . R	O4      pW4,          V8X  g   Q h\	        V^,
          V ^ ^.4      p\	        V^,          ^,           V ^ ^ ^.4      p\	        ^^V,          ,
          V^,          ,           V ^ ^ ^.4      pW4,          V8X  g   Q h\	        W^,          ,          ^,           ^V,          ,           V ^ ^ ^.4      p\	        V ^,
          RV ,          ^,           V,          ,           W^,          ,          ,           V ^ ^ ^.4      pW4,          V8X  g   Q h\	        V^,          ^,
          V ^ ^^.4      p\	        V^,          ^,           V ^ . R
O4      p\	        ^V,          ^V^,          ,          ,           ^V^,          ,          ,           ^V^,          ,          ,
          V^,          ,           V ^ . RO4      pW4,          V8X  g   Q h\        \        V 4      4      p\        ^V ,          ^R7      p\	        V ^V,          ,           W^,          ,          ,           V ^\        ^4      \        ^4      ) \        ^4      ,           .4      pW4,          V8X  g   Q h\        \        V 4      4      p\        \        V ^,          V ,
          4      4      pW4,          P                  4       pV\        V ,          \        V ) ^,           4      ,          8X  g   Q hR# )r)   r*   i  rG   N)rL   rK   rM   rL   )rL   rL   rL   r>   r8   )rL   rL   rL   rN   ii)rP   rM   rP   )rP   rO         rY   rM   )r   r   r&   r+   r   r   r   r   r   r   rS   r   r?   s         r4   %test_multiplication_initial_conditionrZ      s   A!""2"21"5t<EA"a%!$,*Aq1a&9A"a%!)QAq62A1a4"QT'>B.!A#b&8AqD1QT6MB<NPRTUPU;UU	
1a4"Q$A!"ABE 1234a	GA5A::"a%!)QAq62A"a%!Q$,1i8A1a4"QT'"*,r!Q$wrz9C1HQJFQPQTRSSAqD!Q$hqj!#$Q3')1"a9=>q!tVc!Q$hqj=PSUVWYZVZSZ[\S\=\_abcefbf_fgh_h=h	!Q$q>ad7>xA.>/13>457=88;<QT6Bq!tGAI;MPRSTVWSWPWXYPY;Y\]^_ab^b\bcd\d;dAqD
<xA&<'(*A;.. 24AqD3q!t8B;1F1a4PQ1QSUVWYZVZSZ[\S\1\_ghkmn_o1oqsuvqv0vw 
1a4"QT'!)	bAgai	'!AqD&	01Xb!_3D	DxPRTU	VXZ\]X]]	^ bddeghdhahkmnoqrnrkrstktathq!nb !ea$	$
 ()!tbAgbj'81QT6!8'ChrSTo'UWY[\W\&\]
 _`abd{	}A 5A::"q&!Q,A"a%!)QAq62A!QrT'BE/1a!Q8A5A::!E'A+",aAq6:A1q5RTAXrM1A!eG;QAq6JA5A::"a%!)QAq62A"a%!)Q95A!B$2q5.1RU72Qr1uW<r1uDaLabA5A::#a&!A!A#!$A!ad(Q1uW,aSVc!fWs1v=M4NOA5A::$q'"A$q!tAv,'A	
A!D!aL    r6   c                  (   \        R 4      p \        \        P                  ! V 4      R4      w  r\	        V^,
          V 4      P                  V ^,          V ,           4      p\	        RV ,          ^,
          V,           V 4      pW48X  g   Q h\	        V^,          ^,           V 4      P                  V ^,          V ^,          ,           ^,           4      p\	        ^}V ^,          ,          ^V ^	,          ,          ,           ^<V ^,          ,          ,           ^V ^,          ,          ,           RV ^,          ,          ^,
          V,          ,           ^V ^,          ,          ^V ,          ,           V^,          ,          ,           V 4      pW48X  g   Q h\	        V^,          V ,          V ,           V 4      P                  ^V ^,          ,          V ^,          ,           ^,           4      p\	        ^V ^	,          ,          RV ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,           ^pV ^,          ,          ,           ^$V ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ,          ,
          ^,
          V,          ,           ^V ^,          ,          ^V ^,          ,          ,           V ^,          ,           ^V ^,          ,          ,           V ,           V^,          ,          ,           V 4      pW48X  g   Q h\	        V^,          ^,           V 4      P                  ^V ^,          ,
          4      p\	        ^V ^,          ,          V,
          W^,          ,          ,           V 4      pW48X  g   Q h\	        V^,          ^,           V 4      P                  V ^V ^,          ^,           ,          ,
          4      p\	        V ^,          ^V ^
,          ,          ,           ^V ^	,          ,          ,           ^V ^,          ,          ,           ^0V ^,          ,          ,           ^DV ^,          ,          ,           ^HV ^,          ,          ,           ^oV ^,          ,          ,           ^pV ^,          ,          ,           ^6V ^,          ,          ,           ^V ,          ,           ^,           ^V ^,          ,          ^ V ^,          ,          ,           ^V ^,          ,          ,           ^,
          V,          ,           V ^,          ^V ^
,          ,          ,           ^V ^	,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ,          ,           ^,           V^,          ,          ,           V 4      pW48X  g   Q hR# )r)   r*   iD  Nr>   i)r   r   r%   r+   r   composition)r)   r1   r*   r2   r:   s        r4   "test_HolonomicFunction_compositionr]      s}   A!""2"21"5t<EA"Q$"..q!tAv6A2a4!8r/1-A6M6"a%'1%11!Q$q!t)A+>A3q"u9s1a4x/"QT'9AadFBs1a4xRS|UWFWW	
1a4!A#r1u	!A6M6"a%'!)Q'33AadF1a4KMBA3q!t8c!Q$h.QT9C1HDs1a4xO
1a4AqD&1Wr!Q$w.1a47!A#=A2EFIJ1a4RSTUWXTXRX	1JAvJJq5I!!"#	%A 6M6"a%'1%11!AqD&9A1QT6R-!E'115A6M6"a%'1%11!aAl2BCA1b51QU7?R1W4r!Q$w>AqDH2aQRd7R
1a4ad( AX&(*1a4024Q$79:;?A!Q$wAqD?P
1a4@@?b51QU7?Qq!tV3bAg=1a4G"QPQT'QTVWXZ[W[T[[
1a4QT'adF#%&qS)+,-.0!e44 67	8A 6M6r6   c                     \        R 4      p \        \        P                  ! V 4      R4      w  r\	        ^^.\        ^^4      .V ^,          ^,          4      p\        ^V ,          ^V ^,          ,          ^,
          V,          ,           V ^,          ^V ,          ,
          V^,          ,          ,           V ^^\        ^4      ,          \        ,          ^	,          R\        ^4      ,          \        ,          ^,          \        ^^4      ,           .4      p\        V4      pWT8X  g   Q h\        \	        ^.\        ^^4      .V ^,          ^,          4      4      p\        V ) V ^,          ) ^,          ^,           V,          ,           W^,          ,          ,           V 4      pRp\        VP                  4      V8X  g   Q hVP                  VP                  8X  g   Q hR# )r)   r*   z?[sqrt(pi)*exp(1/4)*erf(1/2), -sqrt(pi)*exp(1/4)*erf(1/2)/2 + 1]Nr8   )r   r   r&   r+   r    r   r   r   r   r   r"   y0annihilator)r)   r1   r*   r2   r3   r:   r_   s          r4   test_from_hyperra      su   A!""2"21"5t<EAq!fx1~&1Q/A1Q31QT6A:r/1QTAaCZQ4FF1qQUVWQXyY[|\]~_abfghbi_ijl_lmo_orz{|~  sA  `A  OB  	CA1A6M65!x1~.1Q78A1"Aa!R//!E'91=A	JB:==AMM)))r6   c                     \        R 4      p \        \        P                  ! V 4      R4      w  r\	        \        . \        ^^4      .3\        P                  .\        P                  ^.3V 4      4      p\        V ^,          \        ^^4      ,
          V ^,          ) V ^,          ,           V,          ,           V ^,          V^,          ,          ,           V ^,          V^,          ,          ,           V ^^\        \        4      ,          ^^\        \        4      ,          ,          R^\        \        4      ,          ,          .4      pW48X  g   Q h\	        \        . . 3^ .. 3V 4      4      p\        ^V,           V ^ ^.4      pW48X  g   Q h\	        \        ^.. 3\        P                  .^ .3V 4      4      p\        V \        P                  ,           V,          W^,          ,          ,           V ^\        \        4      \        ^4      ,          \        R4      .4      pW48X  g   Q h\	        \        ^ .^.3^ .. 3^V ^,          ,          4      4      p\        ^V ^,          ,          ^,
          V,          V ^,          V^,          ,          ,           V ^\        \        R^4      4      ) ^,           \        \        R^4      4      ) .4      pW48X  g   Q hR# )r)   r*   NrJ   )r   r   r&   r+   r   r!   r   r   rR   r   r   r   r   r   r0   s        r4   test_from_meijergrc      s$   A!""2"21"5t<EAWb8Aq>"23qvvh5LaPQA!A#A.1a4%!A#+r1AAAqDQJNQRTUQUVXZ[V[Q[[]^`a	
48Q$r(
^R48_5	7A6M6Wb"XRy!45A!b&!Q,A6M6Wqc2Y!&&A3;<A1qvv:r/A!eG3QDHSVOSQSW;UVA6M6WqcA3Z1#rAadF;<A1QT6A:r/AqDQJ61HRQRO@T?TWX?X[^_ghjlm_n[oZo>pqA6M6r6   c                     \        R 4      p \        \        P                  ! V 4      R4      w  r\        RRR7      p\	        \        P                  ! V4      R4      w  rE\        V ^,          V^,          ,          V ,           V,           V 4      P                  4       p\        ^V^,           V^,          ,          ,           V^,          ^V^,          ,          ,           ^V^,          ,          ,           ^V,          ,           V^,          ,          ,           4      ^ ^3.pWg8X  g   Q h\        V ^,          V^,          ,          V ^,          ,           V^,          ,           V 4      P                  4       p\        ^V^,          ^V^,          ,          ,           ^HV^,          ,          ,           ^V,          ,           ^,           V^,          ,          ,           4      ^ ^ 3.pWg8X  g   Q h\        V ^,          V^,          ,          ^,           V^,          ,           V 4      P                  4       p\        ^V^,          ^V^,          ,          ,
          V^,          ,
          ^V,          ,           V,          ,           V^,          ^V,          ,           ^,           V^,          ,          ,           4      ^ ^ 3.pWg8X  g   Q h\        ^V ^,          ,          V^,          ,          ^V ,          V,          ,           W^,          ,          ,           V 4      P                  4       p\        ^V,          ^V^,          ,          ^V^,          ,          ,
          ^V^,          ,          ,
          ^V,          ,           V,          ,           V^,          ^V^,          ,          ,           ^V,          ,           V^,          ,          ,           4      ^ ^3.pWg8X  g   Q hR# )r)   r*   rD   TintegerSnN)r   r   r%   r+   r	   r   to_sequencer
   )r)   r1   r*   rD   _rg   r2   r3   s           r4   test_to_Sequencerj      s   A!""2"21"5t<EAT"A 0 0 3T:EA!Q$r1u*q.2-q1==?A
AQA-A!Q$AqD1H1Q31NPRTUPU0UU
VXY[\	]^A6M6!Q$r1u*q!t+b!e3Q7CCEA
AA1a4"QT'!9CE!AC!GQ NN
OQRTU	VWA6M6!Q$r1u*q.2q50!4@@BA
AA!Q$A!5!!;R ??1a4!A#:PQ>SUWXSXBXX
Y[\^_	`aA6M6!AqD&Q,1R/!E'91=IIKA
AaC1QT6AadF?Qq!tV#;ac#A2"EEAPQRSUVRVPVYZ[\Y\I\^`bc^cHcc
dfgij	klA6M6r6   c                  <   \        R 4      p \        \        P                  ! V 4      R4      w  r\        RRR7      p\	        \        P                  ! V4      R4      w  rE\        V^,
          V ^ ^.4      P                  4       p\        RV^,           V,          ,           ^4      ^ 3.pWg8X  g   Q h\        V^,          ^,           V ^ ^ ^.4      P                  4       p\        ^V^,          ^V,          ,           ^,           V^,          ,          ,           ^ ^.4      ^ 3.pWg8X  g   Q h\        V^,          ^,           V ^,          V,          ,           V ^ ^^.4      P                  4       p\        W5^,          ,           V^,          ^V,          ,           ^,           V^,          ,          ,           ^^R\        R^4      \        ^^4      .4      ^3.pWg8X  g   Q h\        V ^,          V^,          ,          ^,           V,           V 4      P                  4       p\        ^V^,           V,          ,           V^,          ^V^,          ,          ,
          ^V,          ,           V^,          ,          ,           4      ^ ^3.pWg8X  g   Q h\        R4      w  rr\        \        ^V ^,          ,           4      4      p\        V^,          V^,          ^V,          ,           V^,          ,          ,           ^ ^ V
.4      ^ ^3.pVP                  4       V8X  g   Q hVP                  4       p\        V^,           V^,           V^,          ,          ,           V^ .4      ^^ 3.pVP                  4       V8X  g   Q h\        \        V 4      V ,           4      P                  4       p\        ^V^,          ,          ^V,          ,
          V^,          ^V^,          ,          ,           V,
          ^,
          V^,          ,          ,           ^ ^^\        \        4      ,          ,           ^ V.4      ^ ^3.pWg8X  g   Q hR# )	r)   r*   rD   Tre   rg   C_0, C_1, C_2, C_3NrJ   )r   r   r&   r+   r	   r   rh   r
   r   r   r   r9   r   r   r   )r)   r1   r*   rD   ri   rg   r2   r3   C_0rI   C_2C_3s               r4   #test_to_Sequence_Initial_Coniditonsrp     sT   A!""2"21"5t<EAT"A 0 0 3T:EA"q&!Q,88:A
B!a%OQ
/	34A6M6"a%!)QAq62>>@A
AA!aQ 66A
?	CDA6M6"a%!)ad2g-q!aV<HHJA
AAIA!b"a%(??!QHUWYZO]efgik]lAm
npq	rsA6M6!Q$r1u*q.2-q1==?A
AQ
NadQq!tVmac.A2q5-HH
I1a	PQA6M6 !56Cc#a1f+&A
AqDAqD1Q3JA#551c{
CQ	JKA==?a	A
QUq1ub!em3c1X
>1	EFA==?a#a&1*%113A
Qq!tVac\adQq!tVma.?!.CRU-JJQPQTUVZ[]V^T^P^`acfLg
hjkmn	opA6M6r6   c                  .   \        R 4      p \        \        P                  ! V 4      R4      w  r\	        V^,          ^V ,          V,          ,           V ^ ^ ^.4      P                  ^
R7      pW ^,          ^,          ,
          V ^,          ^
,          ,           V ^,          ^*,          ,
          V ^	,          ^,          ,           \        V ^
,          4      ,           pW48X  g   Q h\	        V^,
          V 4      P                  V ^,          ^ ^.4      p\	        V^,          ^,           V ^ ^^ .4      pW4,          P                  ^
R7      p^V ^,          ^,          ,           V ^,          ^,          ,           ^V ^,          ,          R,          ,
          ^V ^,          ,          R,          ,
          \        V ^
,          4      ,           pWV8X  g   Q h\	        ^V ,           V^,          ,          V,           V ^ ^ ^.4      pW7,          V,           P                  ^
R7      p^V ,           V ^,          ,
          ^V ^,          ,          ^,          ,           ^V ^,          ,          ^,          ,
          ^V ^,          ,          ^,          ,           RV ^,          ,          R,          ,
          ^GV ^,          ,          ^i,          ,           RV ^,          ,          R,          ,
          RV ^	,          ,          R	,          ,           \        V ^
,          4      ,           pWV8X  g   Q h\	        ^^V ,          ,           ^V ^,          ,          ,
          ^
V ,          ^V ^,          ,          ,
          ^V ^,          ,          ,
          V,          ,
          ^^V ^,          ,          ,
          ^V ^,          ,          ,           V^,          ,          ,           V ^ ^ ^.4      P                  ^R7      pW ^,          ^,          ,           ^V ^,          ,          ^,          ,
          V ^,          ^,          ,           ^V ^,          ,          R
,          ,
          \        V ^,          4      ,           pW48X  g   Q h\	        ^^V ,          ,           ^V ^,          ,          ,
          ^
V ,          ^V ^,          ,          ,
          ^V ^,          ,          ,
          V,          ,
          ^^V ^,          ,          ,
          ^V ^,          ,          ,           V^,          ,          ,           V ^ ^^ .4      P                  ^R7      p^^V ^,          ,          ^,          ,
          V ^,          ^,          ,
          ^V ^,          ,          ^ ,          ,
          ^V ^,          ,          ^(,          ,
          ^V ^,          ,          R,          ,
          \        V ^,          4      ,           pW48X  g   Q h\        \        V 4      V ,           4      P                  ^
R7      p\        R4      p\        V 4      V ,           P                  ^
R7      pVP                  VR^\        \        4      ,          ,          4      V8X  g   Q h\        \        V ^,          V ,           4      4      P                  ^
R7      \        V ^,          V ,           4      P                  ^
R7      8X  g   Q h\        ^V ,          ^V ^,          ,          ,
          \        ^^4      ,          4      P                  4       ^V ,          ^V ^,          ,          ,
          \        ^^4      ,          P                  4       8X  g   Q h\        \        V ^,          V ,
          4      4      P                  4       \        V ^,          V ,
          4      P                  4       8X  g   Q h\        \        V 4      ^,          V ^,          ,          R. RO/R7      P                  ^
R7      \        V 4      ^,          V ^,          ,          P                  ^
R7      8X  g   Q h\        \        V 4      ^,          V ^,          ,          ^R7      P                  ^
R7      P                  4       \        V 4      ^,          V ^,          ,          P                  ^
^R7      P                  4       8X  g   Q h\        \        V ^,
          4      ^,          V ^,
          ^,          ,          ^R. RO/R7      P                  ^
R7      \        V ^,
          4      ^,          V ^,
          ^,          ,          P                  ^^
R7      8X  g   Q hR# )r)   r*   rD   i  i  i  iN  i  i{  iH  i  i  ro   r_   rG   )rD   rH   )rH   r_   )rH   rD   Nr>   rP   r   rJ   )r   r   r%   r+   r   seriesr#   r\   r   r   rT   r   r   r   r   together)	r)   r1   r*   r2   r3   r:   r;   tro   s	            r4   test_seriesrx   $  s   A!""2"21"5t<EA"a%!A#b&.!QA7>>>DA	qDF
QT"Wq!tBw&Ac1AaeH<A6M6"q&!$00Aq1#>A"a%!)QAq62A	
A	AqDF
QT"Wr!Q$ws{*SAXd]:Qq"uXEA6M61q5"a%-",aQF;A	
R A	A1qAvax"QT'"*,r!Q$wrz9C1HSLH1WS[	A:e#	$&)!Q$hsl	356q"uX	>A6M61QqS51a4<BqD1a4K!Q$,>+BB	
1QT6!AqD&"a% !"A1v	//5vv{ 	qDF
Qq!tVBYAb(2ad73;61a4@A6M61QqS51a4<BqD1a4K!Q$,>+BB	
1QT6!AqD&"a% !"A1v	//5vv{ 	AadF1Hq!tAv!Q$r	)AadF2I51a4Ca1gMA6M6#a&1*%,,r,2A
%.C	Q!b!A66#r1T"X:'1,,,T!Q$(^,33b39T!Q$(^=R=RUW=R=XXXXacAadFlXa^;<CCE1Q3QRSTVWSWQW<ZbcdfgZhJhIpIpIrrrrd1a46l+224ad1f8M8M8OOOOSVQYq!t^Z0@AHH2HNSVWXSY[\S\]^`a]aSaRiRilnRiRooooSVQYq!t^299B9?HHJsSTvWXyYZ\]Y]~NeNehjopNeNqNzNzN||||S1Xq[!A#1aR<LMTTWYTZ!HaK1q ((A(45 5 5r6   c                  n	   \        R 4      p \        \        P                  ! V 4      R4      w  r\	        ^V ,           V^,          ,          V,           V ^ ^ ^.4      p. ROpRp\        VP                  VRR7      R,          4      V8X  g   Q hRR\        ,          ,           .p\        ^	4       F4  pVP                  VR,          R,           R\        ,          ,           4       K6  	  \        ^
4       F4  pVP                  VR,          R,           R\        ,          ,
          4       K6  	  Rp\        VP                  VRR7      R,          4      V8X  g   Q h\	        V^,          ^,           V ^ ^ ^.4      pRp\        VP                  VRR7      R,          4      V8X  g   Q hR.p\        ^4       F"  pVP                  VR,          R,           4       K$  	  VP                  \        ^,          4       Rp\        VP                  VRR7      R,          4      V8X  g   Q hR\        ,          .p\        ^	4       F-  pVP                  VR,          R\        ,          ,           4       K/  	  \        ^4       F"  pVP                  VR,          R,           4       K$  	  VP                  \        ^,          \        ,           4       \        ^
4       F-  pVP                  VR,          R\        ,          ,
          4       K/  	  R	p\        VP                  VRR7      R,          4      V8X  g   Q h\	        V^,          ^,           V ^ ^^ .4      pR
.p\        ^=4       F"  pVP                  VR,          R
,           4       K$  	  VP                  \        4       Rp\        VP                  VRR7      R,          4      V8X  g   Q hR\        ,          .p\        ^	4       F-  pVP                  VR,          R\        ,          ,           4       K/  	  \        ^4       F"  pVP                  VR,          R,           4       K$  	  \        ^
4       F-  pVP                  VR,          R\        ,          ,
          4       K/  	  \	        V^,          ^,           V ^ ^^.4      P                  VRR7      pRp\        VR,          4      V8X  g   Q hR# )r)   r*   皙?z0.699525841805253Euler)methodz'1.07530466271334 - 0.0251200594793912*Iz*0.905546532085401 - 6.93889390390723e-18*Iz1.08016557252834z*0.976882381836257 - 1.65557671738537e-16*I皙?z-1.08140824719196z*0.501421652861245 - 3.88578058618805e-16*IN
rz   g?g333333?g?g      ?g333333?gffffff?g?g?rP   rJ   )r   r   r&   r+   r   r"   evalfr   rB   appendr   r)   r1   r*   r2   r:   r;   is          r4   test_evalf_eulerr   H  s   A!""2"21"5t<EA 	1q5"a%-",aQF;A 	9AA'*2./1444 
s1uA1X	2s3q5! 2Y	2s3q5! 
 	2A'*2./1444 	"a%!)QAq62A4A'*2./1444 
A2Y	2 HHRTNA'*2./1444 
QA1X	2s1u 2Y	2s HHRT!V2Y	2s1u  	5A'*2./1444 	"a%!)QAq62A	A2Y	2t HHRLA'*2./1444 
QA1X	2s1u 2Y	2s 2Y	2s1u  	"a%!)QAa5177'7JA4A";!r6   c                  p	   \        R 4      p \        \        P                  ! V 4      R4      w  r\	        ^V ,           V^,          ,          V,           V ^ ^ ^.4      p. ROpRp\        VP                  V4      R,          4      V8X  g   Q hRR\        ,          ,           .p\        ^	4       F4  pVP                  VR,          R,           R\        ,          ,           4       K6  	  \        ^
4       F4  pVP                  VR,          R,           R\        ,          ,
          4       K6  	  Rp\        VP                  V4      R,          P                  ^4      4      V8X  g   Q h\	        V^,          ^,           V ^ ^ ^.4      pRp\        VP                  V4      R,          4      V8X  g   Q hR.p\        ^4       F"  pVP                  VR,          R,           4       K$  	  VP                  \        ^,          4       Rp\        VP                  V4      R,          4      V8X  g   Q hR\        ,          .p\        ^	4       F-  pVP                  VR,          R\        ,          ,           4       K/  	  \        ^4       F"  pVP                  VR,          R,           4       K$  	  VP                  \        ^,          \        ,           4       \        ^
4       F-  pVP                  VR,          R\        ,          ,
          4       K/  	  Rp\        VP                  V4      R,          4      V8X  g   Q h\	        V^,          ^,           V ^ ^^ .4      pR.p\        ^=4       F"  pVP                  VR,          R,           4       K$  	  VP                  \        4       R	p\        VP                  V4      R,          4      V8X  g   Q hR\        ,          .p\        ^	4       F-  pVP                  VR,          R\        ,          ,           4       K/  	  \        ^4       F"  pVP                  VR,          R,           4       K$  	  \        ^
4       F-  pVP                  VR,          R\        ,          ,
          4       K/  	  \	        V^,          ^,           V ^ ^^.4      P                  V4      pR
p\        VR,          4      V8X  g   Q hR# )r)   r*   rz   z0.693146363174626z1.098616 + 1.36083e-7*Iz)0.90929463522785 + 1.52655665885959e-16*Iz0.999999895088917z)1.00000003415141 + 6.11940487991086e-16*Ir}   z-0.999999993238714z*0.493152791638442 - 1.41553435639707e-15*INr~   rJ   )r   r   r&   r+   r   r"   r   r   rB   r   rD   r   r   s          r4   test_evalf_rk4r     s   A!""2"21"5t<EA 	1q5"a%-",aQF;A 	9AA
21$$$ 
s1uA1X	2s3q5! 2Y	2s3q5! 
 	"A
2  #$))) 	"a%!)QAq62A3A
21$$$ 
A2Y	2 HHRTNA
21$$$ 
QA1X	2s1u 2Y	2s HHRT!V2Y	2s1u  	4A
21$$$ 	"a%!)QAq62A	A2Y	2t HHRLA
21$$$ 
QA1X	2s1u 2Y	2s 2Y	2s1u  	"a%!)QAa5177:A4A";!r6   c                     \        R 4      p \        \        P                  ! V 4      R4      w  r\	        \        V 4      V ,          ^,          4      p\        ^V ,          ^V ^,          ,          ^,           V,          ,           ^V ,          V^,          ,          ,           V ^,          V^,          ,          ,           V ^ ^^ \        R	^4      .4      pW48X  g   Q h\	        ^^V ^,          ,           ^,          ,          4      p\        ^V ,          V ^,          ^,           V,          ,           V ^ ^.4      pW48X  g   Q h\	        \        V 4      \        V 4      ,          V \        ^V ,           4      ,          ,           4      p\        ^V ^,          ,          ^
V ^,          ,          ,           ^V ,          ,           ^,           R	V ^,          ,          ^
V ^,          ,          ,
          ^V ^,          ,          ,
          ^V ,          ,
          V,          ,           ^V ^,          ,          ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,           ^
V ,          ,           ^,
          V^,          ,          ,           R	V ^,          ,          ^V ^,          ,          ,
          ^V ^,          ,          ,
          ^V ^,          ,          ,           V ,
          ^,           V^,          ,          ,           V ^,          ^V ^,          ,          ,           V ^,          ,
          ^V ^,          ,          ^,          ,
          V ,           \        ^^4      ,           V^,          ,          ,           V ^ . R
O4      pW48X  g   Q h\	        V \        V 4      ,          \        V 4      ,           ^,           4      p\        V ) ^,
          V,          V ^,           V^,          ,          ,           V ) ^,
          V^,          ,          ,           V ^,           V^,          ,          ,           V ^ . RO4      pW48X  g   Q hV \        V 4      ,          \        V 4      ,           ^,           P                  ^
R7      VP                  ^
R7      8X  g   Q h\	        \        ^V ,           4      ^,          ^,           4      p\        V^V ,          ^,           V^,          ,          ,           V ^,          ^V ,          ,           ^,           V^,          ,          ,           V ^ . RO4      pW48X  g   Q h\	        \        V 4      ^,          V ,           4      p\        ^V ^,          ,          ^V ^,          ,          ,
          ^,           V^,          ,          ^V ^,          ,          V ^,          ,
          V^,          ,          ,           V ^,          \        ^^4      ,           V^,          ,          ,           V ^ ^ ^^\        ,          ^ .4      pW48X  g   Q h\	        \        V 4      V ,          4      p\        V ^,          ) ^,           ^V ,          V,          ,
          V ^,          V^,          ,          ,           V ^ ^ ^.4      pW48X  g   Q h\	        \        ^V 4      4      p\        V ^,          ^,
          W,          ,           V ^,          V^,          ,          ,           V ^ ^ ^ .4      pW48X  g   Q h\	        \        ^ V 4      \        V 4      ,           4      p\        V ^,          ) V ^,          ,
          \         P"                  ,           V ^,          V ^,          ,
          \        ^^4      ,
          V,          ,           V ^,          ) V ^,          ,           ^,           V^,          ,          ,           V ^,          V ^,          ,           V^,          ,          ,           V ^ ^^\         P"                  .4      pW48X  g   Q h\	        \        V 4      ^,          V ,          4      p\        ^^V ,          V,          ,           ^V^,          ,          ,           W^,          ,          ,           V ^ . RO4      pW48X  g   Q h\	        \        V 4      ^,          V ,          ^R7      p\        ^^V ,          V,          ,           ^V^,          ,          ,           W^,          ,          ,           V ^\        ^4      ^,          ^,          \        ^4      \        ^4      ,          \        ^4      ^,          ^,          ,
          R\        ^4      ^,          ,          ^,          \        ^4      ^,          ,           \        ^4      \        ^4      ,          ,
          .4      pW48X  g   Q h\	        \        V 4      ^,          \%        ^V ,          4      ^,          ,
          \%        ^4      ^,          ,           4      p\        ^V,          ^V ,          V^,          ,          ,           ^V^,          ,          ,           W^,          ,          ,           V ^ \        ^4      ) ^,          \&        ^,          ,
          \%        ^4      ^,          ,           ^ ^^ .4      pW48X  g   Q hVP)                  4       p\        V 4      ^,          \%        ^V ,          4      ^,          ,
          \%        ^4      ^,          ,           pW48X  g   Q h\	        V \         P"                  ,          ^R7      p\        W,          \         P"                  ,
          V ^^.4      pW48X  g   Q h\	        \+        ^V ^,          ,           4      4      p\        V ) V ^,          ^,           V,          ,           V ^ ^.4      pW48X  g   Q h\	        \+        V 4      \+        ^V ,          4      ,           4      P)                  4       \+        V 4      \+        ^V ,          4      ,           ,
          P-                  4       ^ 8X  g   Q h\	        ^V ,          ^\+        V 4      ,          ,           4      P)                  4       ^V ,          ^\+        V 4      ,          ,           8X  g   Q h\	        V ^,          V ^,          ,           ^V ^,          ,          ,           ^V ,          ,           ^,           V ^,          ,          ^R7      p\        R	V ^,          ,          V ^,          ,
          ^V ,          ,           ^,           V ^,          V ^,          ,           ^V ^,          ,          ,           ^V ^,          ,          ,           ^V ,          ,           V,          ,           V ^ R	. RO/4      pW48X  g   Q h\	        ^V ^,
          ^,          ,          ^^R7      p\        ^V ^,
          V,          ,           V ^R	. RO/4      pW48X  g   Q h\        R4      p\	        \+        WP,          4      V R7      pVP)                  4       \+        V4      \+        V 4      ,          8X  g   Q hR# )r)   r*   rr   rG   lenics)r   rH   rC   r)   Nr>   )r   rP   rK   rJ   )rM   rP   rP   rL   )rP   r   rM   )r   rP   r   rN   )rM   rL   rO   rP   r   r   )r   r   r&   r+   r   r   r   r   r   r   r   ru   r   r   r   r   r   rR   r   r   rS   r   simplify)r)   r1   r*   r2   r3   rC   s         r4   test_expr_to_holonomicr     s	   A!""2"21"5t<EA3q6!8a-(A!A#1a4!R/!A#b!e);ad2q5jH!Q	
AxA	!A6M6!Qq!tVaK-(A!A#A2-q!aS9A6M6#a&Q-#ac(
23A1QT6Bq!tG+bd2R7Bq!tGbAg<MPRSTVWSWPW<W
Q$=< 1fqAvo!Q$.1a47"Q$>BBEIJ	AqD1QT6	AadF	"Qq!tV	+a	/!	3RU:;>?dQq!tVmaQRd>R	!Q$q??1~?&')1u>-- /0M	CA 6M6!CF(3q6/!+,AA262+QA5!aQF!a%QSUVQVVXY	<	A6M6c!fHSVOA%%%+qxx"x~===#a!e*a-!+,A"!aQ.!Q$1*q."a%1GGAyYA6M6#a&!)a-(A1QT6AadF?Q.A51a4!A#r1u8LL	
Ax1~	r1u$%&'Q1R4O	=A6M6$q'!)$AAqD519qs2v-1RU
:Aq1a&IA6M6'!Q-(A1a4!8qt+ad2q5j8!QAGA6M6'!Q-#a&01AAqD51Q3;/AqD1Q3J!Q4OQS3SSXY[\X\W\_`ab_bWbefWfhjlmhmVmm	
A!RUq1a.	2A6M6#a&!)A+&A!ac"f*qQw.q58!Q	JA6M6#a&!)A+!,A11b.Ar1u94q5y@!QQQRSTAs1vA	!#RA	\!^c!fai%?#a&Q-%OIQ 	RA6M6#a&(R!WQY.Aq89A!B$1RU*Qr1uW4qQw>1
a&Z\	!BqE!G	+Q15	7A6M6			AAq2ac719r!uQw&A6M6!QVV)*A!$-As3A6M6$q1a4x.)AA2!Q$(B.1qc:A6M6d1gQqS	12::<	a4!9	'xzQ/ / /QqS47]+3351QtAwYFFF1a419Qq!tV+AaC/11a47BA2ad7QT>AaC/!31q!ta1f8LqQRTUQUv8U	!98 AI	(A6M6!QqS1H*Q15A1Q
*Aq2y/BA6M6A$qs)q)A99;$q'$q'/)))r6   c                     \        R 4      p \        \        P                  ! V 4      R4      w  r\	        V^,
          V ^ ^.4      P                  4       p^\        . . ^V ,          4      ,          pW48X  g   Q h\        \	        ^V ,           V,          ^,
          V ^ ^.4      P                  4       4      P                  4       p^V ^,          ,          ^V ^,          ,          ,           ^V ,          ,           ^,           pW48X  g   Q h\	        ^V ,           V^,          ,          V,           V ^ ^ ^.4      P                  4       pV ^,          ) \        RRV ) 4      ,          ^,          V ,           pW48X  g   Q h\	        ^V ,          V,          V^,          ,           V ^ ^ ^\        \        4      ,          .4      P                  4       p^V ,          \        \        P                  3\        ^^4      3V ^,          ) 4      ,          \        \        4      ,          pW48X  g   Q h\        \	        ^V ,          V,          V^,          ,           V ^ ^R\        \        4      ,          .4      P                  4       4      p\        V 4      pVP                  \        4      V8X  g   Q h\        \	        V ^,          ^,
          W,          ,           V ^,          V^,          ,          ,           V ^ ^ \        P                  .4      P                  4       4      p\!        ^V 4      pW48X  g   Q h\        \	        W^,          ,          V,           V ,           V ^ ^^ .4      P                  4       4      p\!        ^ V 4      pW48X  g   Q hR# )r)   r*   N)rM   rM   rP   )rL   rM   r>   )r   r   r&   r+   r   to_hyperr    r$   rU   r   r   r   rR   r   r   rewriter   r0   s        r4   test_to_hyperr   +  s   A!""2"21"5t<EA"q&!Q,557A	E"b!A#A6M6%q1ulQ&61qcBKKMNUUWA	!Q$1a4!A#!A6M61q5"a%-",aQF;DDFA	
AeIvr**1,q0A6M6!A#b&2q5.!QAd2hJ@IIKA	!E166)hq!n.A66tBx?A6M6%ac"fr1unaQ48<LMVVXYAQA99T?a	&1qAD'81a4A:'E	1q!&&k#8:
'A1A6M6%aAglQ&61q!fENNPQA1A6M6r6   c                     \        R 4      p \        \        P                  ! V 4      R4      w  r\	        V^,
          V ^ ^.4      P                  4       p\        V 4      pW48X  g   Q h\	        V^,          ^,           V ^ ^^ .4      P                  4       p\        V 4      pW48X  g   Q h\	        V^,          ^,
          V ^ ^^ .4      P                  4       p\        V 4      pW48X  g   Q h\	        ^^V ,          ^,
          V,          ,           V ^,          V ,
          V^,          ,          ,           V ^ ^^.4      P                  4       P                  4       p^V ^,          ^V ,          ,
          ^,           ,          pW48X  g   Q h\        \        V 4      ^,          V ,          4      P                  V ^ V 34      P                  4       p\        V 4      ^,          V ,          P                  V ^ V 34      pW48X  g   Q h\        R4      w  rVrx\        \        ^V ^,          ,           4      4      P                  4       pV\        V ^,          ^,           4      ,          pW48X  g   Q h\        \        ^V ^,          ,           4      4      P                  4       P                  4       pWP,          V ^,          ^,           ,          pW48X  g   Q h\        \        V 4      V ,           4      P                  4       p^V,          V ,          ^\!        \"        4      ,          V,          \        V 4      ,          ^,          ,
          V ,           ^V ,          \!        \"        4      ,          ,           pW48X  g   Q h\        \!        V 4      ^R7      P                  4       pV\!        V 4      8X  g   Q h\        \!        V 4      4      P                  4       \!        V 4      8X  g   Q h\        \!        ^V ^,          ,           4      4      P                  4       pV\!        ^V ^,          ,           4      8X  g   Q h\        ^V ^,          ,          ^,           \%        ^^4      ,          4      P                  4       pV^V ^,          ,          ^,           \%        ^^4      ,          8X  g   Q h\        \!        V ^,          ) ^V ,          ,           4      4      P                  4       pV\!        V 4      \!        V ) ^,           4      ,          8X  g   Q h\        R
V ^,          ,          ^V ,          ,           \%        ^^4      ,          4      P                  4       pV \%        ^^4      ,          R
V ^,          ,          ^,           \%        ^^4      ,          ,          pW48X  g   Q h\'        \)        R\*        P,                  3V 4      4      p\/        \)        R\*        P,                  3V 4      4      p	\1        R4      p
\1        R4      pVP                  4       P3                  V^V
^ /4      V	,
          P5                  4       ^ 8X  g   Q h^ ^.\*        P,                  ^ ./Vn        VP                  4       V	8X  g   Q h\        V ^,          4      P                  4       V ^,          8X  g   Q h\        ^V ^,          ,          ^V ^,          ,          ,
          4      P                  4       P                  4       ^V ^,          ,          ^V ^,          ,          ,
          8X  g   Q h\        R4      p\        RV ,          4      \        W,          V 4      ,          P                  4       pRV,          V ^,          ,          pW48X  g   Q h\        RV ,          4      \        W,          V 4      ,           P                  4       pWR,           ,          pW48X  g   Q h\        RV ,          4      \        V 4      ,           P                  4       pVRV ,          8X  g   Q hR	# )r)   r*   rl   rG   D_0rm   rC   gffffff?g333333@Nr>   r>   rN   )r   r   r%   r+   r   rS   r   r   r   rU   r   r   	integrater   r9   r   r   r   r   r   r    r   rR   r$   r   rT   r   r_   )r)   r1   r*   r2   r3   rm   rI   rn   ro   r;   r   rC   s               r4   test_to_exprr   E  s   A!""2"21"5t<EA"q&!Q,446AAA6M6"a%!)QAq62::<AAA6M6"a%!)QAq62::<AQA6M6!qsQwl*	
A2q5QA	((/	&&( 	1a4!A#:>A6M6#a&!)A+&00!Q;CCEA	Q1Aq	*A6M6 !56Cc#a1f+&..0AC1qMA6M6#a1f+&++-557Aq!taxA6M6#a&1*%--/A	#a!DH*S.Q'))A-!DH<A6M6$q'a(002AQ<<T!W%--/47:::$q1a4x.)113AQq!tV1QT6A:A67??AA1a4!hq!n,,,,$1uQqSy/*224AQaR!V$$$$2ad71Q3;!Q78@@BA	8Aq>2ad7Q;!Q77A6M65AFF:q12AE(QVVJ23A
-C
-CIIKc!S^,q0::<AAAsAFFQC AD99;!QT"**,1444Qq!tVAadF]+335<<>	!Q$qAv  A	3q5	!"3AC";	;DDFAAad
A6M6	3q5	!"3AC";	;DDFA	s7A6M6	3q5	!"3A"6	6??AAA::r6   c                  `   \        R 4      p \        \        P                  ! V 4      R4      w  r\	        \        V 4      ^,          V ,          ^R7      P                  V ^^34      pRp\        V4      V8X  g   Q h\	        \        V 4      4      P                  V ^ V 34      P                  4       p^\        V 4      ,
          pW48X  g   Q h\	        \        V 4      4      P                  V ^ ^34      p^\        ^4      ,
          pW48X  g   Q h\	        \        V 4      V ,          ^R7      P                  V ^^34      pRp\        V4      V8X  g   Q h\	        \        V 4      ^,          V ,          ^R7      P                  V ^^ 34      pRp\        V4      V8X  g   Q h\	        \        V 4      V ,          4      pVP                  V 4      P                  4       \        V 4      8X  g   Q hVP                  V ^ ^34      \        ^4      8X  g   Q h\	        \        V 4      ^,          V ,          4      pVP                  4       pVP                  V 4      P                  4       VP                  V ^ V 34      8X  g   Q hVP                  V ^ ^34      VP                  V ^ ^34      8X  g   Q h\	        ^V ,          ^R7      P                  V 4      P                  4       \        V 4      8X  g   Q h\	        V ^,           ^,          \        V ) 4      ,          R	R7      P                  V 4      P                  4       pV ^,          ) ^V ^,          ,          ,
          ^V ,          ,
          ^\        V ^,           4      ,          ,           ^,
          \        V ) 4      ,          pW48X  g   Q h\	        \        V 4      ^,          V ^,          ,          R
. RO/R7      P                  V 4      P                  4       p\        ^V ,          4      ) \        V 4      ^,          V ,          ,
          pW48X  g   Q h\	        \        V ^,          V ,           4      4      P                  V 4      P                  4       pV \        ^^4      ,          ^V ^,          ,          ^V ,          ,           ^,           ,          V \        V ^,           4      ,          \        \        V 4      4      ,          ,
          ^V ,          \        V ^,           4      ,          ,          pW48X  g   Q h\	        \        V ^,          ^,           4      4      P                  V 4      P                  4       p\        V ^,          ^,           4      P                  V 4      pW4,
          P!                  4       ^ 8X  g   Q h\	        ^V ^,          ,          R
. RO/R7      p\	        ^V ^,          ,          ^R7      pW58X  g   Q h\	        \        V 4      ^,          4      pWT,          P                  V 4      P                  4       \        ^V ,          4      ) \        V 4      ^,          V ,          ,
          8X  g   Q hR# )r)   r*   rG   z0.166270406994788z0.659329913368450z-0.423690480850035rs   r   NrJ   r>   rt   r   )r   r   r%   r+   r   r   r   r"   rS   r   r   r   r   r   r   r   r   r?   s         r4   test_integrater     s   A!""2"21"5t<EA#a&!)A+!,661ayAAA7a<<#a&!++Q1I6>>@A	CF
A6M6#a&!++Q1I6A	CF
A6M6#a&(q)33Q1I>AA7a<<#a&!)A+!,661ayAAA7a<<#a&(#A;;q>!!#r!u,,,;;1ay!RU***#a&!)A+&A			A;;q>!!#q{{Aq!9'====;;1ay!Q[[!Q%;;;;QqSQ'11!4<<>#a&HHH1q51*S!W,4>>qAIIKA
Q$1a4"Q$	3q1u:	-	2CG;A6M6#a&!)AqD.b*-=>HHKSSUA	AaC3q619Q;A6M6$q!tAv,'11!4<<>A	
HQN	AadFQqSL1,	-$q1u+eDGn0L	LqQRsSWXY\]X]S^_A6M6$q!tAv,'11!4<<>A	ad1f  #AC>>q   !AqD&b^4A!AqD&+A6M6#a&!)$AC??1%%'BqsG8c!faik+AAAAr6   c                     \        R 4      w  r\        \        P                  ! V 4      R4      w  r#\	        W^,          ,          ^,           V ^ ^ ^.4      pVP                  4       P                  4       VP                  4       P                  4       P                  4       8X  g   Q h\	        V^,          ^,
          V ^ ^^ .4      pVP                  V ^4      P                  4       VP                  4       8X  g   Q h\        \        V 4      4      pVP                  4       P                  4       \        V 4      V ,          8X  g   Q hVP                  V4      ^ 8X  g   Q h\        R4      w  rVrx\        V 4      p	VP                  V 4      P                  4       V	P                  4       8X  g   Q hVP                  V ^4      P                  4       P                  V\        R^4      4      P                  4       V	P                  V ^4      P                  4       8X  g   Q hVP                  V ^4      P                  4       P                  V\        R^4      V^ /4      V	P                  V ^4      P                  4       8X  g   Q hR# )zx, yr*   rl   NrJ   )r   r   r%   r+   r   r9   rS   r   r   r   r   rT   r   cancelru   )
r)   yr1   r*   r2   rm   rI   rn   ro   r3   s
             r4   	test_diffr     s   6?DA!""2"21"5t<EA!E'A+q!aV4A668!1!1!3!<!<!>>>>"a%!)QAq62A66!Q<!QYY[000"Q% A668Q)))66!9>> !56Cc
1A66!9!&&(***66!Q<!&&sHRO<CCEPQSTI\I\I^^^^66!Q< %%sHROS!&DEPQSTI\I\I^^^^r6   c                     \        R 4      p \        R\        RV ,          4      ,          4      pVP                  4       R\        RV ,          4      ,          8X  g   Q h\	        VP                  V 4      P                  4       4      R8X  g   Q h\        \        P                  ! V 4      R4      w  r#\        RV ,          4      p\        RRV ,          V,          ,           V ^ ^R./4      pW8X  g   Q hVP                  4       RV ,          8X  g   Q h\	        VP                  V ^^34      4      \	        RV ,          P                  V ^^34      4      8X  g   Q h\        R4      w  rV\        \        W,          V,          4      V R7      pVP                  4       \        W,          V,          4      8X  g   Q hVP                  V 4      P                  4       \        W,          V,          4      ) ^,           WV,          ,          8X  g   Q h\        \        W,          V,           4      V R7      P                  V 4      P                  4       p\        V4      \        W,          V,           4      ,
          V,          pW8X  g   Q h\        R4      p\        Wp,          V 4      pVP                  4       Wp,          8X  g   Q hVP                  V 4      P                  4       Wp^,          ,          ^,          8X  g   Q h\        R	4      w  r\        V 4      \        R\        V 4      ,          4      ,           pVP                  4       P                  V^ 4      pW,
          R\        R
V ,          4      ,          ,
          ^ 8X  g   Q h\        V 4      \        R\        V 4      ,          4      ,          pVP                  4       P                  V	^ 4      pVRV ,          \        R
V ,          4      ,          ,
          ^ 8X  g   Q hR# )r)   g333333?g@z0.387096774193548*sin(3.1*x)r*   gݷ:>j ?zy, zr   rC   zD_2, C_1g      ?Ngݷ:>j )r   r   r   rS   r"   r   r   r'   r+   r   r   rT   )
r)   r2   ri   r*   r3   r   zrC   D_2rI   s
             r4   )test_extended_domain_in_expr_to_holonomicr     s   A#c#a%j.)A99;#c#a%j.(((A&&()-KKKK!""2"21"5t<EA,q.)A=\!^R,??AL>GZ[A6M699;,q.(((Q1I&'4a0J0JAqRS90U+VVVV6?DA#ac!e**A99;#ac!e*$$$;;q>!!#QSUa!#'>>>>#acAg,!,66q9AACA	Q#acAg,	!A6M6A!#q!A99;!#;;q>!!#qAvax///z"HC!0SV<<A			a A53s3q5z>!Q&&&!0SV<<A			a As1uSQZ1$$$r6   c            	      	   \        R 4      p \        \        \        V 4      4      P	                  4       4      \        V 4      8X  g   Q h\        \        \        V 4      4      P	                  4       4      \        V 4      8X  g   Q h\        \        \        V 4      4      P	                  4       4      \        V 4      8X  g   Q h\        \        \        V 4      4      P	                  4       4      P                  4       \        V 4      8X  g   Q h\        ^V ^,          ,          ^,          ^,           4      P	                  4       ^V ^,          ,          ^,          ^,           8X  g   Q h\        \        \        ^V 4      ^R7      P	                  4       4      \        ^V 4      8X  g   Q h\        \        R	^4      R
3RV 4      p\        V4      P	                  4       \        V4      8X  g   Q h\        \        P                  \        ^4      3\        ^4      3V 4      p\        \        V4      P	                  4       4      \        V4      ,
          P                  4       ^ 8X  g   Q h\        \        R\        P                   3V 4      4      p\        \        R\        P                   3V 4      4      p\#        R4      p\#        R4      p\#        R4      p\        VP	                  4       4      P%                  V^V^ /4      V,
          P                  4       ^ 8X  g   Q h^ ^.\        P                   ^ ./Vn        \        VP	                  4       4      V,
          P                  4       ^ 8X  g   Q h\        \        \        P                   V 4      RR7      pVP)                  4       V\        V 4      ,          V\        V 4      ,          ,           V\        V 4      ,          ,           \+        V 4      ,          ,
          P                  4       ^ 8X  g   Q h\        \        \        P                   V 4      \        R	^4      \+        ^4      \+        \,        4      ,          \+        ^4      \+        \,        4      ,          ./R7      pVP)                  4       \        \        P                   V 4      ,
          \        \        R	^4      V 4      ,
          P                  4       ^ 8X  g   Q hR# )r)   r   rm   rI   r   F)initcondrs   NrJ   rN    r   )r   r$   r   r   
to_meijergr   r   r   r   r   r    r   r   r   OnerU   rR   r   rT   r_   rS   r   r   )r)   r2   r;   rm   rI   r   s         r4   test_to_meijergr     sl   A(Q0;;=>#a&HHH(Q0;;=>#a&HHH(Q0;;=>#a&HHH(Q0;;=>GGISQRVSSSQq!tVAX\*5571QT6!8a<GGG(AqALLNOSZ[\^_S````xA#R+Aa=##%Q777quuadmadXq)A
10023k!nDLLNRSSSS5AFF:q12AE(QVVJ23A
-C
-C
-C',,c!S^<q@JJLPQQQQsAFFQC AD'!+5571<<<'!&&!,u=AIIK3s1v:CF
2SQZ?aHHRRTXYYYY'!&&!,(2q/DGDQSHDTVZ[\V]^bce^fVfCg1hiAIIK'!&&!,,wxA/JJTTVZ[[[[r6   c                     \        R 4      w  r\        RRR7      p\        W3,          P                  4       p\        ^4      \	        V ) V,           ^,          ) ^V^,          ,          ,          4      ,          ^\        \
        4      ,          V,          ,          p\        WAVR7      p\        VP                  V4      R4      w  rg\        V ) V^,          ,          W^,          ,          ,           ^V,          ,           V4      pWX8X  g   Q hR# )zmu xsdTpositivedomainr*   N)
r   r&   	get_fieldr   r   r   r   r   r+   r   )	mur)   r   Qeh1ri   r*   h2s	            r4   test_gaussianr     s    FOEB		%B
26
AQrcAg\M1RU7+,,aRjm<A	1	*B!!//!"4d;EA	RCAIa%/Ar691	=B8O8r6   c                     \        R RR7      w  rpW ^,
          ,          V) ^,           V^,
          ,          ,          \        W4      ,          p\        W3,          P                  4       p\	        W2VR7      p\        VP                  V4      R4      w  rg\        WV ) V,
          ^,           ,          ,           ^,
          V^,          V,
          V,          ,           V4      pWX8X  g   Q hR# za b xTr   r   r*   N)r   r   r&   r   r   r   r+   r   	rC   br)   r   r   r   ri   r*   r   s	            r4   	test_betar     s    g-GA!	E
QBFa!e$$T!Z/A
14A	1	*B!!//!"4d;EA	AA26A:.2q!taxmCQ	GB8O8r6   c                     \        R RR7      w  rpW) ,          W ^,
          ,          ,          \        V) V,          4      ,          \        V 4      ,          p\        W3,          P	                  4       p\        W2VR7      p\        VP                  V4      R4      w  rg\        V ) ^,           W!,          ,           W',          ,           V4      pWX8X  g   Q hR# r   )	r   r   r   r&   r   r   r   r+   r   r   s	            r4   
test_gammar     s    g-GA!	BE
3r!t9$U1X-A
14A	1	*B!!//!"4d;EA	QBFQSLQF2A	6B8O8r6   c                     \        R 4      w  r\        V,          P                  4       p\        VP	                  V 4      R4      w  r4\        RW,          ,           V 4      V) ,          p\        WV,          ,           V 4      pWV8X  g   Q hR# zx nr*   NrJ   r   r&   r   r   r+   r   r)   rD   r   ri   r*   r   r   s          r4   test_symbolic_powerr     sh    5>DA
1A!!//!"4d;EA	B1&=!	,	2B	AR<	+B8O8r6   c                      \        R 4      p \        \        P                  ! V 4      R4      w  r\	        RW,          ,           V 4      R,          p\	        ^W,          ,           V 4      pW48X  g   Q hR# )r)   r*   NrJ   r>   )r   r   r&   r+   r   )r)   ri   r*   r   r   s        r4   test_negative_powerr   )  sV    A!""2"21"5t<EA	B1&=!	,	2B	A!<	+B8O8r6   c                  0   \        R 4      w  r\        V,          P                  4       p\        VP	                  V 4      R4      w  r4\        RW,          ,           V 4      V^,
          ,          p\        V) ^,           W,          ,           V 4      pWV8X  g   Q hR# r   r   r   s          r4   test_expr_in_powerr   2  sq    5>DA
1A!!//!"4d;EA	B1&=!	,Q	7B	QBFqf,a	0B8O8r6   c                  R   \        R RR7      p \        \        P                  ! V 4      R4      w  r\	        V ^,          VP
                  P                  VP
                  P                  .V4      p\	        V ^,          ^V .V4      pW48X  d   Q hVP                  ^ ,          pWE8X  d   Q hR# r)   Tre   r*   N)r   r   r&   r+   r   r-   r.   
listofpoly)r)   r1   r*   dodo2p2s         r4   test_DifferentialOperatorEqPolyr   <  s    T"A!""2"21"5t<EA	q!tQVV[[!&&++>	BB
1a|Q
/Cy= 
	By=yr6   c                  t   \        R RR7      p \        \        P                  ! V 4      R4      w  r\	        V ^,          VP
                  P                  VP
                  P                  .V4      p\	        VP
                  P                  .V4      p\        ^
4       F  pWCV,          8X  g   Q hWC,          pK  	  R# r   )	r   r   r&   r+   r   r-   r.   r/   rB   )r)   r1   ri   r   rC   rD   s         r4   test_DifferentialOperatorPowr   L  s    T"A !1!1!!4d;DA	q!tQVV[[!&&++>	BBaffjj\1-A2YEzz	 r6   N)Wsympy.holonomicr   r   r   r   r   r   sympy.holonomic.recurrencer	   r
   
sympy.corer   sympy.core.numbersr   r   r   sympy.core.singletonr   sympy.core.symbolr   r   &sympy.functions.elementary.exponentialr   r   %sympy.functions.elementary.hyperbolicr   r   (sympy.functions.elementary.miscellaneousr   (sympy.functions.elementary.trigonometricr   r   sympy.functions.special.besselr   &sympy.functions.special.beta_functionsr   'sympy.functions.special.error_functionsr   r   r   r   'sympy.functions.special.gamma_functionsr   sympy.functions.special.hyperr    r!   sympy.printing.strr"   sympy.series.orderr#   sympy.simplify.hyperexpandr$   sympy.polys.domains.integerringr%   !sympy.polys.domains.rationalfieldr&   sympy.polys.domains.realfieldr'   r5   r<   r@   rE   rV   rZ   r]   ra   rc   rj   rp   rx   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r6   r4   <module>r      s   > > N ! 0 0 " / = ? 9 ? 2 7 G G 9 : #   2 . 0 ,"D4 1F#!L4*$&8"5HKZK\F*P4;|*BZ_"%>\4
		 r6   