Codes: part 20

a(n)=if(n=0),1+(trace(vecextract(factor(a(n-1)),”1″,”1″)))*(a(n-1)-1))
for(n=0,4,print(trace(vecextract(factor(a(n)),”1″,”1″))))

Advertisements