function A=ar9(c1_x,c1_y, cj_x,cj_y, c2_x,c2_y, d, pt_x,pt_y,G1_y)
%global G1_y
bf = angle(cj_x,cj_y, pt_x,pt_y) - 1.5*pi;
cf = abs(cos(bf-0.5*(angle(cj_x,cj_y,c1_x,c1_y)+angle(c2_x,c2_y,cj_x,cj_y))));
	gc = cj_y - G1_y;
	gw = 5.04 - G1_y;
	A = 2.2 * (gc/gw) + 2.0;
	A = A * cf * d * pi/2.0;