code: def from_diff_fit(cls, strains, stresses, eq_stress=None,
                        tol=1e-10, order=3):
          
          c_list = diff_fit(strains, stresses, eq_stress, order, tol)
          return cls(c_list)
docs: Generates an elastic tensor expansion via the fitting function
          defined below in diff_fit