Hi, I got a problem &nbsp;about the series function,<div>I used the following test code(see below or the attached file),</div><div>the s1 output only gives Order(ep^(-1)), while s2 output indicates that</div><div>there is a 1/ep^2 term, so the s1 is not correct or I used the series in a wrong way?</div><div><br></div><div>Thanks a lot!</div><div><br></div><div><div>#include &lt;ginac/ginac.h&gt;</div><div><br></div><div>using namespace GiNaC;</div><div>using namespace std;</div><div><br></div><div>int main(int argn, char** args) {</div><div>&nbsp; &nbsp; symbol x1("x1"),x2("x2"),x3("x3"),x4("x4");</div><div>&nbsp; &nbsp; symbol ep("ep");</div><div>&nbsp; &nbsp; ex expr = -2*pow(x1,-1 + ep)*(exp(3*ep*Euler)*pow(ep,-1)*pow(x2,-3*ep)*</div><div>&nbsp; &nbsp; &nbsp; pow(1 + x3 + x4 + x3*x4,-3*ep)*</div><div>&nbsp; &nbsp; &nbsp; pow(1 + x2 + x3 + x2*x3 + x4 + x2*x4 + x3*x4 + x2*x3*x4,-2 + 4*ep)*tgamma(3*ep) -</div><div>&nbsp; &nbsp; &nbsp;exp(3*ep*Euler)*pow(ep,-1)*pow(x2,-3*ep)*pow(1 + x3 + x4 + x3*x4,-3*ep)*</div><div>&nbsp; &nbsp; &nbsp; pow(1 + x2 + x3 + x1*x3 + x2*x3 + x4 + x2*x4 + x3*x4 + x1*x3*x4 + x2*x3*x4,</div><div>&nbsp; &nbsp; &nbsp; &nbsp;-2 + 4*ep)*tgamma(3*ep));</div><div>&nbsp; &nbsp;&nbsp;</div><div>&nbsp; &nbsp; auto s1 = expr.series(ep,-1);</div><div>&nbsp; &nbsp; cout &lt;&lt; s1.normal() &lt;&lt; endl;</div><div>&nbsp; &nbsp;&nbsp;</div><div>&nbsp; &nbsp; auto s2 = expr.series(ep,0);</div><div>&nbsp; &nbsp; cout &lt;&lt; s2.normal() &lt;&lt; endl;</div><div>&nbsp; &nbsp;&nbsp;</div><div>&nbsp; &nbsp; return 0;</div><div>}</div></div>