[x,y]=meshgrid(-3:.1:3);
f=sqrt(20-x.^2-y.^2);
figure(1);
surf(x,y,f);
hold on;
F=sqrt(20-1-2^2);
L = F -1/F*(x-1)-2/F*(y-2);
surf(x,y,L);
hold off;