it seems that when plotting steady state P wrt A,eps for the two notebooks below gives different results when using the same paramters (paramters used:
delta = 1.0 * 2 * np.pi # qubit sigma_x coefficient
f = 5
w = f * 2 * np.pi # driving frequency
omega=w
# T = 2 * np.pi / w # driving period
gamma1 = 10/w # relaxation rate
gamma2 = 0.1/w # dephasing rate
# omega = 0.05 * 2 * np.pi # qubit sigma_z coefficient
A_vec = np.linspace(0.0, 20.0, 50) * 2 * np.pi # driving amplitude
eps0_vec = np.linspace(-20.0, 20.0, 50) * 2 * np.pi
results from notebook 1)