P3points := proc(a,b,c) local ex1,ex2,n,_x,_y,_z; n := crossprod(b-a, c-a); r := vector([_x,_y,_z]); ex1 := innerprod(n, r-a); ex2 := collect(ex1,[_x,_y,_z]); ex2 = 0; end;