Current Path : /compat/linux/proc/self/root/usr/local/share/examples/ruby18/bigdecimal/ |
FreeBSD hs32.drive.ne.jp 9.1-RELEASE FreeBSD 9.1-RELEASE #1: Wed Jan 14 12:18:08 JST 2015 root@hs32.drive.ne.jp:/sys/amd64/compile/hs32 amd64 |
Current File : //compat/linux/proc/self/root/usr/local/share/examples/ruby18/bigdecimal/nlsolve.rb |
#!/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