euclid_mullin(n):=if n=1 then 2 else ifactors(1+product(euclid_mullin(i),i,1,n-1))[1][1]$

/* Test case */
makelist(euclid_mullin(k),k,16);