ホイン法による数値計算
yi = 1.0 h = 0.001 a = 1 func = Proc.new{|x| 1.0 + a * x } yi_1 = Num4DiffLib::heunMethod(yi, 0.0, h, func)