code: def test_method_of_moments():
      # example from https://en.wikipedia.org/wiki/Method_of_moments_(statistics)
      np.random.seed(1234)
      x = [0, 0, 0, 0, 1]
      a = 1 / 5 - 2 * np.sqrt(3) / 5
      b = 1 / 5 + 2 * np.sqrt(3) / 5
      # force use of method of moments (uniform.fit is overridden)
      loc, scale = super(type(stats.uniform), stats.uniform).fit(
          x, method="MM"
      )
      npt.assert_almost_equal(loc, a, decimal=4)
      npt.assert_almost_equal(loc + scale, b, decimal=4)
  
imports: import sys
  import numpy
  import numpy.testing
  import pytest
  from pytest import raises
  from scipy.integrate import IntegrationWarning
  import itertools
  from scipy import stats
  from .common_tests import check_normalization, check_moment, check_mean_expect, check_var_expect, check_skew_expect, check_kurt_expect, check_entropy, check_private_entropy, check_entropy_vect_scale, check_edge_support, check_named_args, check_random_state_property, check_meth_dtype, check_ppf_dtype, check_cmplx_deriv, check_pickling, check_rvs_broadcast, check_freezing, check_munp_expect
  from scipy.stats._distr_params import distcont
  from scipy.stats._distn_infrastructure import rv_continuous_frozen
label_desc: Use image node_id 15 for calling a global function with example usage: test_method_of_moments() without return types