Solve for S
•
SVD is made to do this.
D is diagonal with non-increasing
values.
U and V have orthonormal rows.
Ignoring values that get set to 0, we
have U(:,1:3) for S, and
D(1:3,1:3)*V(1:3,
:) for P.