The spec of math/big.modInverse requires the modulus to be prime. But the inverse is
well defined as long as the number being inverted and the modulus are relatively prime.
Perhaps more importantly, the current code seems to work for all of those cases.
Should we relax the spec to include these cases?
Any reason why the current code wouldn't always work with this relaxation?