TPISDC17 := (1+16*k-16*m)*z+(-64-632*k+224*m)*z^2+(1240+5376*k-2520*m)*z^3+(-6328+5440*m-17068*k)*z^4+(14332-9264*m+21232*k)*z^5+(-11760+496*m-15864*k)*z^6+(1432*k+9128-6124*m)*z^7+(-6154*k-4896*m+2240)*z^8+(-4496*m+5526-4896*k)*z^9+(-4496*k+4696-2880*m)*z^10+(-2880*k+3688-1464*m)*z^11+(-1464*k-544*m+2172)*z^12+(-544*k+1004-144*m)*z^13+(344-24*m-144*k)*z^14+(-2*m-24*k+84)*z^15+(13-2*k)*z^16+z^17: TPISDC17:=expand(TPISDC17/z): st:=time(): ANS:=resultant(TPISDC17,diff(TPISDC17,z),z): time()-st;