.PS include(libcct.m4) cct_init # Usual defs... qrt=dimen_/4; hlf=dimen_/2; dim=dimen_; mm=1/25.4; pi=atan2(0,-1); 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 2*pi 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 2*pi by 0.01 do { psi = pi/2*cos(theta) + delta; e=abs(factor*(sin(n*psi/2)/sin(psi/2))); line to L1.center+(e*cos(theta),e*sin(theta)) } "n=9" at L1.center+(0,1.3*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 2*pi by 0.01 do { psi = pi/2*cos(theta) + delta; e=abs(factor*(sin(n*psi/2)/sin(psi/2))) line to L1.center+(e*cos(theta),e*sin(theta)) } "n=12" at L1.center+(0,1.3*dim) "\tiny UniformEnd1" at (5.5*dim,-dim) .PE