After I learned all the details in grad school, I found out how to implement it in a real world setting:

  [U,S,V] = svd(X)