C++で楕円体の体積を計算するプログラムの作成方法

3つの半径 r1、r2、r3 が与えられたとき、楕円体（だえんたい）の体積を求めるのが本記事の目的です。

楕円体とは、2次曲面の一種であり、「3つの変数を含む2次多項式のゼロ集合として定義できる曲面」として表されます。2次曲面の中でも、楕円体は以下のいずれかの性質によって特徴づけられます。

• 中心を通る任意の平面による断面が楕円（または円）になる
• 3つの互いに垂直な軸方向に異なる半径を持つ球面を伸縮した形状である

楕円体の体積を求める公式

楕円体の体積は、次の公式を使って計算できます。

楕円体の体積 = (4/3) × π × r1 × r2 × r3

これは球の体積公式 (4/3)πr³ を一般化したもので、r1 = r2 = r3 の場合には通常の球の体積と一致します。

入出力例

入力：r1 = 6.3, r2 = 43.4, r3 = 3.7
出力：楕円体の体積は : 4224.87

アルゴリズム

開始
ステップ1 → マクロを定義する
    #define pi 3.14
ステップ2 → 楕円体の体積を計算する関数を宣言する
    float volume(float r1, float r2, float r3)
        return 1.33 * pi * r1 * r2 * r3
ステップ3 → main() 内で処理を行う
    変数を宣言する：float r1 = 6.3, r2 = 43.4, r3 = 3.7
    volume(r1, r2, r3) を呼び出す
終了

C++での実装例

#include <bits/stdc++.h>
#define pi 3.14
using namespace std;

// 楕円体の体積を求める関数
float volume(float r1, float r2, float r3){
    return 1.33 * pi * r1 * r2 * r3;
}

int main(){
    float r1 = 6.3, r2 = 43.4, r3 = 3.7;
    cout << "楕円体の体積は : " << volume(r1, r2, r3);
    return 0;
}

実行結果

楕円体の体積は : 4224.87

補足：より高精度な計算について

上記のコードでは簡略化のため π を 3.14、4/3 を 1.33 としていますが、より正確な結果が必要な場合は、M_PI（<cmath> ヘッダで定義）や分数表現の (4.0 / 3.0) を使用することをおすすめします。例えば次のように書けます。

float volume(float r1, float r2, float r3){
    return (4.0f / 3.0f) * M_PI * r1 * r2 * r3;
}

このように、楕円体の体積計算は非常にシンプルな公式で実装でき、3Dグラフィックスや物理シミュレーションなど、さまざまな分野で応用されています。