Some new GP features
polylogmult
The function polylogmult is an extension of zetamult to
compute multiple polylogarithm values.
? polylogmult([2,2,2],[1,-1,1])
%47 = 0.071635745321742507017850817313956412547