A = randn(5,3)

U, σ, V = svd(A)

σ.^2    # the σ² values

eigvals(A'*A) # AᴴA has the same eigenvals!

eigvals(A*A') # the same *nonzero* eigenvalues

eigvecs(A'A)

V