code: def get_extrapolated_diffusivity(temps, diffusivities, new_temp):
      
      Ea, c, _ = fit_arrhenius(temps, diffusivities)
      return c * np.exp(-Ea / (const.k / const.e * new_temp))
docs: Returns (Arrhenius) extrapolated diffusivity at new_temp
  
      Args:
          temps ([float]): A sequence of temperatures. units: K
          diffusivities ([float]): A sequence of diffusivities (e.g.,
              from DiffusionAnalyzer.diffusivity). units: cm^2/s
          new_temp (float): desired temperature. units: K
  
      Returns:
          (float) Diffusivity at extrapolated temp in mS/cm.