Sha256: 3ed65975752ac771a60b3ace4bc670d2623ac2a7f1503e84868e7d7fb3e2a8f9
Contents?: true
Size: 613 Bytes
Versions: 22
Compression:
Stored size: 613 Bytes
Contents
#!/usr/bin/env ruby require("gsl") if ARGV.size != 1 puts("Usage: integrate n\nn is the number of quadrature points.") exit end Q = ARGV[0].to_i if Q < 1 puts("Usage: integrate n\nn is the number of quadrature points.") exit end def fun(x) Math::cos(3.0*x) end def dfun(x) -3.0*Math::sin(3.0*x) end f = GSL::Vector.alloc(4*Q) quad = Jac::Quadrature.alloc(Q) quad.zwd(Jac::GJ, 0.0, 0.0) for i in 0...Q do f[i] = fun(quad.x[i]) end d = quad.differentiate(f) printf("X \t Derivative \t Error\n"); for i in 0...Q do printf("%f \t %f \t %e\n", quad.x[i], d[i], d[i] - dfun(quad.x[i])); end
Version data entries
22 entries across 22 versions & 4 rubygems