model ("inv-v uda-yagi") {

real dist, height, length, angle, eps, gauge, segments, a, b, r_a, r_b, r_coeff;
element center, left, right, r_center, r_left, r_right ;

length = 5.10;
// 4.9 - gain: 11.48, F/B: 9.85
// 4.5 - gain: 11.57, F/B: 9.80
// 4.0 - gain: 11.68, F/B: 9.60
dist = -4.5;
angle = 125;
r_coeff = 1.04;
height = 10;

eps = 0.01;
gauge = #20;
segments = 21;

b = length*sin(pi*angle/2/180);
a = length*sin(pi*(90 - angle/2)/180);
r_b = b*r_coeff;
r_a = a*r_coeff;


center = wire(0,-eps, height,
              0, eps, height,
              gauge, segments);
left = wire(0, -b, height - a,
            0,-eps, height,
            gauge, segments);
right = wire(0, b, height - a,
             0,eps,height,
             gauge, segments);
             
// reflector
r_center = wire(dist,-eps, height,
              dist, eps, height,
              gauge, segments);
r_left = wire(dist, -r_b, height - r_a,
            dist,-eps, height,
            gauge, segments);
r_right = wire(dist, r_b, height - r_a,
             dist,eps,height,
             gauge, segments);
             
voltageFeed(center, 1.0, 0.0);

azimuthPlotForElevationAngle(31);
elevationPlotForAzimuthAngle(0);

frequencySweep(14.0, 14.35, 10);
//frequencySweep(21.0, 21.45, 10);
//setFrequency(14.070);

// poorGround();
averageGround();
// goodGround();
}