for 0<p<150 plot Im((a^p)%p)