code: def plot_s_bcr(fignum, Bcr, S, sym):
      
      plt.figure(num=fignum)
      plt.plot(Bcr, S, sym)
      plt.xlabel('Bcr')
      plt.ylabel('Mr/Ms')
      plt.title('Squareness-Bcr Plot')
      bounds = plt.axis()
      plt.axis([0, bounds[1], 0, 1])
docs: function to plot Squareness,Coercivity of remanence
  
      Parameters
      __________
      fignum : matplotlib figure number
      Bcr : list or array coercivity of remenence values
      S : list or array of ratio of saturation remanence to saturation
      sym : matplotlib symbol (e.g., 'g^' for green triangles)