]> gitweb.michael.orlitzky.com - numerical-analysis.git/commitdiff
src/Integration/Simpson.hs: fix monomorphism restriction warnings.
authorMichael Orlitzky <michael@orlitzky.com>
Fri, 7 Dec 2018 16:04:36 +0000 (11:04 -0500)
committerMichael Orlitzky <michael@orlitzky.com>
Fri, 7 Dec 2018 16:04:36 +0000 (11:04 -0500)
src/Integration/Simpson.hs

index f0f57f29245db1a7db7bd10bd05a39a87c2347df..afa932b672cee04075a43c52ccacba1e9f56be2c 100644 (file)
@@ -1,4 +1,5 @@
 {-# LANGUAGE NoImplicitPrelude #-}
+{-# LANGUAGE ScopedTypeVariables #-}
 {-# LANGUAGE RebindableSyntax #-}
 
 module Integration.Simpson (
@@ -40,16 +41,16 @@ import qualified Algebra.ToRational as ToRational ( C )
 --   >>> simpson_1 f 0 1
 --   0.25
 --
-simpson_1 :: (RealField.C a, ToRational.C a, RealField.C b)
+simpson_1 :: forall a b. (RealField.C a, ToRational.C a, RealField.C b)
             => (a -> b) -- ^ The function @f@
             -> a       -- ^ The \"left\" endpoint, @a@
             -> a       -- ^ The \"right\" endpoint, @b@
             -> b
-simpson_1 f a b =
-  coefficient * ((f a) + 4*(f midpoint) + (f b))
+simpson_1 f x y =
+  coefficient * ((f x) + 4*(f midpoint) + (f y))
   where
-    coefficient = fromRational' $ (toRational (b - a)) / 6
-    midpoint = (a + b) / 2
+    coefficient = fromRational' $ (toRational (y - x)) / 6 :: b
+    midpoint = (x + y) / 2
 
 
 -- | Use the composite Simpson's rule to numerically integrate @f@