haskell中的Heron方法

Question

我被除以零例外，在这个代码中，赫伦的方法，我有点迷失在这里。

epsilon:: Integral a => a
epsilon = 1

heron:: Integral a => a -> a
heron r = help 0
  where
    help x
      | abs (heron' x - heron' (x + 1)) < epsilon = heron' (x + 1)
      | otherwise                                 = help (x + 1)

    heron' 0 = 1
    heron' x = (1 `div` 2) * (heron' (x-1) + (r `div` heron' (x-1)))

有什么建议吗?在这段代码中，我需要寻找哪些方法来解决这个问题？

(1 `div` 2)绝对是个问题，但是我需要写什么呢？

Answer 1

如果您需要这类除法，您可能希望使用(/)代替div，使用Fractional代替Integral。所以：

epsilon:: Fractional a => a
epsilon = 1

heron:: (Fractional a, Ord a) => a -> a
heron r = help 0
  where
    help x
      | abs (heron' x - heron' (x + 1)) < epsilon = heron' (x + 1)
      | otherwise                                 = help (x + 1)

    heron' 0 = 1
    heron' x = (1 / 2) * (heron' (x-1) + (r / heron' (x-1)))