add naive generator function primegen, uses trial division on ints of form 6*k(+/-)1