Por que o GHCi fornece respostas incorretas abaixo?
GHCi
λ> ((-20.24373193905347)^12)^2 - ((-20.24373193905347)^24)
4.503599627370496e15
Python3
>>> ((-20.24373193905347)**12)**2 - ((-20.24373193905347)**24)
0.0
ATUALIZAÇÃO Eu implementaria a função de Haskell (^) da seguinte maneira.
powerXY :: Double -> Int -> Double
powerXY x 0 = 1
powerXY x y
| y < 0 = powerXY (1/x) (-y)
| otherwise =
let z = powerXY x (y `div` 2)
in if odd y then z*z*x else z*z
main = do
let x = -20.24373193905347
print $ powerXY (powerXY x 12) 2 - powerXY x 24 -- 0
print $ ((x^12)^2) - (x ^ 24) -- 4.503599627370496e15
Embora minha versão não pareça mais correta do que a fornecida abaixo por @WillemVanOnsem, estranhamente fornece a resposta correta para este caso em particular, pelo menos.
Python é semelhante.
def pw(x, y):
if y < 0:
return pw(1/x, -y)
if y == 0:
return 1
z = pw(x, y//2)
if y % 2 == 1:
return z*z*x
else:
return z*z
# prints 0.0
print(pw(pw(-20.24373193905347, 12), 2) - pw(-20.24373193905347, 24))
2.243746917640863e31 - 2.2437469176408626e31
que possui um pequeno erro de arredondamento que é amplificado. Parece um problema de cancelamento.
a^24
é aproximadamente2.2437e31
e, portanto, há um erro de arredondamento que produz isso.