第37講 円周率を計算する△
第4話 
オイラーの公式による円周率の計算コード例
オイラーの公式
πの2乗/6=1/1*1+1/2*2+1/3*3+1/4*4+・・・
コード例
#pragma once
#include<math.h>
namespace スレッドの起動 {
       ・
       ・
       ・
#pragma endregion
  private: System::Void button1_Click(System::Object^ sender, System::EventArgs^ e) {
           int k;
           long j=1;
           double i,pi;
           double h=-1;
           pi=0;
           k=long::Parse(textBox1->Text);
           if(k>1){
             for(i=0;i<k-1;i++){
               j*=10;
             }
           }
           for(i=1;i<j+1;i++){
             pi+=1/(i*i);
           }
           pi=sqrt(6*pi);
           textBox2->Text=pi.ToString();
        }
};
}
実行例
オイラー
計算速度はあまり変わりませんが、精度は少し上がりました。

π = 3.14159 26535 89793 23846 26433 83279 50288 …


そこで、確かこれもオイラーが発見した公式だと思いますが、
πの4乗/90=1/1の4乗+1/2の4乗+1/3の4乗+1/4の4乗+・・・
を利用したらどうでしょうか。



第3話へ 第5話へ 第38講第1話へ

戻る

VC++講義第1部へ
vb講義へ
VB講義基礎へ
初心者のための世界で一番わかりやすいVisual C++入門基礎講座
初心者のための世界で一番わかりやすいVisual Basic入門基礎講座
初心者のための世界で一番わかりやすいVBA入門講義(基礎から応用まで)