アルゴリズムのチカラ

昨日の社内勉強会で詰まって悔しかったのでリベンジ。 x^3 + x^2 + x = 1 となるx を求めよ（近似値 = 0.000001とする）という三次方程式の解を求めるプログラム。 まず、ループを使ったベタな線形探索。 解答に至るまで543690ステップの計算を要しています。 EPSILON = 0.000001 # 近似値イプシロン FUNC_EQUAL = 1 # 三次方程式func の値 # 三次方程式func def func(x) return x**3 + x**2 + x end x = 0 # 仮定解x i = 0 # 計算回数カウンタ while func(x) <= FUNC_EQUAL do # FUNC_EQUAL に到達するまでぶん回し x += EPSILON i += 1 end puts "#{x}:[#{i}]" # 0.543689999994796