diff options
Diffstat (limited to 'ext/bigdecimal')
-rw-r--r-- | ext/bigdecimal/lib/bigdecimal/nlsolve.rb | 38 | ||||
-rw-r--r-- | ext/bigdecimal/sample/linear.rb | 2 | ||||
-rw-r--r-- | ext/bigdecimal/sample/nlsolve.rb | 2 |
3 files changed, 2 insertions, 40 deletions
diff --git a/ext/bigdecimal/lib/bigdecimal/nlsolve.rb b/ext/bigdecimal/lib/bigdecimal/nlsolve.rb deleted file mode 100644 index 7f729e6aaa..0000000000 --- a/ext/bigdecimal/lib/bigdecimal/nlsolve.rb +++ /dev/null @@ -1,38 +0,0 @@ -#!/usr/local/bin/ruby - -# -# nlsolve.rb -# An example for solving nonlinear algebraic equation system. -# - -require "bigdecimal" -require "bigdecimal/newton" -include Newton - -class Function - def initialize() - @zero = BigDecimal::new("0.0") - @one = BigDecimal::new("1.0") - @two = BigDecimal::new("2.0") - @ten = BigDecimal::new("10.0") - @eps = BigDecimal::new("1.0e-16") - end - def zero;@zero;end - def one ;@one ;end - def two ;@two ;end - def ten ;@ten ;end - def eps ;@eps ;end - def values(x) # <= defines functions solved - f = [] - f1 = x[0]*x[0] + x[1]*x[1] - @two # f1 = x**2 + y**2 - 2 => 0 - f2 = x[0] - x[1] # f2 = x - y => 0 - f <<= f1 - f <<= f2 - f - end -end - f = BigDecimal::limit(100) - f = Function.new - x = [f.zero,f.zero] # Initial values - n = nlsolve(f,x) - p x diff --git a/ext/bigdecimal/sample/linear.rb b/ext/bigdecimal/sample/linear.rb index 7c8ca263eb..88a62ffa71 100644 --- a/ext/bigdecimal/sample/linear.rb +++ b/ext/bigdecimal/sample/linear.rb @@ -11,7 +11,7 @@ # require "bigdecimal" -require "ludcmp" +require "bigdecimal/ludcmp" # # NOTE: diff --git a/ext/bigdecimal/sample/nlsolve.rb b/ext/bigdecimal/sample/nlsolve.rb index 08f17f9ecd..7f729e6aaa 100644 --- a/ext/bigdecimal/sample/nlsolve.rb +++ b/ext/bigdecimal/sample/nlsolve.rb @@ -6,7 +6,7 @@ # require "bigdecimal" -require "newton" +require "bigdecimal/newton" include Newton class Function |