chern.lib
LIB "chern.lib"; ring r=0, (c(1..5), t), ds; list l=c(1..5); // get the Chern polynomial 1 + c(1)*t + c(2)*t^2 + ... + c(5)*t^5 print( chernPoly(l, t) ); ==> 1+c(1)*t+c(2)*t^2+c(3)*t^3+c(4)*t^4+c(5)*t^5