def normal(U,V,N): p1=0;p2=0 for i in range(3): p1=p1+U[i]*N[i] p2=p2+V[i]*N[i] rep=(p1==0 ... p2==0) return(rep)