theorem :: XPRIMES1:4001

4001 is prime