// 08.06.05 // 16.12.29 ( for too short segment ) function Out=Anglemark(varargin) Nargs=length(varargin); Eps=10^(-3); PA=varargin(1); PB=varargin(2); PC=varargin(3); r = 0.5; if Nargs>=4 r = varargin(4)*r; end Out=[]; // 16.12.29 if r>min(norm(PA-PB),norm(PC-PB)) then return; end; Cir=Circledata(PB,r); Tmp=IntersectcrvsPp(Cir,Listplot(PA,PB)); P1=Op(2,Op(1,Tmp)); Tmp=IntersectcrvsPp(Cir,Listplot(PC,PB)); P2=Op(2,Mixop(1,Tmp)); if abs(P1-P2)