.PS cct_init # Usual defs... qrt=dimen_/4; hlf=dimen_/2; dim=dimen_; L1:line right_ 2*dim delta = -pi_/2; # Endfire n = 4; factor=0.1; # d = lambda/4. beta d = pi_/2 # get rid of first point. move to L1.center for theta = 0.01 to twopi_ by 0.01 do { psi = pi_/2*cos(theta) + delta; e = abs(factor*(sin(n*psi/2)/sin(psi/2))); if theta==0.01 then{ move to L1.center+(e*cos(theta),e*sin(theta)) }else{ line to L1.center+(e*cos(theta),e*sin(theta)) } } "n=4" at L1.center+(0,1.3*dim) move to L1.end+(dim,0) L1:line right_ 2*dim n = 9; factor=0.075; move to L1.center for theta = 0.01 to twopi_ by 0.01 do { psi = pi_/2*cos(theta) + delta; e=abs(factor*(sin(n*psi/2)/sin(psi/2))); if theta==0.01 then{ move to L1.center+(e*cos(theta),e*sin(theta)) }else{ line to L1.center+(e*cos(theta),e*sin(theta)) } } "n=9" at L1.center+(0,1.3*dim) "$d=\lambda/4; \beta d=\pi/2; \delta=-\pi/2$" at L1.center+(0,-dim) move to L1.end+(1.5*dim,0) L1:line right_ 2*dim n = 12; factor=0.075; move to L1.center for theta = 0.01 to twopi_ by 0.01 do { psi = pi_/2*cos(theta) + delta; e=abs(factor*(sin(n*psi/2)/sin(psi/2))) if theta==0.01 then{ move to L1.center+(e*cos(theta),e*sin(theta)) }else{ line to L1.center+(e*cos(theta),e*sin(theta)) } } "n=12" at L1.center+(0,1.3*dim) "\tiny UniformEnd1" at (7.5*dim,-dim) .PE