find where 2 lines cross given hough transform parameters in matlab

main.m
% Read the image and convert it to grayscale
img = imread('path/to/image.jpg');
bw = rgb2gray(img);

% Apply edge detection
edges = edge(bw);

% Apply Hough Transform to obtain the parameters of the lines
[H,T,R] = hough(edges);
P = houghpeaks(H,2);
lines = houghlines(edges,T,R,P);

% Find the intersection point of the two lines
x1 = lines(1).point1(1);
y1 = lines(1).point1(2);
x2 = lines(1).point2(1);
y2 = lines(1).point2(2);
slope1 = (y2-y1)/(x2-x1);
intercept1 = y1-slope1*x1;

x3 = lines(2).point1(1);
y3 = lines(2).point1(2);
x4 = lines(2).point2(1);
y4 = lines(2).point2(2);
slope2 = (y4-y3)/(x4-x3);
intercept2 = y3-slope2*x3;

x_intersect = (intercept2-intercept1)/(slope1-slope2);
y_intersect = slope1*x_intersect + intercept1;

% Plot the lines and the intersection point on the image
imshow(img), hold on
line([x1 x2],[y1 y2],'Color','r','LineWidth',2)
line([x3 x4],[y3 y4],'Color','g','LineWidth',2)
plot(x_intersect,y_intersect,'bo','MarkerSize',10,'LineWidth',2)
981 chars
36 lines