convert line to hesse normal form and plot it in matlab

main.m
% Input values
m = 2; % slope
b = 5; % y-intercept

% Convert to Hesse normal form
theta = atan(-m); % angle between x-axis and line
r = b / cos(theta); % perpendicular distance from origin to line
s = b / sin(theta); % perpendicular distance from origin to line
d = r * cos(theta); % signed distance from origin to line
hnf = [r, s, d]; % Hesse normal form [r, s, d] 

% Plot the line
x = linspace(-10, 10); % create range of x values
y = m * x + b; % calculate corresponding y values
plot(x, y, 'LineWidth', 2) % plot the line
hold on % hold current plot
x = linspace(-10, 10); % create range of x values
y = (d - r*cos(theta))/sin(theta)*ones(1,length(x)); % calculate corresponding y values for Hesse normal form
plot(x, y, '--', 'LineWidth', 2) % plot the line in Hesse normal form
axis equal % set aspect ratio to 1:1
824 chars
21 lines