Hi, I got a problem  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 <ginac/ginac.h></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>    symbol x1("x1"),x2("x2"),x3("x3"),x4("x4");</div><div>    symbol ep("ep");</div><div>    ex expr = -2*pow(x1,-1 + ep)*(exp(3*ep*Euler)*pow(ep,-1)*pow(x2,-3*ep)*</div><div>      pow(1 + x3 + x4 + x3*x4,-3*ep)*</div><div>      pow(1 + x2 + x3 + x2*x3 + x4 + x2*x4 + x3*x4 + x2*x3*x4,-2 + 4*ep)*tgamma(3*ep) -</div><div>     exp(3*ep*Euler)*pow(ep,-1)*pow(x2,-3*ep)*pow(1 + x3 + x4 + x3*x4,-3*ep)*</div><div>      pow(1 + x2 + x3 + x1*x3 + x2*x3 + x4 + x2*x4 + x3*x4 + x1*x3*x4 + x2*x3*x4,</div><div>       -2 + 4*ep)*tgamma(3*ep));</div><div>    </div><div>    auto s1 = expr.series(ep,-1);</div><div>    cout << s1.normal() << endl;</div><div>    </div><div>    auto s2 = expr.series(ep,0);</div><div>    cout << s2.normal() << endl;</div><div>    </div><div>    return 0;</div><div>}</div></div>