new_specification : string * string list * thm -> thm
new_specification (name, ["c1",...,"cn"], |- ?x1...xn. t)
|- t[c1,...,cn/x1,...,xn]
th |- ?MOD DIV. !n. 0 < n ==> !k. (k = (DIV k n * n) + MOD k n) /\ MOD k n < n
new_specification ("DIVISION", ["MOD","DIV"], th)