mathpack 0.3.2 → 0.4.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 4ff99464e0f1a06b1ba3c2d3423096577a3649d1
-  data.tar.gz: e41163fd0942911d66072c0d1b579d83ba9cec75
+  metadata.gz: 78c7a6ef1eb501cbcc71788e0afe54a97c18b466
+  data.tar.gz: 52b242cd86918bbf029762f6191b86674af05b2f
 SHA512:
-  metadata.gz: 73f068a79ab03282563afd02416d48e1d54fb70b1f88d8286e0ac918215070ee8f3cc674cb731dc12878073b5fd4b0f1bce9ada7818a4271fb9748209c0c7880
-  data.tar.gz: 8287a6d59751b79d6e6ed2df1cebb1fcea04763618f104c838bd747e01c8e835fb657b47e815a246007e1988cf288c87c63690d4e3da6ff32657ee4149e4b75b
+  metadata.gz: 6c25535090593e3f413123b78a63c40140eaa0ef58094059e4b2c321d65d143d320e5d7bedce7e75b3c55a86b4c62f0dc6d9a215e3e75993fa41a15fb12781ee
+  data.tar.gz: 85627190173a1ef07beec6be3c7ceb559c19867e2fb3edbd293c9e7a3bdd47ffc5dd0cfa87f42ab01e0547f24846c694c9fa3ac166f8a6af8246041d5faddbfa

data/lib/mathpack/integration.rb

@@ -0,0 +1,60 @@
+module Mathpack
+  module Integration
+    STEP = 0.0625
+    INTEGRATION_LIMIT = 1e3
+
+    def self.integrate(params = {}, &f)
+      if params[:from] == -Float::INFINITY && params[:to] == Float::INFINITY
+        result = solve_ni_2(&f)
+      elsif params[:to] == Float::INFINITY || params[:from] == -Float::INFINITY
+        if params[:from] == -Float::INFINITY
+          params[:from] = -params[:to]
+          g = ->(x) { f.call(-x) }
+        else
+          g = f
+        end
+        result = solve_ni_1(params[:from], &g)
+      else
+        result = solve_oi(params[:from], params[:to], &f)
+      end
+      result
+    end
+
+    def self.basic_integrate(a, b, &f)
+      result = 5.0/9.0 * f.call((a + b)/2.0 + (b - a)/2.0 * (-(3.0/5.0)**0.5))
+      result += 8.0/9.0 * f.call((a + b)/2.0)
+      result += 5.0/9.0 * f.call((a + b)/2.0 + (b - a)/2.0 * (3.0/5.0)**0.5)
+      result * (b - a)/2.0
+    end
+
+    def self.solve_oi(from, to, &f)
+      result = 0.0
+      nodes = Mathpack::Approximation.generate_nodes(start: from, end: to, step: STEP)
+      for i in 0...nodes.length - 1  do
+        result += basic_integrate(nodes[i], nodes[i + 1], &f)
+      end
+      result
+    end
+
+    def self.solve_ni_1(from, &f)
+      to = [from, INTEGRATION_LIMIT].max
+      result = solve_oi(from, to, &f)
+      coefficients_array = [0.4589646740, 0.4170008308, 0.1133733821, 0.0103991974, 0.0002610172, 0.0000008985]
+      nodes_array = [0.2228466042, 1.8889321017, 2.9927363261, 5.7751435691, 9.8374674184, 15.9828739806]
+      nodes_array.each_index do |i|
+        result += coefficients_array[i] * Math.exp(nodes_array[i]) * f.call(to + nodes_array[i])
+      end
+      result
+    end
+
+    def self.solve_ni_2(&f)
+      result = 0.0
+      coefficients_array = [0.004530009906, 0.1570673203, 0.7246295952, 0.7246295952, 0.1570673203, 0.004530009906]
+      nodes_array = [-2.3506049737, -1.3358490740, -0.4360774119, 0.4360774119, 1.3358490740, 2.3506049737]
+      nodes_array.each_index do |i|
+        result += coefficients_array[i] * Math.exp(nodes_array[i]**2) * f.call(nodes_array[i])
+      end
+      result
+    end
+  end
+end

data/lib/mathpack/version.rb

@@ -1,3 +1,3 @@
 module Mathpack
-  VERSION = '0.3.2'
+  VERSION = '0.4.0'
 end

data/lib/mathpack.rb

@@ -3,3 +3,4 @@ require 'mathpack/statistics'
 require 'mathpack/equation'
 require 'mathpack/sle'
 require 'mathpack/approximation'
+require 'mathpack/integration'

data/spec/integration_spec.rb

@@ -0,0 +1,42 @@
+describe 'Integration' do
+  require 'mathpack/integration'
+
+  context '#a-b integrals' do
+    it 'should integrate e^(-x^2)' do
+      expect(Mathpack::Integration.integrate(from: 0, to: 10){ |x| Math.exp(-x**2) }).to be_between(0.8862269254525, 0.8862269254534)
+    end
+
+    it 'should integrate sin(x^2)' do
+      expect(Mathpack::Integration.integrate(from: 1, to: 7.3){ |x| Math.sin(x**2) }).to be_between(0.384315, 0.384324)
+    end
+
+    it 'should integrate sin(x)/x' do
+      expect(Mathpack::Integration.integrate(from: 0, to: 3.6){ |x| Math.sin(x) / x }).to be_between(1.821945, 1.821954)
+      expect(Mathpack::Integration.integrate(from: -1, to: 1){ |x| Math.sin(x) / x }).to be_between(1.89216614073435, 1.89216614073444)
+    end
+
+    it 'should integrate 1/lnx' do
+      expect(Mathpack::Integration.integrate(from: 1.1, to: 4.7){ |x| 1 / Math.log(x) }).to be_between(5.120315, 5.120324)
+    end
+  end
+
+  context '#ni-1 integrals' do
+    it 'should integrate e^(-x)/(x+1)' do
+      expect(Mathpack::Integration.integrate(from: 0, to: Float::INFINITY){ |x| Math.exp(-x) / (x + 1) }).to be_between(0.596346, 0.596348)
+    end
+
+    it 'should count gamma function' do
+      expect(Mathpack::Integration.integrate(from: 0, to: Float::INFINITY){ |x| Math.exp(-x) * x**1.5 }).to be_between(1.329339, 1.329341)
+    end
+
+    it 'should count Laplace function' do
+      expect(Mathpack::Integration.integrate(from: -Float::INFINITY, to: 1.65){ |x| 1.0 / (2 * Math::PI)**0.5 * Math.exp(-x**2/2) }).to be_between(0.950450, 0.950549)
+    end
+  end
+
+  context '#ni-2 integrals' do
+    it 'should integrate e^(-x^2)*cos(x)' do
+      expect(Mathpack::Integration.integrate(from: -Float::INFINITY, to: Float::INFINITY){ |x| Math.exp(-x**2) * Math.cos(x) }).to be_between(1.3803, 1.3804)
+    end
+  end
+end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: mathpack
 version: !ruby/object:Gem::Version
-  version: 0.3.2
+  version: 0.4.0
 platform: ruby
 authors:
 - maxmilan
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2015-02-12 00:00:00.000000000 Z
+date: 2015-02-21 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -54,12 +54,14 @@ files:
 - lib/mathpack.rb
 - lib/mathpack/approximation.rb
 - lib/mathpack/equation.rb
+- lib/mathpack/integration.rb
 - lib/mathpack/sle.rb
 - lib/mathpack/statistics.rb
 - lib/mathpack/version.rb
 - mathpack.gemspec
 - spec/approximation_spec.rb
 - spec/equation_spec.rb
+- spec/integration_spec.rb
 - spec/sle_spec.rb
 - spec/statistics_spec.rb
 homepage: ''