code: def plot_intrusion_curve(self, fig=None, inv_points=None):
          r
          data = self.get_intrusion_data(inv_points)
          if fig is None:
              fig = plt.figure()
          a = fig.add_subplot(111)
          a.plot(data.Pcap, data.S_pore, 'r*-', label='pore')
          a.plot(data.Pcap, data.S_throat, 'b*-', label='throat')
          a.plot(data.Pcap, data.S_tot, 'g*-', label='total')
          a.legend(bbox_to_anchor=(0, 1.02, 1, 0.102),
                   loc=3, ncol=3, borderaxespad=0)
          a.set_xlabel("Capillary Pressure [Pa]")
          a.set_ylabel("Saturation")
          a.set_ybound(lower=0.0, upper=1.0)
          return fig
docs: r"""
          Plot a simple drainage curve