? 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)
%41 = 1
? m=3981;s=m;for(X=1,6,s=s^4-s^2+1);print(m,” “,round(log(s)/log(10)))
3981 14746
? ispseudoprime(s,30)
%42 = 1
