こちらでも利用できるのですね！！こちらに、
[n : n in [2..500] | IsPrime(n)];
を放りこんでみたりしました。
（コメント） 以下が、出力結果です。２〜５００の範囲での素数一覧ですね。
[ 2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101,
103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197,
199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271, 277, 281, 283, 293, 307, 311,
313, 317, 331, 337, 347, 349, 353, 359, 367, 373, 379, 383, 389, 397, 401, 409, 419, 421, 431,
433, 439, 443, 449, 457, 461, 463, 467, 479, 487, 491, 499 ]
Dengan kesaktian Indukmu さんからのコメントです。（令和６年１２月２３日付け）
[(n+1)*(n+2)-(n+1)*(-1)^n: n in [43..44]]
をほうりこむと、[ 2024, 2025 ] が得られますね。
kuiperbelt さんからのコメントです。（令和６年１２月２３日付け）
「ＭＡＧＭＡ」は私も時々使いますが、いろいろな機能がありますね。
○ G2型リー代数の表現の次元数について、
R := RootDatum("G2");
for i:=0 to 4 do;
for j:=0 to i do;
i,j,RepresentationDimension(R, [j,i-j]);
end for;
end for;
○ PSL(2,9)の指標表について
G := PSL(2,9);
CharacterTable(G);
○ 方程式 x^6+2*x^5+3*x^4+4*x^3+5*x^2+6*x+7=0 のガロア群について
P<x> := PolynomialRing(Integers());
f:=x^6+2*x^5+3*x^4+4*x^3+5*x^2+6*x+7;
G, L, S := GaloisGroup(f);
G;
S;
○ 3-進数体での-7/2の平方根について
K := pAdicField(3,40);
_<x> := PolynomialRing(Integers(K)); // printing
HasRoot(2*x^2+7);
K`SeriesPrinting := true;
Sqrt(K!(-7/2));
○ 二次体Q(√-31)でのノルムが100未満の素イデアルについて
R<x> := PolynomialRing(Integers());
K := NumberField(x^2+x+8);
FactorBasis(K, 100);
以下、工事中！
