PARI GP ispseudoprime


PRTWINSB=b

aa=11*13*17*19*101*103*107*109*191*193*197*199

default(realprecision,20000)
b=floor(aa^456.84);

To be completed...