Matlab - Dado dos puntos encontrar la recta y dibujar

%% Datos

% Puntos instertados por usuario

x1 = input('Inserte x1 = ');

y1 = input('Inserte y1 = ');

x2 = input('Inserte x2 = ');

y2 = input('Inserte y2 = ');

% Vector X e Y

x = [x1, x2];

y = [y1, y2];

%% Programa

% Coeficientes de regresión

k = polyfit(x,y,1);

% Vector X preciso para gráfico

x_fit = linspace(min(x),max(x));

% Ecuación de regresión

y_fit = k(1)*x_fit + k(2);

%% Gráfico

figure;

plot(x_fit,y_fit,'LineWidth',2);

hold on;

plot(x1,y1,'x','LineWidth',2);

plot(x2,y2,'x','LineWidth',2);

grid on;

xlabel('\itx');

ylabel('\ity');

legend('y = m*x + k','(x_1, y_1)','(x_2, y_2)');