Sha256: d13900667589a91be2b41216db3491f5a34994f8ec8f81eec7ec95c836957a8d

Contents?: true

Size: 738 Bytes

Versions: 16

Compression:

Stored size: 738 Bytes

Contents

#!/usr/bin/env ruby
require("gsl")
include GSL::Sf

l = 2
x = 0.01
File.open("sphbessel.dat", "w") do |file|
  while x <= 12.0
    j0 = bessel_j0(x)
    j1 = bessel_j1(x)
    j2 = bessel_jl(2, x)
    j3 = bessel_jl(3,x)
    y0 = bessel_y0(x)
    y1 = bessel_y1(x)
    y2 = bessel_y2(x)
    i0 = bessel_i0_scaled(x)
    i1 = bessel_i1_scaled(x)
    i2 = bessel_i2_scaled(x)
    k0 = bessel_k0_scaled(x)
    k1 = bessel_k1_scaled(x)
    k2 = bessel_k2_scaled(x)
    file.printf("%e %e %e %e %e %e %e %e %e %e %e %e %e %e\n", 
                x, j0, j1, j2, j3, y0, y1, y2, i0, i1, i2, k0, k1, k2)
    x += 0.1
  end
end

#puts("sphbessel.dat created. Try sphbessel.gp.")
system("gnuplot -persist sphbessel.gp")
File.delete("sphbessel.dat")

Version data entries

16 entries across 16 versions & 4 rubygems

Version Path
rb-gsl-1.16.0.4 examples/sf/sphbessel.rb
rb-gsl-1.16.0.3 examples/sf/sphbessel.rb
rb-gsl-1.16.0.3.rc1 examples/sf/sphbessel.rb
gsl-nmatrix-1.15.3.2 examples/sf/sphbessel.rb
gsl-nmatrix-1.15.3.1 examples/sf/sphbessel.rb
rb-gsl-1.16.0.2 examples/sf/sphbessel.rb
rb-gsl-1.16.0.1 examples/sf/sphbessel.rb
rb-gsl-1.16.0 examples/sf/sphbessel.rb
rb-gsl-1.15.3.2 examples/sf/sphbessel.rb
rb-gsl-1.15.3.1 examples/sf/sphbessel.rb
blackwinter-gsl-1.15.3.2 examples/sf/sphbessel.rb
gsl-nmatrix-1.15.3.0 examples/sf/sphbessel.rb
gsl-1.15.3 examples/sf/sphbessel.rb
gsl-1.14.7 examples/sf/sphbessel.rb
gsl-1.14.6 examples/sf/sphbessel.rb
gsl-1.14.5 examples/sf/sphbessel.rb