FArcEllipse.m
33bba04159d6master

FArcEllipse.m
View Options

			classdef FArcEllipse < FEllipse
			properties
			angles
			end
			methods
			function obj=FArcEllipse(a,b,angles,C)
			obj = obj@FEllipse(a,b,C);

			obj.angles = zeros(1,2);

			% shift form -pi to pi
			obj.angles(1) = rem(angles(1), 2*pi);
			obj.angles(2) = rem(angles(2), 2*pi);

			while obj.angles(2) < obj.angles(1)
			obj.angles(2) = obj.angles(2) + 2 * pi;
			end

			% adjust boundaries
			if obj.angles(2) < obj.angles(1)
			obj.angles(2) = obj.angles(2) + 2 * pi;
			end
			end

			function flag=inside(obj,x,y)
			flag = inside@FEllipse(obj,x,y);
			ang = atan2((y - obj.Cy) / obj.b, (x - obj.Cx) / obj.a);
			while ang < obj.angles(1)
			ang = ang + 2 * pi;
			end
			flag = flag & (ang < obj.angles(2));
			end
			end
			end