mirror of
https://github.com/ruby/ruby.git
synced 2022-11-09 12:17:21 -05:00
c4fdfabcc8
g -L frozen_string_literal ext/**/*.rb|xargs ruby -Ka -e'ARGV.each{|fn|puts
fn;open(fn,"r+"){|f|s=f.read.sub(/\A(#!.*\n)?(#.*coding.*\n)?/,"\\&#
frozen_string_literal: false\n");f.rewind;f.write s}}'
git-svn-id: svn+ssh://ci.ruby-lang.org/ruby/trunk@53143 b2dd03c8-39d4-4d8f-98ff-823fe69b080e
40 lines
871 B
Ruby
40 lines
871 B
Ruby
#!/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
|