use (r+1)*(r+1) so that radius=0 works