W:=Matrix([[1,1,1,1],[1,1,1,1],[1,1,1,1],[1,1,1,1]]);
P1:=Matrix(4, 4, {(1, 1) = [-3., -6., 0.], (1, 2) = [-1., -6., 4.], (1, 3) = [1., -6., 4.], (1, 4) = [3., -6., 0.], (2, 1) = [-3., -2., 0.], (2, 2) = [-1., -2., 4.], (2, 3) = [1., -2., 4.], (2, 4) = [3., -2., 0.], (3, 1) = [-3., 2., 0.], (3, 2) = [-1., 2., 4.], (3, 3) = [1., 2., 4.], (3, 4) = [3., 2., 0.], (4, 1) = [-3., 6., 0.], (4, 2) = [-1., 6., 4.], (4, 3) = [1., 6., 4.], (4, 4) = [3., 6., 0.]});
P2:=Matrix(4, 4, {(1, 1) = [-6., -3., 0.], (1, 2) = [-2., -3., 0.], (1, 3) = [2., -3., 0.], (1, 4) = [6., -3., 0.], (2, 1) = [-6., -1., 4.], (2, 2) = [-2., -1., 4.], (2, 3) = [2., -1., 4.], (2, 4) = [6., -1., 4.], (3, 1) = [-6., 1., 4.], (3, 2) = [-2., 1., 4.], (3, 3) = [2., 1., 4.], (3, 4) = [6., 1., 4.], (4, 1) = [-6., 3., 0.], (4, 2) = [-2., 3., 0.], (4, 3) = [2., 3., 0.], (4, 4) = [6., 3., 0.]});