用 Ruby 求定积分

最新推荐文章于 2024-06-09 06:30:00 发布

aabbcc456aa

最新推荐文章于 2024-06-09 06:30:00 发布

阅读量718

点赞数

分类专栏：其他算法（algorithm）

其他同时被 2 个专栏收录

27 篇文章

订阅专栏

算法（algorithm）

6 篇文章

订阅专栏

galeki posted @ 2008年6月16日 14:23 in RoR , 2220 阅读

纯属无聊……

曾经以为求积分之类的是个很复杂的过程，对那些可以求出积分值的计算程序佩服不以，昨日脑子不知道为什么忽然想到这个问题，翻了翻书，发现这个问题很简单～

用最简单的矩形法，Ruby 代码如下：

 
       def integral
       (a, b, n = 
       100
       )
      
           sum = 
       0.0
      
           dx = 
       (b - a
       )/n.
       to_f
      
       n.times 
       do
      
               a += dx
      
               sum = sum + 
       (
       yield a
       )*dx
      
       end
      
           sum
      
       end

n 为精度，Ruby 的 block 真方便，如果要求：

$\int_0^1 (x^2-x^3) \mathrm{d}x$

那么就是：

如果要求：

$\int_0^\pi sin(x) \mathrm{d}x$

那么就是：

 
    
       integral
       (
       0, 
       Math::PI
       ) 
       {|x| 
       Math.
       sin
       (x
       )
       }
      
 
   

如果函数很复杂，也可以在外面定义一个 f(x)，然后写成 integral(a, b) {|x| f(x)}，相当美观～

看看结果：

 
    
       puts integral
       (
       0, 
       Math::PI
       ) 
       {|x| 
       Math.
       sin
       (x
       )
       }
      

       >> 
       1.99983550388745
      
 
   