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

RandomState:=proc(n,K)
local ra,i,k,L:
ra:=rand(-K..K):
L:=[seq(ra()+I*ra(),i=1..n)]:
k:=sqrt(add(abs(L[i])^2,i=1..n)):
[seq(L[i]/k,i=1..n)]:
end:

GeneralPauli:=proc(n) local L,i:
L:=PM():
evalm(n[1]*matrix(L[1])+n[2]*matrix(L[2])+n[3]*matrix(L[3])):
end:

#PM(): the three famous Pauli matrices sigma_x,sigma_y, sigma_z
PM:=proc():
[
[[0,1],[1,0]],
[[0,-I],[I,0]],
[[1,0],[0,-1]]
]:
end: