import matplotlib.pyplot as plt
import numpy as np

if __name__=='__main__':
    physical_error_rate=np.array([10**-3,10**-2.5,10**-2])
    logical_failure_rate=np.array([0.00042/5,0.005584/5,0.080752/5])
    fig, ax = plt.subplots(figsize=(8, 5))
    ax.plot(physical_error_rate, logical_failure_rate)
    ax.plot(physical_error_rate, physical_error_rate, color='grey')
    ax.set_xscale('log')
    ax.set_yscale('log')
    plt.show()