mata function iseven(n) { return(mod(n,2)==0) } function isodd(n) { return(mod(n,2)==1) } end