code: def gradient(self, mu, dist):
          
          return dist.levels/(mu*(dist.levels - mu))
docs: derivative of the link function wrt mu
  
          Parameters
          ----------
          mu : array-like of legth n
          dist : Distribution instance
  
          Returns
          -------
          grad : np.array of length n