summaryrefslogtreecommitdiff log msg author committer range
path: root/ext/bigdecimal/sample/nlsolve.rb
blob: b1dd08e0a370801224bf7d882a644b94f1633fa6 (plain)
 ```1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 ``` ``````#!/usr/local/bin/ruby # frozen_string_literal: false # # nlsolve.rb # An example for solving nonlinear algebraic equation system. # require "bigdecimal" require "bigdecimal/newton" include Newton class Function # :nodoc: all 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 ``````