Some new GP features foreach, parforeach foreach iterate over the elements of a vector. parforeach is the parallel version. ? V = binomial(10) %1 = [1,10,45,120,210,252,210,120,45,10,1] ? foreach(V,n,print1(sigma(n)," ")) 1 18 78 360 576 728 576 360 78 18 1 ? my(s=0);parforeach(V,n,sigma(n),S,s+=S);s %3 = 2794