X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FIntegration%2FTrapezoid.hs;h=df4da781b6cea30e7fbe763ee20da9ba85549db7;hb=7638f1a379035c045185e25cacf0ea2d9d257265;hp=444d92d2eb9a0ed2e6911970f54cd88aafd71b3a;hpb=2f54e89d36e835c58efcc281741632d457859b20;p=numerical-analysis.git diff --git a/src/Integration/Trapezoid.hs b/src/Integration/Trapezoid.hs index 444d92d..df4da78 100644 --- a/src/Integration/Trapezoid.hs +++ b/src/Integration/Trapezoid.hs @@ -1,15 +1,19 @@ +{-# LANGUAGE NoImplicitPrelude #-} {-# LANGUAGE RebindableSyntax #-} -module Integration.Trapezoid +module Integration.Trapezoid ( + trapezoid, + trapezoid_1 ) where -import Misc (partition) +import Misc ( partition ) + +import NumericPrelude hiding ( abs ) +import qualified Algebra.Field as Field ( C ) +import qualified Algebra.RealField as RealField ( C ) +import qualified Algebra.ToInteger as ToInteger ( C ) +import qualified Algebra.ToRational as ToRational ( C ) -import NumericPrelude hiding (abs) -import qualified Algebra.Field as Field -import qualified Algebra.RealField as RealField -import qualified Algebra.ToInteger as ToInteger -import qualified Algebra.ToRational as ToRational -- | Use the trapezoid rule to numerically integrate @f@ over the -- interval [@a@, @b@]. @@ -38,8 +42,9 @@ trapezoid_1 :: (Field.C a, ToRational.C a, Field.C b) -> a -- ^ The \"right\" endpoint, @b@ -> b trapezoid_1 f a b = - (((f a) + (f b)) / 2) * (fromRational' $ toRational (b - a)) - + (((f a) + (f b)) / 2) * coerced_interval_length + where + coerced_interval_length = fromRational' $ toRational (b - a) -- | Use the composite trapezoid rule to numerically integrate @f@ -- over @n@ subintervals of [@a@, @b@].