RosettaCodeData/Task/Fermat-numbers/Maxima/fermat-numbers.maxima

fermats(n) := makelist(2^(2^i)+1, i, 0, n)$
fermats(9);
map(ifactors, fermats(7));