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);