pnatural is a Prolog module that implements natural number arithmetic by using Peano axioms. The natural numbers are represented like s(s(z)) ("successor form"), where s(X) means "successor of X" and z means "zero". s(s(z)) represents the natural number two.
For more details and download, visit its github repo.