#!conda deactivate
!conda update anaconda-navigator

!conda update conda

!conda update conda
!conda update anaconda

import numpy as np

a = np.array([[ 1.,  0.,  0.,  0.],
       [ 0.,  1.,  0.,  0.],
       [ 0.,  0.,  0.,  1.],
       [ 0.,  0.,  1.,  0.]])

p = q =r= s = 1

b = np.array([[ p,  0.,  q,  0.],
       [ 0.,  0.,  0.,  0.],
       [ r,  0.,  s,  0.],
       [ 0.,  0.,  0.,  0.]])

a@b

a@b@a

