def test_powerlaw_stats():
 cases = [(1.0, (0.5, (1.0 / 12), 0.0, (-1.2))), (2.0, ((2.0 / 3), (2.0 / 36), (-0.5656854249492473), (-0.6)))]
 for (a, exact_mvsk) in cases:
  mvsk = stats.powerlaw.stats(a, moments='mvsk')
  assert_array_almost_equal(mvsk, exact_mvsk)

def test_powerlaw_stats ( ) : cases = [ ( 1.0 , ( 0.5 , ( 1.0 / 12 ) , 0.0 , ( - 1.2 ) ) ) , ( 2.0 , ( ( 2.0 / 3 ) , ( 2.0 / 36 ) , ( - 0.5656854249492473 ) , ( - 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 .