Lucas-Lehmer test: Difference between revisions
Content deleted Content added
→{{header|Pascal}}: add example |
adding maxima |
||
Line 618:
3 7 31 127 8191 131071 524287</lang>
=={{header|Maxima}}==
<lang maxima>lucas_lehmer(p) := block([s, n, i],
if not primep(p) then false elseif p = 2 then true else
(s: 4,
n: 2^p - 1,
for i: 2 thru p - 1 do s: mod(s*s - 2, n),
is(s = 0))
)$
sublist(makelist(i, i, 1, 200), lucas_lehmer);
/* [2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127] */</lang>
=={{header|Modula-3}}==
| |||