def METHOD_NAME(self, stddev, axes):
    low_slip, high_slip = stddev
    axes.fill_between(self.impulse_coords[:,1] / 1.0e+3, low_slip[:,1], high_slip[:,1], color=self.COLOR_MODEL, lw=0, alpha=0.5, label="median$\pm$stddev")