def test_simpson(self):
"""simpson should match hand-calculated values"""
c = array([1,0,2,5,2])
self.assertFloatEqual(simpson(c), .66)
d = array([5])
self.assertFloatEqual(simpson(d), 0)
def test_simpson_e(self):
"""simpson e should match hand-calculated value"""
c = array([1,2,3,1])
s = simpson(c)
self.assertEqual((1/s)/4, simpson_e(c))
