function pt = intersect(pt1,pt2,style) A = [diff(pt1)' -diff(pt2)']; b = pt2-pt1; q = A\b(1,:)'; pt = pt1(1,:) + q(1)*A(:,1)'; %pt = pt2(1,:) - q(2)*A(:,2)'; % solution based on second line if nargin>2 hold on plot(pt(1),pt(2),style); hold off end
function pt = intersect(pt1,pt2,style) A = [diff(pt1)' -diff(pt2)']; b = pt2-pt1; q = A\b(1,:)'; pt = pt1(1,:) + q(1)*A(:,1)'; %pt = pt2(1,:) - q(2)*A(:,2)'; % solution based on second line if nargin>2 hold on plot(pt(1),pt(2),style); hold off end