Calculating the Posterior

# Calculating the posterior for prior > 0.5
# Specify the prior probability and the modifier
prior = 0.6
modifier = 1.2
qc = QuantumCircuit(4)
# Apply prior to qubit 0
qc.ry(prob_to_angle(prior), 0)
# Apply 0.3*prior to qubit 1
qc.x(0)
qc.cry(prob_to_angle(0.3*prior/(1-prior)), 0, 1)
# Apply the modifier to qubit 2
qc.cry(prob_to_angle((modifier-1)/0.3), 1,2)
# Make qubit 3 represent the posterior
qc.x(0)
qc.cx(0,3)
qc.cx(2,3)
run_circuit(qc)

qc.cry(prob_to_angle((modifier-1)/0.3), 1,2)