V=matrix(3, 3, [1/6*sqrt(6), 2/6*sqrt(6), 1/6*sqrt(6), 1/2*sqrt(2), 0, -1/2*sqrt(2), 1/3*sqrt(3), -1/3*sqrt(3), 1/3*sqrt(3)])
U=matrix(2, 2, [1/2*sqrt(2), -1/2*sqrt(2), 1/2*sqrt(2), 1/2*sqrt(2)])
print V.transpose()*F.augment(y).transpose()*U.transpose()
|
[1/36*sqrt(2)*sqrt(3)*sqrt(6) - 1/2 1/36*sqrt(2)*sqrt(3)*sqrt(6) + 1/2]
[ 1/18*sqrt(2)*sqrt(3)*sqrt(6) 1/18*sqrt(2)*sqrt(3)*sqrt(6)]
[1/36*sqrt(2)*sqrt(3)*sqrt(6) + 1/2 1/36*sqrt(2)*sqrt(3)*sqrt(6) - 1/2]
[1/36*sqrt(2)*sqrt(3)*sqrt(6) - 1/2 1/36*sqrt(2)*sqrt(3)*sqrt(6) + 1/2]
[ 1/18*sqrt(2)*sqrt(3)*sqrt(6) 1/18*sqrt(2)*sqrt(3)*sqrt(6)]
[1/36*sqrt(2)*sqrt(3)*sqrt(6) + 1/2 1/36*sqrt(2)*sqrt(3)*sqrt(6) - 1/2]
|