from qiskit import *
import numpy as np

circuito = QuantumCircuit(2,2)
circuito.x(0)
circuito.z(0)
circuito.h(1)
# circuito.measure([0,1],[0,1])
circuito.draw(output = 'mpl')

# Podemos encontrar a matriz unitária que representa as operacoes no circuito

backend = Aer.get_backend('unitary_simulator')
job = execute(circuito, backend)
resultado = job.result()
unitaria = resultado.get_unitary()

print(unitaria)