import numpy as np

import numpy.linalg as LA

import matplotlib.pyplot as plt

rng = np.random.default_rng()

```

Having imported these modules, let's do the example.

```{code-cell} ipython3

import numpy as np

X = np.random.rand(5,2)

X = rng.random((5, 2))

U, S, V = np.linalg.svd(X,full_matrices=True) # full SVD

Uhat, Shat, Vhat = np.linalg.svd(X,full_matrices=False) # economy SVD

print('U, S, V =')

```{code-cell} ipython3

import numpy as np

X = np.random.rand(2,5)

X = rng.random((2, 5))

U, S, V = np.linalg.svd(X,full_matrices=True) # full SVD

Uhat, Shat, Vhat = np.linalg.svd(X,full_matrices=False) # economy SVD

print('U, S, V = ')