function [x_i,y_i] = line_intersect_a(m1,b1,m2,b2)
% function [x_i,y_i] = line_intersect_a(m1,b1,m2,b2)
% This function checks to see if two lines,
% both in the plane z = 0, intersect.
% If they do, it returns the coordinates of the
% point of intersection and plots the point.
% If the lines don't intersect, it returns values of
% NaN ("Not a Number") for x and y, and no plot
% is produced
% Input parameters
% m1 = slope of line 1
% b1 = y-intercept of line 1
% m2 = slope of line 2
% b2 = y-intercept of line 2
%
% The equations of the lines are
% y = m1.*x + b1
% y = m2.*x + b2
% Example
% [x_i,y_i] = line_intersect_a(1,1,2,2)
% Name: Steve Martel
% Date: 9/15/04

 
% Check to see if lines intersect.
% The lines do not intersect if the slopes are the same.
if m1 == m2
    % The lines are parallel and do not intersect
    x = NaN;
    y = NaN;
    return
else
    % The lines are not parallel and therefore intersect
    % Find the coordinates of the intersection
    x_i = (b2 - b1)./(m1 - m2);
    y_i = m1.*x_i + b1;

    % Plot the intersection, extending the lines a bit from that point
    % Prepare two subplots, the first showing the x- and y- intercepts
    % of the lines, the second showing their intersection
    x_int_1 = -b1./m1;
    x_int_2 = -b2./m2;

    figure(1)
    clf
    subplot(1,2,1)
    plot([x_int_1 0],[0 b1],'b',[x_int_2 0],[0 b2],'r')
    title (['Lines intersect at x = ',num2str(x_i),' y = ',num2str(y_i)])
    xlabel('x')
    ylabel('y')
    
    subplot(1,2,2)
    plot(x_i,y_i,'+')
    title (['Lines intersect at x = ',num2str(x_i),' y = ',num2str(y_i)])
    xlabel('x')
    ylabel('y')
end