f(x):=(x-h)^2+k;

f(a);

[a,f(a)];

m(x):= diff(f(x),x)$ m(x);

m(x):=2*(x-h);

y-f(a) = m(a)*(x-a);

first(ratsimp(solve(%,y)));