eg1:
eg2:
x1=mux1_mcmc[0]
y1=muy1_mcmc[0]
x2=mux2_mcmc[0]
y2=muy2_mcmc[0]
x1err=sigmax1_mcmc[0]
y1err=sigmay1_mcmc[0]
x2err=sigmax2_mcmc[0]
y2err=sigmay2_mcmc[0]
x=samplers[:,0]
y=samplers[:,1]
print(len(x))
print(np.min(x),np.max(x),np.min(y),np.max(y))
left, width = 0.1, 0.65
bottom, height = 0.1, 0.65
spacing = 0.005
rect_scatter = [left, bottom, width, height]
rect_histx = [left, bottom + height + spacing, width, 0.2]
rect_histy = [left + width + spacing, bottom, 0.2, height]
# start with a rectangular Figure
plt.figure(figsize=(8, 8))
ax_scatter = plt.axes(rect_scatter)
ax_scatter.tick_params(direction='