#Please do not post homework
#Jeffrey Tang, 2/8/25, Assignment 4

with(linalg):

#UM:= inputs alpha, beta, gamma, delta and outputs a unitary matrix
UM := (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/2), 0, 0, exp(I*delta/2)])
):