mk-mode BLOG

このブログは自作の自宅サーバに構築した Debian GNU/Linux で運用しています。
PC・サーバ構築等の話題を中心に公開しております。(クローンサイト: GitHub Pages

ブログ開設日2009-01-05
サーバ連続稼働時間
Reading...
Page View 合計
Reading...
今日
Reading...
昨日
Reading...

Ruby - 数値積分(台形則による定積分)!

[ プログラミング, 数学 ] [ Ruby ]

こんばんは。

先日は、数値積分(台形則による定積分)の C++ によるアルゴリズムを紹介しました。

今日は、同じアルゴリズムを Ruby で実現してみました。 アルゴリズムについては、上記リンクの記事を参照してください。

実際、ほとんど同じです。

以下、Ruby によるサンプルスクリプトです。

記録

0. 前提条件

  • Cygwin 1.7.15
  • Ruby 1.9.3-p194

1. Ruby スクリプト作成

今回作成した Ruby ソースは以下の通りです。

  • a も b も 0 だったら何もせず終了するようにしている。
  • 関数 f(x) はメソッドにした。
  • 被積分関数は とした。

【 ファイル名: definite_integral_trapzoid.rb 】

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
# -*- coding: utf-8 -*-
#*********************************************
# 台形則による定積分
#*********************************************
class DefiniteIntegralTrapzoid
  # 積分区間分割数
  M  = 100

  # 計算クラス
  class Calc
    # 被積分関数
    def f(x)
      return Math.sqrt(4 - x * x)
    end

    # 定積分計算
    def generate_integral(a, b)
      # 1区間の幅
      h = (b - a) / M

      # 初期化
      x = a  # X 値を a で初期化
      s = 0  # 面積初期化

      # 計算
      1.upto(M - 1) do |k|
        x = x + h
        s = s + f(x)
      end
      s = h * ((f(a) + f(b)) / 2 + s)

      # 結果表示
      printf("  /%f\n",b);
      printf("  |  f(x)dx = %f\n", s);
      printf("  /%f\n",a);
    end
  end

  # メイン処理
  begin
    # データ入力
    print "積分区間 A :"
    a = gets.chomp.to_f
    print "         B :"
    b = gets.chomp.to_f
    exit if a == 0 && b == 0

    # 計算クラスインスタンス化
    obj_calc = Calc.new

    # 定積分計算
    obj_calc.generate_integral(a, b)
  rescue => e
    # エラーメッセージ
    puts "[例外発生] #{e}"
  end
end

2. 実行

実際に実行して検証してみる。

1
2
3
4
5
6
$ ruby definite_ingegral_trapezoid.rb
積分区間 A :0
         B :2
  /2.000000
  |  f(x)dx = 3.140417
  /0.000000

積分記号(インテグラル)をイメージした出力にしています。


当然、C++ 版と同じ結果となりました。

Ruby に限らず、プログラミング言語で色々計算してみるのは、実におもしろいです。

以上。

Comments