Some Rabin-Miller pseudoprime primality tests

? m=55;s=m;for(X=1,6,s=s^4-s^2+1);print(m,” “,round(log(s)/log(10)))55 7128 ? ispseudoprime(s,30)%39 = 1 ? m=1948;s=m;for(X=1,6,s=s^4-s^2+1);print(m,” “,round(log(s)/log(10)))1948 13474 ? ispseudoprime(s,30)%40 = 1 ? m=3269;s=m;for(X=1,6,s=s^4-s^2+1);print(m,” “,round(log(s)/log(10)))3269 14395 ? ispseudoprime(s,30)…