Sha256: 957491e768e2d011bb5795e7b951f64f93eb664daccb4cdda27a11f6658bd1ad

Contents?: true

Size: 390 Bytes

Versions: 22

Compression:

Stored size: 390 Bytes

Contents

#!/usr/bin/env ruby
require("gsl")

N = 100

y0 = 1.0
A = 2.0
x0 = 3.0
B = 4.0

r = GSL::Rng.alloc()
x = GSL::Vector.linspace(-10, 20, N)
y =  y0 + A/(GSL::pow_2(x-x0)+B) + 0.03*GSL::Ran::gaussian(r, 1, N)

coef, err, chi2, dof = GSL::MultiFit::FdfSolver.fit(x, y, "lorentzian")
y0 = coef[0]
amp = coef[1]
x0 = coef[2]
b = coef[3]
p coef
p err
GSL::graph(x, y, y0+amp/(GSL::pow_2(x-x0)+b))

Version data entries

22 entries across 22 versions & 4 rubygems

Version Path
gsl-2.1.0.3 examples/fit/lorentzfit.rb
gsl-2.1.0.2 examples/fit/lorentzfit.rb
gsl-2.1.0.1 examples/fit/lorentzfit.rb
gsl-2.1.0 examples/fit/lorentzfit.rb
gsl-1.16.0.6 examples/fit/lorentzfit.rb
rb-gsl-1.16.0.5 examples/fit/lorentzfit.rb
rb-gsl-1.16.0.4 examples/fit/lorentzfit.rb
rb-gsl-1.16.0.3 examples/fit/lorentzfit.rb
rb-gsl-1.16.0.3.rc1 examples/fit/lorentzfit.rb
gsl-nmatrix-1.15.3.2 examples/fit/lorentzfit.rb
gsl-nmatrix-1.15.3.1 examples/fit/lorentzfit.rb
rb-gsl-1.16.0.2 examples/fit/lorentzfit.rb
rb-gsl-1.16.0.1 examples/fit/lorentzfit.rb
rb-gsl-1.16.0 examples/fit/lorentzfit.rb
rb-gsl-1.15.3.2 examples/fit/lorentzfit.rb
rb-gsl-1.15.3.1 examples/fit/lorentzfit.rb
blackwinter-gsl-1.15.3.2 examples/fit/lorentzfit.rb
gsl-nmatrix-1.15.3.0 examples/fit/lorentzfit.rb
gsl-1.15.3 examples/fit/lorentzfit.rb
gsl-1.14.7 examples/fit/lorentzfit.rb