function [numVfacets] = Facetchecker1(x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x15, x16, x17, x18, x19, x20, n1, n2, n3, n4, n5, n6, n7, n8, n9, n10, n11, n12, n13, n14, n15, n16, n17, n18, n19, n20, A) A1 = A(:, 1:20); b = 1*A(:,21); p = zeros(size(A,1),1); obspoint = [x1/n1; x2/n2; x3/n3; x4/n4; x5/n5; x6/n6; x7/n7; x8/n8; x9/n9; x10/n10; x11/n11; x12/n12; x13/n13; x14/n14; x15/n15; x16/n16; x17/n17; x18/n18; x19/n19; x20/n20]; i = 1; G = A1*obspoint; while i<=size(A,1) if (G(i,1) > b(i,1)); p(i,1) = 1; else p(i,1) = 0; end i = i + 1; end p; numVfacets = sum(p);