Intersects

It finds intersect points between lines using arithmetic geometry formulas. to find all intersects run it like following;

here, x1,x2,x3,x4 are defined as;

x1=rand(1,100); y1=rand(1,100); x2=rand(1,100); y2=rand(1,100);

x=[];
y=[];
tic
for i = 1:100
x2a = circshift(x2,i);
y2a = circshift(y2,i);
[xa,ya]=Intersects(x1,y1,x2a,y2a);
x=[x,xa];
y=[y,ya];
end
toc
plot(x1,y1,'g',x2,y2,'b',x,y,'ko'),