code: def plot_bcr(fignum, Bcr1, Bcr2):
      
      plt.figure(num=fignum)
      plt.plot(Bcr1, Bcr2, 'ro')
      plt.xlabel('Bcr1')
      plt.ylabel('Bcr2')
      plt.title('Compare coercivity of remanence')
docs: function to plot two estimates of Bcr against each other