X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FIntegration%2FSimpson.hs;h=c9ad414b4fa7cc85705af2604d2295524ceeb01e;hb=2f54e89d36e835c58efcc281741632d457859b20;hp=c3d59ff797c0538c377a08196065a63b43178003;hpb=2d0ecad8695e129443e53311d6494b2465f1a672;p=numerical-analysis.git diff --git a/src/Integration/Simpson.hs b/src/Integration/Simpson.hs index c3d59ff..c9ad414 100644 --- a/src/Integration/Simpson.hs +++ b/src/Integration/Simpson.hs @@ -6,10 +6,7 @@ where import Misc (partition) import NumericPrelude hiding (abs) -import Algebra.Absolute (abs) -import qualified Algebra.Field as Field import qualified Algebra.RealField as RealField -import qualified Algebra.RealRing as RealRing import qualified Algebra.ToInteger as ToInteger import qualified Algebra.ToRational as ToRational @@ -26,6 +23,7 @@ import qualified Algebra.ToRational as ToRational -- >>> simpson_1 f (-1) 1 -- 0.0 -- +-- >>> import Algebra.Absolute (abs) -- >>> let f x = x^2 -- >>> let area = simpson_1 f (-1) 1 -- >>> abs (area - (2/3)) < 1/10^12 @@ -56,6 +54,7 @@ simpson_1 f a b = -- -- Examples: -- +-- >>> import Algebra.Absolute (abs) -- >>> let f x = x^4 -- >>> let area = simpson 10 f (-1) 1 -- >>> abs (area - (2/5)) < 0.0001 @@ -64,6 +63,7 @@ simpson_1 f a b = -- Note that the convergence here is much faster than the Trapezoid -- rule! -- +-- >>> import Algebra.Absolute (abs) -- >>> let area = simpson 10 sin 0 pi -- >>> abs (area - 2) < 0.00001 -- True