File: qagp.rb

package info (click to toggle)
ruby-gsl 2.1.0.1%2Bdfsg1-1
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 4,892 kB
  • ctags: 5,459
  • sloc: ansic: 61,660; ruby: 15,656; sh: 19; makefile: 10
file content (20 lines) | stat: -rwxr-xr-x 400 bytes parent folder | download | duplicates (10) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 
#!/usr/bin/env ruby
require 'gsl'
include GSL
include Math

f454 = Function.alloc{ |x|
  x2 = x*x
  x3 = x2*x
  x3*log(((x2-1)*(x2-2)).abs)
}
exp_result = 5.274080611672716401E+01
exp_abserr = 1.755703848687062418E-04
pts = [0, 1, sqrt(2), 3]
result = f454.qagp(pts, 0.0, 1e-3)
p result

puts("exp_result: #{exp_result}")
puts("exp_abserr: #{exp_abserr}")

p Integration.qagp(f454, pts, [0.0, 1e-3])