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.