function vizualizeEpiLines(I0,I1,F)

figure(1); imshow(I0); hold on;
figure(2); imshow(I1); hold on;
X1 = size(I0,2); % number of columns = size in X (horizontal)
Y1 = size(I0,1); % number of rows = size in Y (vertical)
X2 = size(I1,2); % number of columns = size in X (horizontal)
Y2 = size(I1,1); % number of rows = size in Y (vertical)
while 1
    figure(1);
    [x,y,b]=ginput(1); % wait for key / mouse
    plot(x,y,'rx'); % plot the point
    figure(2);
    l = F*[x;y;1]; % compute the epipolar line equation   
    % intersect line l with the left and right boundaries of the image
    qleft = [0;-l(3)/l(2)];
    qright = [X2;-(l(3)+X2*l(1))/l(2)];
    plot([qleft(1) qright(1)],[qleft(2) qright(2)],'r-'); % plot it
end;