/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 13.04.2 ] */
/* [wxMaxima: comment start ]
Un
altro numero di Mersenne
[wxMaxima: comment end ] */
/* [wxMaxima: comment start ]
Questo primo : 74207281 che ha come piu' piccola
radice primitiva il numero 65, e' un numero speciale
perche' (-1+2^74207281) e' anche lui un numero primo,
come e' stato trovato recentissimamente, nel gennaio 2016.
[wxMaxima: comment end ] */
/* [wxMaxima: input start ] */
buonprimo:74207281$
print("(-1+2^74207281) dicono che e' primo!")$
print("Verifico che 74207281 sia un numero primo: ",
primep(buonprimo))$
print("74207281 ha 65 come radice primitiva: ",
zn_primroot(buonprimo))$
/* [wxMaxima: input end ] */
/* [wxMaxima: comment start ]
I numeri p che sono primi e sono usabili per definire altri
primi nella forma (-1+2^p) dove anche questa espressione e'
un numero primo, sono gli esponenti dei cosi' detti
primi di Mersenne. Vedere il sito:
http://www.mersenne.org/primes/
[wxMaxima: comment end ] */
/* [wxMaxima: input start ] */
pesp:[2,3,5,7,13,17,19,31,61,89,107,127,521,607,1279,
2203,2281,3217,4253,4423,9689,9941,11213,19937,
21701,23209,44497,86243,110503,132049,216091,
756839,859433,1257787,1398269,2976221,3021377,
6972593,13466917,20996011,24036583,25964951,
30402457,32582657,37156667,42643801,43112609,
57885161,74207281];
/* [wxMaxima: input end ] */
/* [wxMaxima: comment start ]
Ecco la verifica che tutti gli esponenti dei numeri
di Mersenne noti sono anche loro numeri primi.
[wxMaxima: comment end ] */
/* [wxMaxima: input start ] */
for j:1 step 1 thru length(pesp) do
(print(j," primo ",pesp[j]," Verifica: ",primep(pesp[j]),
" Radice primitiva : ",zn_primroot(pesp[j])));
/* [wxMaxima: input end ] */
/* [wxMaxima: comment start ]
https://it.wikipedia.org/wiki/Numero_perfetto
https://it.wikipedia.org/wiki/Numero_primo_di_Mersenne
[wxMaxima: comment end ] */
/* [wxMaxima: comment start ]
[wxMaxima: comment end ] */
/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$