W:=Matrix([[1.01,0.99,1,1.1],[1,1,1,0.98],[1,1,1,1],[1,1,1,1]]);
P1:=Matrix(4, 4, {(1, 1) = [1, -1, -1], (1, 2) = [1, -1, -1/3], (1, 3) = [1, -1, 1/3], (1, 4) = [1, -1, 1], (2, 1) = [-1/3, 1, -1], (2, 2) = [-1/3, 1, -1/3], (2, 3) = [-1/3, 1, 1/3], (2, 4) = [-1/3, 1, 1], (3, 1) = [-1/3, -1, -1], (3, 2) = [-1/3, -1, -1/3], (3, 3) = [-1/3, -1, 1/3], (3, 4) = [-1/3, -1, 1], (4, 1) = [1, 1, -1], (4, 2) = [1, 1, -1/3], (4, 3) = [1, 1, 1/3], (4, 4) = [1, 1, 1]});
P2:=Matrix(4, 4, {(1, 1) = [1, -1, -1], (1, 2) = [1, -1/3, -1], (1, 3) = [1, 1/3, -1], (1, 4) = [1, 1, -1], (2, 1) = [-1/3, -1, 1], (2, 2) = [-1/3, -1/3, 1], (2, 3) = [-1/3, 1/3, 1], (2, 4) = [-1/3, 1, 1], (3, 1) = [-1/3, -1, -1], (3, 2) = [-1/3, -1/3, -1], (3, 3) = [-1/3, 1/3, -1], (3, 4) = [-1/3, 1, -1], (4, 1) = [1, -1, 1], (4, 2) = [1, -1/3, 1], (4, 3) = [1, 1/3, 1], (4, 4) = [1, 1, 1]});