TPISDC13 := (1+12*k-12*m)*z+(-36-262*k+80*m)*z^2+(1072*k+386-526*m)*z^3+(-883-1632*k+280*m)*z^4+(628*k-576*m+1043)*z^5+(-312*m-598*k-82)*z^6+(-312*k-318*m+456)*z^7+(-318*k-168*m+315)*z^8+(-168*k-68*m+243)*z^9+(-68*k-16*m+118)*z^10+(42-16*k-2*m)*z^11+(9-2*k)*z^12+z^13: TPISDC13:=expand(TPISDC13/z): st:=time(): ANS:=resultant(TPISDC13,diff(TPISDC13,z),z): time()-st;