p=input("p premier = ")
s=4;M=2^p-1
for k=1:(p-2)
    s=reste(s^2-2,M)
end
if s==0 then afficher("Mp premier")
else afficher("Mp non premier")
end