p := vector([1,1,1]): a:=vector([1,1,0]): n := vector([0,0,1]): dist := innerprod(n,p-a)/sqrt(innerprod(n,n)); p2plane := proc(p,a,n) innerprod(n,p-a)/sqrt(innerprod(n,n)); end; p2plane(p,a,n); p2plane([-1,sqrt(2),5],[-1,0,123],[-1,3,0]); evalf(");