Neil Gunther has pointed out on his blog that the prime number decomposition R package schoolmath contains mistakes in the function primes, listing 1 as a prime number but also including decomposable numbers like 133 in its list of prime numbers: __ primes(100,140) 101 107 111 113 123 129 131 137 __ primes(50,140)