function F = sevenp(x1,x2) Q1 = x1'; Q2 = x2'; Q = [Q1(:,1).*Q2(:,1) , ... Q1(:,2).*Q2(:,1) , ... Q1(:,3).*Q2(:,1) , ... Q1(:,1).*Q2(:,2) , ... Q1(:,2).*Q2(:,2) , ... Q1(:,3).*Q2(:,2) , ... Q1(:,1).*Q2(:,3) , ... Q1(:,2).*Q2(:,3) , ... Q1(:,3).*Q2(:,3) ] ; %EE = null( Q); [U,S,V] = svd(Q,0); F = V(:,9); F = F/norm(F);