16/01/03 15:15:54.32 gIFHPruM.net
>>255 Haskell
適当にやったけどそれなりの値に収束したぞ
main = print $ f 3 1 !! 10000
f x d = x : f (x - g x * d) (d * 0.9)
g x = sqrt(1-((x^2+10^2-8^2)/(x*10*2))^2)*(x+1)-7
--
7.958400097301163
16/01/03 15:15:54.32 gIFHPruM.net
>>255 Haskell
適当にやったけどそれなりの値に収束したぞ
main = print $ f 3 1 !! 10000
f x d = x : f (x - g x * d) (d * 0.9)
g x = sqrt(1-((x^2+10^2-8^2)/(x*10*2))^2)*(x+1)-7
--
7.958400097301163