Matlab - funcion ecuacion recta

f='y=m*x+n'

subs(f,'m',tan(alfa))

a partir de aqui estoy perdido, he probado con subs para x0 e y0 y con solve pero no consigo hacerlo bien.

syms x y

xo=input('ingrese coordenada xo = ');

yo=input('ingrese coordenada yo = ');

alpha=input('ingrese angulo alpha: ');

m=tand(alpha); % TANGENTE EN SEXAGESIMALES

y=m*(x-xo)+yo;

fprintf('la ecuacion de la recta es: y = %s\n',char(y));

clear all

syms y1 x

y1=2;

y2=x+1;

figure(gcf)

plot([-5 5],[ 2 2],[-5 5],[-5 5]+1)

hold on

X=solve(y1==y2);

plot(X,2,'*r')

hold off

EJECUCION

>> interseccion_recta

clear all; clc

syms x

y1=x+1;

y2=2;

if isnumeric(y1)==1

    y1=num2str(y1);

end

if isnumeric(y2)==1

    y2=num2str(y2);

end

xp=str2double(char(solve(y1==y2)));

ec=inline(y2);

yp=ec(xp);

fplot(char(y1),[xp-5, xp+5],'b')

hold on

fplot(char(y2),[xp-5, xp+5],'r')

grid

plot(xp,yp,'*g'),text(xp+0.2,yp+0.2,['(',num2str(xp),',',num2str(yp),')'])

legend(['y = ',char(y1)],['y = ',char(y2)])

hold off