function ang=beta4(i1,i2,i3,g,f,dir,si,ni,pn)
%計算i2處的左夾角
[dir1,kdir1]=finddir(i2,i1,g,f,dir,si,ni);%判斷該往測方向觀測值是否存在并返回之
[dir2,kdir2]=finddir(i2,i3,g,f,dir,si,ni);%判斷該反測方向觀測值是否存在并返回之
k=kdir1+kdir2;
if k==2
	ang=dir2-dir1;
	if ang>2*pi
        ang=ang-2*pi;
	elseif ang<0
        ang=ang+2*pi;
	end
 else
     ang=[];
end
return;