正弦関数の再帰計算は役に立たない

Question

私はハスケルを学習しています。私は自分自身でいくつかの関数を実行していますが、この関数には再帰を使った正弦の計算が含まれていますが、

そして、私のコードはこれです：：

私は正弦を計算するために使用している数式は、この一つです

--Returns n to power p 
pow :: Float->Integer->Float 
pow n p = 
    if p == 0 then 
     1 
    else 
     if p == 1 then 
     n 
     else 
     n * (pow n (p-1)) 

--Finds a number's factorial 
f :: Integer->Integer 
f n = 
    if n == 1 then 
     n 
    else 
     n * (f (n-1)) 

--TODO: Trigonometric functions (:v I'll do diz 2) 
sinus :: Float->Char->Float 
sinus n deg = 
    if(deg == 'd')then 
     sinusr 0 (normalize (torad n)) 0 
    else 
     sinusr 0 (normalize n) 0 

--Get the value equivalent to radians of the presented degrees 
torad :: Float->Float 
torad v = ((v * pi)/180) 

--Recursive to get the value of the entering radians 
sinusr :: Integer->Float->Float->Float 
sinusr k x result = 
    if k == 130 then 
    result + (((pow (-1) k) * ((pow x ((2*k)+1)))/(fromIntegral (f ((2*k)+1))))) 
    else 
    result + (sinusr (k+1) x (((pow (-1) k) * ((pow x ((2*k)+1)))/(fromIntegral (f ((2*k)+1)))))) 

--Subtracts pi/2 the necessary times to get a value minor or equals to pi/2 :v 
normalize :: Float->Float 
normalize a = a - (fromIntegral (truncate (a/(pi*2)))*(pi*2)) 

は、例えば、出力は、それがこのです：

*Main> sinus 1 'd' 
1.7452406e-2 
*Main> sinus 1 's' 
0.84147096 
*Main> sinus 2 's' 
NaN 
*Main> sinus 2 'd' 
3.4899496e-2 

誰かが私にそれを示している理由を教えてもらえますか？ 私はLispと同じロジックを使いましたが、完璧に動作しました。私はHaskellの構文を理解しなければなりませんでしたが、わかるように、動作していません。 あらかじめ、ありがとうございます。

Answer 1

三角演算は、三角関数を計算するのに十分ではありません。指数部には、大きな中間の数値のビット数がsinusrで十分ではありません。または、鈍いように、以下の番号はFloatに適合しない：

ghci> 2^130 :: Float 
Infinity 

とすぐに、あなたは通常、それらのかNaNのいずれかで終わる浮動小数点数（-Infinity、Infinity）の境界を打つよう。

代わりにDoubleを使用してください。 lisp の実装ではおそらくは倍精度浮動小数点数も使用します。さらに、すべてのステップで分数全体を再計算しないで、代わりに指名子と分母を更新すると、Floatの値が大きくなりすぎることはありません。