def test_powerlaw_stats():
    cases = [(1.0, (0.5, 1./12, 0.0, -1.2)),
             (2.0, (2./3, 2./36, -0.56568542494924734, -0.6))]
    for a, exact_mvsk in cases:
        mvsk = stats.powerlaw.stats(a, moments="mvsk")
        assert_array_almost_equal(mvsk, exact_mvsk)
Test the powerlaw stats function.
dataset/ETHPy150Open scipy/scipy/scipy/stats/tests/test_distributions.py/test_powerlaw_stats