TODO
 - write manual
 - change init seed (incl. PHI)
 - make function which chooses a,b,c based on period