polar plot r= 1+ sin(pi^2 / phi), phi from 0 to 2pi