src/Integration/Trapezoid.hs: fix monomorphism restriction warning.

[numerical-analysis.git] / src / Integration / Trapezoid.hs

diff --git a/src/Integration/Trapezoid.hs b/src/Integration/Trapezoid.hs
index 444d92d2eb9a0ed2e6911970f54cd88aafd71b3a..807bbc5839cb70e303dbc47d1bbf986305780b19 100644 (file)
--- a/src/Integration/Trapezoid.hs
+++ b/src/Integration/Trapezoid.hs
@@ -1,15 +1,20 @@
+{-# LANGUAGE NoImplicitPrelude #-}
+{-# LANGUAGE ScopedTypeVariables #-}
 {-# 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@].
@@ -32,14 +37,15 @@ import qualified Algebra.ToRational as ToRational
 --   >>> trapezoid_1 f (-1) 1
 --   2.0
 --
-trapezoid_1 :: (Field.C a, ToRational.C a, Field.C b)
+trapezoid_1 :: forall a b. (Field.C a, ToRational.C a, Field.C b)
             => (a -> b) -- ^ The function @f@
             -> a       -- ^ The \"left\" endpoint, @a@
             -> a       -- ^ The \"right\" endpoint, @b@
             -> b
-trapezoid_1 f a b =
-  (((f a) + (f b)) / 2) * (fromRational' $ toRational (b - a))
-
+trapezoid_1 f x y =
+  (((f x) + (f y)) / 2) * coerced_interval_length
+  where
+    coerced_interval_length = fromRational' $ toRational (y - x) :: b
 
 -- | Use the composite trapezoid rule to numerically integrate @f@
 --   over @n@ subintervals of [@a@, @b@].