#Please do not post homework
#Salman Manzoor, 02/09/25, Assignment 4

with (linalg):
with (combinat):
UM:=proc(alpha,beta,gamma,delta) 
evalm( exp(I*alpha) *
matrix(2,2,[exp(-I*beta/2),0,0,exp(I*beta/2)])  &*
matrix(2,2,[cos(gamma/2) , -sin(gamma/2) , sin(gamma/2) , cos(gamma/2)]) &*
matrix(2,2,[exp(-I*delta) , 0,0,exp(I*delta)]) ):
end: