---

C言語 math.hとは？主要関数と使い方を網羅

C言語で数値計算や科学技術計算を行う際、欠かせないのが標準ライブラリの数学関数群です。これらの関数は主に <math.h> ヘッダーファイルに定義されており、三角関数、指数関数、対数関数、べき乗、平方根といった基本的な計算から、より高度な特殊関数まで、幅広い機能を提供します。

この記事では、C言語の <math.h> ライブラリに焦点を当て、その概要から主要な関数、使用方法、そして利用する上での注意点までを詳細に解説します。C言語で数値を取り扱うプログラムを開発するすべての方にとって、この記事が <math.h> を理解し、効果的に活用するための手助けとなれば幸いです。

1. はじめに：C言語と数値計算、標準ライブラリの役割

C言語は、その高速性と低レベルな操作能力から、オペレーティングシステムの開発や組み込みシステムなど、ハードウェアに近い分野で広く利用されています。一方で、科学技術計算やデータ処理といった数値計算の分野でも強力なツールとなり得ます。しかし、基本的なC言語の機能だけでは、サインやコサイン、平方根のような複雑な数学的演算を簡単に行うことはできません。

ここで登場するのが「標準ライブラリ」です。標準ライブラリは、C言語の仕様で標準化された、あらかじめ用意されている関数の集まりです。これらの関数を利用することで、ファイル入出力、文字列操作、メモリ管理、そして数学計算など、様々なタスクを効率的に行うことができます。

<math.h> は、C言語の標準ライブラリの一部として、浮動小数点数に対する数学関数を提供するために存在します。高度なアルゴリズムの実装や物理シミュレーション、データ解析など、数値的な処理が必要なあらゆるCプログラムにおいて、<math.h> は不可欠なツールとなります。

この記事では、<math.h> が提供する豊富な関数群の中から、特に利用頻度の高いものを中心に、その機能、使い方、そして関連する概念（浮動小数点数表現、エラー処理など）を深掘りしていきます。

2. math.hとは何か

<math.h>（C++では <cmath> に対応しますが、本記事はC言語に限定します）は、C言語の標準ライブラリヘッダーファイルの一つです。国際標準化機構（ISO）によって定められたC言語の仕様（ANSI C、ISO Cなど）の一部として標準化されています。このヘッダーファイルをインクルードすることで、プログラム内で数学関数を使用できるようになります。

<math.h> に含まれる関数のほとんどは、浮動小数点数を引数に取り、浮動小数点数を返します。対象となる浮動小数点数の型は、C言語の標準で定められている float, double, long double です。

歴史的に、C言語の初期の標準であるC89/C90では、<math.h> の関数の多くは double 型の引数を受け取り、 double 型の戻り値を返すものだけが標準とされていました。しかし、より高い精度や効率が求められるようになり、C99標準では float 型用の関数（関数名の末尾に f がつく、例: sinf）と long double 型用の関数（関数名の末尾に l がつく、例: sinl）が追加されました。これにより、プログラマは使用したい浮動小数点数の型に応じて、適切な関数を選択できるようになりました。さらに、C99では引数の型に応じて自動的に適切な関数を選択する「型総称マクロ」が <tgmath.h> として導入されました（これについては後述します）。

<math.h> を利用する上で重要な概念の一つが、浮動小数点数の表現です。コンピュータは実数を厳密に表現することが難しいため、近似的な表現である浮動小数点数（IEEE 754標準など）を使用します。これにより、計算に伴う「丸め誤差」や、表現できる値の範囲（オーバーフロー、アンダーフロー）、特殊な値（無限大 INF、非数 NaN – Not a Number）といった問題が発生する可能性があります。<math.h> の関数はこれらの特殊な値や誤差の可能性を考慮して設計されており、多くの場合、エラーが発生した際には特定の戻り値を返したり、標準で定められたエラーフラグを設定したりします。

<math.h> に含まれる主な関数群は以下のカテゴリに分けられます。

• 三角関数: サイン、コサイン、タンジェントなど
• 双曲線関数: 双曲線サイン、コサイン、タンジェントなど
• 指数関数と対数関数: 自然対数、常用対数、指数関数など
• べき乗と平方根: べき乗、平方根、立方根など
• 絶対値: 浮動小数点数の絶対値
• 丸め関数: 切り上げ、切り捨て、四捨五入など
• 剰余関数: 浮動小数点数の剰余
• 浮動小数点操作関数: 浮動小数点数を構成要素に分解したり、操作したりする関数
• 分類・比較関数: 浮動小数点数の状態（有限、無限、非数など）を判定したり、安全に比較したりする関数
• 特殊関数: 誤差関数、ガンマ関数など（C99以降）

これらの関数群を、それぞれのカテゴリごとに詳しく見ていきましょう。

3. 主要関数群の詳細な解説と使用例

3.1. 三角関数 (Trigonometric Functions)

角度に対するサイン、コサイン、タンジェントなどを計算します。引数はすべてラジアン単位です。度に変換するには「度 * π / 180」の計算が必要です。

sin, sinf, sinl

概要: 引数のサイン（正弦）を計算します。

書式:
c
double sin(double x);
float sinf(float x); // C99以降
long double sinl(long double x); // C99以降

引数:
* x: ラジアン単位の角度。任意の実数値。

戻り値:
* x のサイン値（-1.0 から 1.0 の範囲）。
* x が NaN の場合、NaN を返します。
* x が ±INF の場合、定義されないため、NaN を返します（または定義域エラーを示す方法）。

詳細:
数学的なサイン関数 sin(x) を計算します。周期関数であり、周期は 2π です。

使用例:
“`c

include

include

int main() {
double angle_deg = 90.0;
double angle_rad = angle_deg * M_PI / 180.0; // 度をラジアンに変換 (M_PIは後述)
double result = sin(angle_rad);
printf(“sin(%f degrees) = %f\n”, angle_deg, result); // sin(90度) ≈ 1.0

float angle_rad_f = 0.785398f; // 約 45度
float result_f = sinf(angle_rad_f);
printf("sinf(%f radians) = %f\n", angle_rad_f, result_f); // sin(45度) ≈ 0.707

return 0;

}
``
**注意点:**M_PIは多くのシステムでやで定義されていますが、標準Cでは定義が保証されていません。移植性を重視する場合は、自分でM_PIを定義するか、atan(1.0) * 4.0` で計算するなどします。

cos, cosf, cosl

概要: 引数のコサイン（余弦）を計算します。

書式:
c
double cos(double x);
float cosf(float x); // C99以降
long double cosl(long double x); // C99以降

引数:
* x: ラジアン単位の角度。任意の実数値。

戻り値:
* x のコサイン値（-1.0 から 1.0 の範囲）。
* x が NaN の場合、NaN を返します。
* x が ±INF の場合、定義されないため、NaN を返します。

詳細:
数学的なコサイン関数 cos(x) を計算します。周期関数であり、周期は 2π です。

使用例:
“`c

include

include

int main() {
double angle_deg = 180.0;
double angle_rad = angle_deg * M_PI / 180.0;
double result = cos(angle_rad);
printf(“cos(%f degrees) = %f\n”, angle_deg, result); // cos(180度) ≈ -1.0

return 0;

}
“`

tan, tanf, tanl

概要: 引数のタンジェント（正接）を計算します。

書式:
c
double tan(double x);
float tanf(float x); // C99以降
long double tanl(long double x); // C99以降

引数:
* x: ラジアン単位の角度。π/2 + nπ (nは整数) 以外の任意の実数値。

戻り値:
* x のタンジェント値。
* x が π/2 + nπ に近い場合、戻り値は非常に大きくなる可能性があります。
* x が π/2 + nπ の場合、定義されないため、実装定義の値を返します（多くの場合 ±INF または NaN）。定義域エラー EDOM が発生する可能性もあります。
* x が NaN の場合、NaN を返します。
* x が ±INF の場合、定義されないため、NaN を返します。

詳細:
数学的なタンジェント関数 tan(x) = sin(x) / cos(x) を計算します。cos(x) がゼロになる点で不連続になります。

使用例:
“`c

include

include

int main() {
double angle_deg = 45.0;
double angle_rad = angle_deg * M_PI / 180.0;
double result = tan(angle_rad);
printf(“tan(%f degrees) = %f\n”, angle_deg, result); // tan(45度) ≈ 1.0

return 0;

}
“`

asin, asinf, asinl

概要: 引数のアークサイン（逆正弦）を計算します。

書式:
c
double asin(double x);
float asinf(float x); // C99以降
long double asinl(long double x); // C99以降

引数:
* x: -1.0 から 1.0 の範囲内の値。

戻り値:
* x のアークサイン値（ラジアン単位、-π/2 から π/2 の範囲）。
* x が定義域外の場合、定義域エラー EDOM が発生し、実装定義の値を返します（多くの場合 NaN）。
* x が NaN の場合、NaN を返します。
* x が ±0 の場合、±0 を返します。

詳細:
sin(y) = x となる y を計算します。sin 関数の逆関数ですが、一意に定まるように値域は -π/2 から π/2 に制限されています（主値）。

使用例:
“`c

include

include

int main() {
double value = 1.0;
double result_rad = asin(value);
double result_deg = result_rad * 180.0 / M_PI;
printf(“asin(%f) = %f radians (%f degrees)\n”, value, result_rad, result_deg); // asin(1) ≈ 90度

double invalid_value = 2.0;
double invalid_result = asin(invalid_value); // 定義域エラー
if (isnan(invalid_result)) { // C99以降のisnan関数
    printf("asin(%f) resulted in NaN (Domain Error)\n", invalid_value);
}

return 0;

}
“`
注意点: 定義域 [-1, 1] を超える値を渡すとエラーになります。エラー処理の方法は後述します。

acos, acosf, acosl

概要: 引数のアークコサイン（逆余弦）を計算します。

書式:
c
double acos(double x);
float acosf(float x); // C99以降
long double acosl(long double x); // C99以降

引数:
* x: -1.0 から 1.0 の範囲内の値。

戻り値:
* x のアークコサイン値（ラジアン単位、0 から π の範囲）。
* x が定義域外の場合、定義域エラー EDOM が発生し、実装定義の値を返します（多くの場合 NaN）。
* x が NaN の場合、NaN を返します。
* x が 1.0 の場合、+0 を返します。

詳細:
cos(y) = x となる y を計算します。cos 関数の逆関数で、値域は 0 から π に制限されています（主値）。

使用例:
“`c

include

include

int main() {
double value = 0.0;
double result_rad = acos(value);
double result_deg = result_rad * 180.0 / M_PI;
printf(“acos(%f) = %f radians (%f degrees)\n”, value, result_rad, result_deg); // acos(0) ≈ 90度

return 0;

}
“`

atan, atanf, atanl

概要: 引数のアークタンジェント（逆正接）を計算します。

書式:
c
double atan(double x);
float atanf(float x); // C99以降
long double atanl(long double x); // C99以降

引数:
* x: 任意の実数値。

戻り値:
* x のアークタンジェント値（ラジアン単位、-π/2 から π/2 の範囲）。
* x が NaN の場合、NaN を返します。
* x が ±0 の場合、±0 を返します。
* x が ±INF の場合、±π/2 を返します。

詳細:
tan(y) = x となる y を計算します。tan 関数の逆関数で、値域は -π/2 から π/2 に制限されています（主値）。atan(y/x) の計算が必要な場合は、後述の atan2 関数を使う方が、xが0に近い場合や象限を正しく判断できるため、より安定で正確です。

使用例:
“`c

include

include

int main() {
double value = 1.0;
double result_rad = atan(value);
double result_deg = result_rad * 180.0 / M_PI;
printf(“atan(%f) = %f radians (%f degrees)\n”, value, result_rad, result_deg); // atan(1) ≈ 45度

return 0;

}
“`

atan2, atan2f, atan2l

概要: 引数 y と x から、原点から点 (x, y) への角度を計算します。

書式:
c
double atan2(double y, double x);
float atan2f(float y, float x); // C99以降
long double atan2l(long double y, long double x); // C99以降

引数:
* y: 点のy座標。
* x: 点のx座標。

戻り値:
* 原点から点 (x, y) への角度（ラジアン単位、-π から π の範囲）。
* x が 0 かつ y が 0 の場合、結果は実装定義（多くの場合 0 または定義域エラー EDOM）。
* x が正の場合、atan(y/x) と同じです。
* x が負の場合、atan(y/x) ± π です（yの符号による）。
* y が 0 かつ x が負の場合、π を返します。
* y が 0 かつ x が正の場合、±0 を返します（yの符号による）。
* y が正、x が 0 の場合、+π/2 を返します。
* y が負、x が 0 の場合、-π/2 を返します。
* どちらかの引数が NaN の場合、NaN を返します。

詳細:
この関数は atan(y/x) よりも強力です。引数 x と y の両方の符号を考慮して、点 (x, y) がどの象限にあるかを正確に判断し、適切な角度を -π から π の範囲で返します。これにより、 y/x の計算で x がゼロになる場合の除算エラーや、 y/x の値だけでは象限を区別できない問題を回避できます。平面上の点やベクトルの角度を計算する際に非常に役立ちます。

使用例:
“`c

include

include

int main() {
// 点 (1, 1) -> 第1象限
printf(“atan2(1.0, 1.0) = %f radians\n”, atan2(1.0, 1.0)); // 約 0.785 (π/4)

// 点 (-1, 1) -> 第2象限
printf("atan2(1.0, -1.0) = %f radians\n", atan2(1.0, -1.0)); // 約 2.356 (3π/4)

// 点 (-1, -1) -> 第3象限
printf("atan2(-1.0, -1.0) = %f radians\n", atan2(-1.0, -1.0)); // 約 -2.356 (-3π/4)

// 点 (1, -1) -> 第4象限
printf("atan2(-1.0, 1.0) = %f radians\n", atan2(-1.0, 1.0)); // 約 -0.785 (-π/4)

// 点 (0, 1) -> y軸正方向
printf("atan2(1.0, 0.0) = %f radians\n", atan2(1.0, 0.0)); // 約 1.570 (π/2)

return 0;

}
“`

3.2. 双曲線関数 (Hyperbolic Functions)

自然対数の底 e を用いて定義される双曲線サイン、コサイン、タンジェントなどを計算します。これらの関数は、懸垂線や特殊相対性理論など、様々な分野で現れます。

sinh, sinhf, sinhl

概要: 引数の双曲線サインを計算します。

書式:
c
double sinh(double x);
float sinhf(float x); // C99以降
long double sinhl(long double x); // C99以降

引数:
* x: 任意の実数値。

戻り値:
* x の双曲線サイン値 ((e^x – e^-x) / 2)。
* 結果が表現可能な範囲を超えた場合、極エラー ERANGE が発生し、±INF を返します。
* x が NaN の場合、NaN を返します。
* x が ±0 の場合、±0 を返します。
* x が ±INF の場合、±INF を返します。

詳細:
sinh(x) = (e^x – e^-x) / 2 です。

使用例:
“`c

include

include

int main() {
double value = 1.0;
printf(“sinh(%f) = %f\n”, value, sinh(value)); // sinh(1) ≈ 1.175

return 0;

}
“`

cosh, coshf, coshl

概要: 引数の双曲線コサインを計算します。

書式:
c
double cosh(double x);
float coshf(float x); // C99以降
long double coshl(long double x); // C99以降

引数:
* x: 任意の実数値。

戻り値:
* x の双曲線コサイン値 ((e^x + e^-x) / 2)。
* 結果が表現可能な範囲を超えた場合、極エラー ERANGE が発生し、+INF を返します。
* x が NaN の場合、NaN を返します。
* x が ±0 の場合、1.0 を返します。
* x が ±INF の場合、+INF を返します。

詳細:
cosh(x) = (e^x + e^-x) / 2 です。常に 1 以上の値を返します。

使用例:
“`c

include

include

int main() {
double value = 0.0;
printf(“cosh(%f) = %f\n”, value, cosh(value)); // cosh(0) = 1.0

return 0;

}
“`

tanh, tanhf, tanhl

概要: 引数の双曲線タンジェントを計算します。

書式:
c
double tanh(double x);
float tanhf(float x); // C99以降
long double tanhl(long double x); // C99以降

引数:
* x: 任意の実数値。

戻り値:
* x の双曲線タンジェント値 ((e^x – e^-x) / (e^x + e^-x))。
* x が NaN の場合、NaN を返します。
* x が ±0 の場合、±0 を返します。
* x が +INF の場合、+1.0 を返します。
* x が -INF の場合、-1.0 を返します。

詳細:
tanh(x) = sinh(x) / cosh(x) です。戻り値は常に -1.0 から 1.0 の範囲です。

使用例:
“`c

include

include

int main() {
double value = 1.0;
printf(“tanh(%f) = %f\n”, value, tanh(value)); // tanh(1) ≈ 0.761

return 0;

}
“`

（逆双曲線関数 asinh, acosh, atanh もありますが、使用頻度はやや低いため、ここでは詳細な説明は省略します。これらもC99以降で float 版と long double 版が追加されています。）

3.3. 指数関数と対数関数 (Exponential and Logarithmic Functions)

自然対数の底 e や 10、2 を底とする指数関数および対数関数を計算します。

exp, expf, expl

概要: 自然対数の底 e を引数で指定したべき乗に計算します (e^x)。

書式:
c
double exp(double x);
float expf(float x); // C99以降
long double expl(long double x); // C99以降

引数:
* x: 任意の実数値。

戻り値:
* e の x 乗の値。
* 結果が表現可能な範囲を超えた場合、極エラー ERANGE が発生し、±INF を返します（符号は x による）。
* 結果が小さすぎて非正規化数としても表現できない場合、極エラー ERANGE が発生し、±0 を返します。
* x が NaN の場合、NaN を返します。
* x が +INF の場合、+INF を返します。
* x が -INF の場合、+0 を返します。

詳細:
数学的な指数関数 e^x を計算します。eは約 2.71828 です。

使用例:
“`c

include

include

int main() {
double value = 1.0;
printf(“exp(%f) = %f\n”, value, exp(value)); // exp(1) = e ≈ 2.718

return 0;

}
“`

log, logf, logl

概要: 引数の自然対数（底が e）を計算します。

書式:
c
double log(double x);
float logf(float x); // C99以降
long double logl(long double x); // C99以降

引数:
* x: 0より大きい実数値。

戻り値:
* x の自然対数の値。
* x が 0 の場合、極エラー ERANGE が発生し、-INF を返します。
* x が負の場合、定義域エラー EDOM が発生し、実装定義の値を返します（多くの場合 NaN）。
* x が NaN の場合、NaN を返します。
* x が +1 の場合、+0 を返します。
* x が +INF の場合、+INF を返します。

詳細:
数学的な自然対数 ln(x) または log_e(x) を計算します。

使用例:
“`c

include

include

int main() {
double value = exp(1.0); // 約 2.718
printf(“log(%f) = %f\n”, value, log(value)); // log(e) ≈ 1.0

double invalid_value = -1.0;
double invalid_result = log(invalid_value); // 定義域エラー
if (isnan(invalid_result)) {
    printf("log(%f) resulted in NaN (Domain Error)\n", invalid_value);
}

return 0;

}
“`
注意点: 引数が非正の場合はエラーになります。

log10, log10f, log10l

概要: 引数の常用対数（底が 10）を計算します。

書式:
c
double log10(double x);
float log10f(float x); // C99以降
long double log10l(long double x); // C99以降

引数:
* x: 0より大きい実数値。

戻り値:
* x の常用対数の値。
* x が 0 の場合、極エラー ERANGE が発生し、-INF を返します。
* x が負の場合、定義域エラー EDOM が発生し、実装定義の値を返します（多くの場合 NaN）。
* x が NaN の場合、NaN を返します。
* x が +1 の場合、+0 を返します。
* x が +INF の場合、+INF を返します。

詳細:
数学的な常用対数 log_10(x) を計算します。音の大きさ（デシベル）や地震の規模（マグニチュード）など、対数目盛が使われる分野でよく利用されます。

使用例:
“`c

include

include

int main() {
double value = 100.0;
printf(“log10(%f) = %f\n”, value, log10(value)); // log10(100) = 2.0

return 0;

}
“`
注意点: 引数が非正の場合はエラーになります。

exp2, exp2f, exp2l (C99以降)

概要: 2を引数で指定したべき乗に計算します (2^x)。

書式:
c
double exp2(double x);
float exp2f(float x);
long double exp2l(long double x);

詳細: 2のべき乗を計算します。特に浮動小数点数の内部表現（仮数と指数部）に関連する計算で便利です。

log1p, log1pf, log1pl (C99以降)

概要: log(1 + x) を計算します。

書式:
c
double log1p(double x);
float log1pf(float x);
long double log1pl(long double x);

詳細: x がゼロに近い場合、 log(1.0 + x) を直接計算するよりも精度良く計算できます。これは、1+x が非常に1に近いと、加算によって有効数字が失われる（桁落ち）可能性があるためです。

log2, log2f, log2l (C99以降)

概要: 引数の2を底とする対数を計算します。

書式:
c
double log2(double x);
float log2f(float x);
long double log2l(long double x);

詳細: 情報理論（エントロピーなど）やコンピュータサイエンスでよく使われます。

3.4. べき乗と平方根 (Power and Root Functions)

任意の底に対するべき乗や、平方根、立方根などを計算します。

pow, powf, powl

概要: x の y 乗 (x^y) を計算します。

書式:
c
double pow(double x, double y);
float powf(float x, float y); // C99以降
long double powl(long double x, long double y); // C99以降

引数:
* x: 底となる実数値。
* y: 指数となる実数値。

戻り値:
* x の y 乗の値。
* x が負かつ y が非整数の場合、定義域エラー EDOM が発生し、実装定義の値を返します（多くの場合 NaN）。
* x が 0 かつ y が負または 0 の場合、定義域エラー EDOM が発生する可能性があります。y が負の場合は極エラー ERANGE も発生し、+INF を返します。
* 結果が表現可能な範囲を超えた場合、極エラー ERANGE が発生し、±INF を返します（符号は結果による）。
* 結果が小さすぎて非正規化数としても表現できない場合、極エラー ERANGE が発生し、±0 を返します。
* その他の特殊なケース（例: pow(1, y) は 1、pow(x, 0) は 1、pow(INF, y) など）は標準で定められています。

詳細:
数学的な x^y を計算します。浮動小数点数のべき乗であるため、指数が整数の場合でも浮動小数点演算として処理されることに注意が必要です。負の数の非整数べき乗は実数の範囲では定義されないため、エラーとなります。

使用例:
“`c

include

include

int main() {
printf(“pow(2.0, 3.0) = %f\n”, pow(2.0, 3.0)); // 2^3 = 8.0
printf(“pow(9.0, 0.5) = %f\n”, pow(9.0, 0.5)); // 9^0.5 = sqrt(9) = 3.0
printf(“pow(-2.0, 3.0) = %f\n”, pow(-2.0, 3.0)); // (-2)^3 = -8.0

// 負の数の非整数べき乗 (定義域エラー)
double invalid_result = pow(-2.0, 0.5);
if (isnan(invalid_result)) {
     printf("pow(-2.0, 0.5) resulted in NaN (Domain Error)\n");
}

return 0;

}
“`

sqrt, sqrtf, sqrtl

概要: 引数の平方根を計算します。

書式:
c
double sqrt(double x);
float sqrtf(float x); // C99以降
long double sqrtl(long double x); // C99以降

引数:
* x: 0以上の実数値。

戻り値:
* x の非負の平方根。
* x が負の場合、定義域エラー EDOM が発生し、実装定義の値を返します（多くの場合 NaN）。
* x が NaN の場合、NaN を返します。
* x が ±0 の場合、±0 を返します。
* x が +INF の場合、+INF を返します。

詳細:
数学的な平方根 √x を計算します。実数の範囲では非負の数に対してのみ定義されます。

使用例:
“`c

include

include

int main() {
double value = 16.0;
printf(“sqrt(%f) = %f\n”, value, sqrt(value)); // sqrt(16) = 4.0

double invalid_value = -4.0;
double invalid_result = sqrt(invalid_value); // 定義域エラー
if (isnan(invalid_result)) {
    printf("sqrt(%f) resulted in NaN (Domain Error)\n", invalid_value);
}

return 0;

}
“`
注意点: 引数が負の場合はエラーになります。

cbrt, cbrtf, cbrtl (C99以降)

概要: 引数の立方根を計算します。

書式:
c
double cbrt(double x);
float cbrtf(float x);
long double cbrtl(long double x);

詳細: 平方根と異なり、立方根は負の数に対しても実数の値が定義されます。例えば、cbrt(-8.0) は -2.0 です。

使用例:
“`c

include

include

int main() {
printf(“cbrt(8.0) = %f\n”, cbrt(8.0)); // 2.0
printf(“cbrt(-8.0) = %f\n”, cbrt(-8.0)); // -2.0

return 0;

}
“`

hypot, hypotf, hypotl (C99以降)

概要: 二つの引数 x と y から、直角三角形の斜辺の長さを計算します (sqrt(x^2 + y^2))。

書式:
c
double hypot(double x, double y);
float hypotf(float x, float y);
long double hypotl(long double x, long double y);

詳細: sqrt(x*x + y*y) と同じ計算ですが、hypot 関数は x や y が非常に大きい場合に発生しやすい中間計算でのオーバーフローを避けるように設計されています。点 (x, y) と原点の間のユークリッド距離を計算するのに便利です。

使用例:
“`c

include

include

int main() {
double x = 3.0, y = 4.0;
printf(“hypot(%f, %f) = %f\n”, x, y, hypot(x, y)); // sqrt(33 + 44) = sqrt(9+16) = sqrt(25) = 5.0

return 0;

}
“`

3.5. 絶対値 (Absolute Value Functions)

浮動小数点数の絶対値を計算します。整数型の絶対値関数 (abs, labs, llabs) は <stdlib.h> に定義されていることに注意してください。

fabs, fabsf, fabsl

概要: 引数の浮動小数点絶対値を計算します。

書式:
c
double fabs(double x);
float fabsf(float x); // C99以降
long double fabsl(long double x); // C99以降

引数:
* x: 任意の実数値。

戻り値:
* x の絶対値 (|x|)。
* x が NaN の場合、NaN を返します。
* x が ±0 の場合、+0 を返します。
* x が ±INF の場合、+INF を返します。

詳細:
引数の符号を取り除き、常に非負の値を返します。

使用例:
“`c

include

include

int main() {
double value = -123.45;
printf(“fabs(%f) = %f\n”, value, fabs(value)); // 123.45

return 0;

}
“`

3.6. 丸め関数 (Rounding Functions)

浮動小数点数を特定の規則に従って整数値に丸めます。戻り値は浮動小数点数型ですが、値は整数になります。

ceil, ceilf, ceill

概要: 引数以上の最小の整数値を計算します（切り上げ）。

書式:
c
double ceil(double x);
float ceilf(float x); // C99以降
long double ceill(long double x); // C99以降

引数:
* x: 任意の実数値。

戻り値:
* x 以上の最小の整数値を浮動小数点数型で返します。
* x が NaN の場合、NaN を返します。
* x が ±INF または ±0 の場合、x そのものを返します。

詳細:
数直線上において、x または x より大きく、最も x に近い整数を返します。

使用例:
“`c

include

include

int main() {
printf(“ceil(4.2) = %f\n”, ceil(4.2)); // 5.0
printf(“ceil(4.0) = %f\n”, ceil(4.0)); // 4.0
printf(“ceil(-4.2) = %f\n”, ceil(-4.2)); // -4.0
printf(“ceil(-4.0) = %f\n”, ceil(-4.0)); // -4.0

return 0;

}
“`

floor, floorf, floorl

概要: 引数以下の最大の整数値を計算します（切り捨て）。

書式:
c
double floor(double x);
float floorf(float x); // C99以降
long double floorl(long double x); // C99以降

引数:
* x: 任意の実数値。

戻り値:
* x 以下の最大の整数値を浮動小数点数型で返します。
* x が NaN の場合、NaN を返します。
* x が ±INF または ±0 の場合、x そのものを返します。

詳細:
数直線上において、x または x より小さく、最も x に近い整数を返します。正の数に対しては小数点以下を切り捨てるように見えますが、負の数に対してはゼロから遠ざかる方向に切り捨てられます。

使用例:
“`c

include

include

int main() {
printf(“floor(4.2) = %f\n”, floor(4.2)); // 4.0
printf(“floor(4.0) = %f\n”, floor(4.0)); // 4.0
printf(“floor(-4.2) = %f\n”, floor(-4.2)); // -5.0
printf(“floor(-4.0) = %f\n”, floor(-4.0)); // -4.0

return 0;

}
“`

round, roundf, roundl (C99以降)

概要: 引数を最も近い整数に丸めます（四捨五入）。中間値 (x.5) はゼロから遠い方の整数に丸められます。

書式:
c
double round(double x);
float roundf(float x);
long double roundl(long double x);

詳細:
例えば、round(4.5) は 5.0 に、round(-4.5) は -5.0 に丸められます。これは学校で習う一般的な四捨五入（小数点以下が0.5以上の場合は切り上げ）とは負の数に対する挙動が異なります。正負どちらも絶対値が大きな方に丸められます。

使用例:
“`c

include

include

int main() {
printf(“round(4.2) = %f\n”, round(4.2)); // 4.0
printf(“round(4.5) = %f\n”, round(4.5)); // 5.0
printf(“round(4.8) = %f\n”, round(4.8)); // 5.0
printf(“round(-4.2) = %f\n”, round(-4.2)); // -4.0
printf(“round(-4.5) = %f\n”, round(-4.5)); // -5.0
printf(“round(-4.8) = %f\n”, round(-4.8)); // -5.0

return 0;

}
“`

trunc, truncf, truncl (C99以降)

概要: 引数の小数部を切り捨て、ゼロ方向への整数値を計算します。

書式:
c
double trunc(double x);
float truncf(float x);
long double truncl(long double x);

詳細:
正の数に対しては floor と同じ、負の数に対しては ceil と同じ挙動をします。つまり、常にゼロに近づく方向に丸められます。

使用例:
“`c

include

include

int main() {
printf(“trunc(4.2) = %f\n”, trunc(4.2)); // 4.0
printf(“trunc(-4.2) = %f\n”, trunc(-4.2)); // -4.0

return 0;

}
“`

nearbyint, rint (C99以降)

これらの関数は、現在の浮動小数点丸めモードに従って引数を丸めます。round は常に「ゼロから遠い方へ」丸めるのに対し、これらの関数は FE_UPWARD, FE_DOWNWARD, FE_TOWARDZERO, FE_TONEAREST といった丸めモード設定（<fenv.h> で制御）に従います。rint は inexact 例外を発生させる可能性がある点で nearbyint と異なります。

3.7. 剰余関数 (Remainder Functions)

浮動小数点数に対する剰余を計算します。

fmod, fmodf, fmodl

概要: x を y で割った浮動小数点数の剰余を計算します。結果の符号は x と同じになります。

書式:
c
double fmod(double x, double y);
float fmodf(float x, float y); // C99以降
long double fmodl(long double x, long double y); // C99以降

詳細:
x - n * y を計算します。ここで n は x / y をゼロ方向に切り捨てた整数です。結果の符号は x の符号と同じになり、絶対値は y の絶対値より小さくなります。

使用例:
“`c

include

include

int main() {
printf(“fmod(5.3, 2.0) = %f\n”, fmod(5.3, 2.0)); // 5.3 = 2 * 2.0 + 1.3 -> 1.3
printf(“fmod(-5.3, 2.0) = %f\n”, fmod(-5.3, 2.0)); // -5.3 = -2 * 2.0 + -1.3 -> -1.3
printf(“fmod(5.3, -2.0) = %f\n”, fmod(5.3, -2.0)); // 5.3 = -2 * -2.0 + 1.3 -> 1.3

return 0;

}
“`

remainder, remainderf, remainderl (C99以降)

概要: x を y で割った浮動小数点数の剰余を計算します。結果の符号は x と同じになり、絶対値は y / 2 以下になります。

書式:
c
double remainder(double x, double y);
float remainderf(float x, float y);
long double remainderl(long double x, long double y);

詳細:
x - n * y を計算します。ここで n は x / y を最も近い整数に丸めた値です。丸める際にちょうど中間 (x.5) の場合は偶数方向へ丸められます。fmod よりも対称性の良い剰余を返します。周期関数などで便利です。

使用例:
“`c

include

include

int main() {
printf(“remainder(5.3, 2.0) = %f\n”, remainder(5.3, 2.0)); // 5.3/2.0=2.65, 最も近い整数は3 -> 5.3 – 32.0 = -0.7
printf(“remainder(5.0, 2.0) = %f\n”, remainder(5.0, 2.0)); // 5.0/2.0=2.5, 最も近い偶数整数は2 -> 5.0 – 22.0 = 1.0
printf(“fmod(5.0, 2.0) = %f\n”, fmod(5.0, 2.0)); // 5.0/2.0=2.5, ゼロ方向へ切り捨て整数は2 -> 5.0 – 2*2.0 = 1.0 (この場合はfmodと同じ)

return 0;

}
“`

3.8. 浮動小数点操作関数 (Floating-Point Manipulation Functions)

浮動小数点数をその構成要素（仮数、指数）に分解したり、それらを操作したりする関数です。

frexp, frexpf, frexpl

概要: 浮動小数点数 value を、正規化された仮数 (mantissa) と 2 を底とする指数 (exp) に分解します。value = mantissa * 2^exp となるようにします。

書式:
c
double frexp(double value, int *exp);
float frexpf(float value, int *exp); // C99以降
long double frexpl(long double value, int *exp); // C99以降

詳細:
mantissa は [0.5, 1.0) または [-1.0, -0.5) の範囲（値が0の場合は0）となり、exp は対応する整数指数です。exp はポインタで渡された整数変数に格納されます。浮動小数点数の内部表現を調べたり、正規化された値を扱ったりするのに便利です。

使用例:
“`c

include

include

int main() {
double value = 12.5; // バイナリで 1100.1 = 1.1001 * 2^3
int exponent;
double mantissa = frexp(value, &exponent);
printf(“frexp(%f) = mantissa %f, exponent %d (mantissa * 2^exponent = %f)\n”,
value, mantissa, exponent, mantissa * pow(2.0, exponent)); // 12.5 = 0.78125 * 2^4

// IEEE 754では仮数は [1.0, 2.0) または [0.0, 1.0)（非正規化数の場合）
// frexpの仕様は実装によって[0.5, 1.0)または[0, 1)など異なる場合があります。
// 上の例では、12.5(1100.1_2) -> 1.1001_2 * 2^3。 frexpは0.11001 * 2^4と分解することが多いです。
// 0.11001_2 = 1/2 + 1/4 + 0/8 + 0/16 + 1/32 = 0.5 + 0.25 + 0.03125 = 0.78125

return 0;

}
“`

ldexp, ldexpf, ldexpl

概要: x * 2^exp を計算します。

書式:
c
double ldexp(double x, int exp);
float ldexpf(float x, int exp); // C99以降
long double ldexpl(long double x, int exp); // C99以降

詳細:
frexp の逆の操作を行います。仮数部と指数部から浮動小数点数を生成するのに使えます。x * 2^exp はビット演算での左シフトに似ていますが、浮動小数点数に対して機能し、指数部のみを変更することで効率的に計算が行われます。

使用例:
“`c

include

include

int main() {
double mantissa = 0.78125;
int exponent = 4;
double value = ldexp(mantissa, exponent);
printf(“ldexp(%f, %d) = %f\n”, mantissa, exponent, value); // 0.78125 * 2^4 = 12.5

return 0;

}
“`

modf, modff, modfl

概要: 引数 value を整数部と小数部に分解します。結果の符号は value と同じになります。

書式:
c
double modf(double value, double *iptr);
float modff(float value, float *iptr); // C99以降
long double modfl(long double value, long double *iptr); // C99以降

詳細:
value の整数部（小数点以下切り捨て、符号は元の値と同じ）が iptr が指す変数に格納され、小数部が関数の戻り値として返されます。戻り値は value から整数部を引いた値になります。

使用例:
“`c

include

include

int main() {
double value = -123.456;
double integer_part;
double fractional_part = modf(value, &integer_part);
printf(“modf(%f) = integer part %f, fractional part %f\n”,
value, integer_part, fractional_part); // integer part -123.000, fractional part -0.456000

return 0;

}
“`

（その他、C99以降で scalbn, scalbln, ilogb, logb, nextafter, nexttoward, copysign などがあります。これらは浮動小数点数の指数部の操作や、隣接する表現可能な数の取得、符号のコピーなど、より低レベルまたは詳細な操作を提供します。）

3.9. 分類・比較関数 (Classification and Comparison Functions) (C99以降)

浮動小数点数が正常な数、無限大、非数（NaN）のいずれであるかなどを判定したり、NaN が含まれる可能性のある浮動小数点数を安全に比較したりするための関数です。これらの関数は通常 <math.h> と共に <fenv.h> で定義される浮動小数点環境フラグと組み合わせて使用されます。

主要な分類マクロ/関数:
* fpclassify(x): x の種類を分類します（FP_NAN, FP_INFINITE, FP_NORMAL, FP_SUBNORMAL, FP_ZERO のいずれか）。
* isfinite(x): x が有限値であれば真を返します（NaN や INF でなければ真）。
* isinf(x): x が無限大 (±INF) であれば真を返します。
* isnan(x): x が非数 (NaN) であれば真を返します。
* isnormal(x): x が正規化数であれば真を返します。
* signbit(x): x の符号ビットが立っていれば（負であれば）真を返します。

安全な比較マクロ (NaNを考慮):
通常の比較演算子 (<, >, <=, >=, ==, !=) は、オペランドに NaN が含まれる場合、常に偽を返します（!= の場合は真）。これにより予期しない結果になることがあります。これらのマクロは、NaN が含まれる場合は例外を発生させることなく、指定された比較の結果を返します。
* isgreater(x, y): x > y なら真
* isgreaterequal(x, y): x >= y なら真
* isless(x, y): x < y なら真
* islessequal(x, y): x <= y なら真
* islessgreater(x, y): x < y または x > y なら真 (x != y と同じだが、NaNの場合は偽)
* isunordered(x, y): x または y が NaN の場合、真を返します (比較が順序付けられない場合)。

使用例:
“`c

include

include

int main() {
double finite_val = 10.0;
double inf_val = 1.0 / 0.0; // 無限大
double nan_val = 0.0 / 0.0; // 非数

printf("%f is finite: %d\n", finite_val, isfinite(finite_val)); // 真 (1)
printf("%f is infinite: %d\n", finite_val, isinf(finite_val));   // 偽 (0)
printf("%f is NaN: %d\n", finite_val, isnan(finite_val));     // 偽 (0)

printf("%f is finite: %d\n", inf_val, isfinite(inf_val));     // 偽 (0)
printf("%f is infinite: %d\n", inf_val, isinf(inf_val));     // 真 (1)
printf("%f is NaN: %d\n", inf_val, isnan(inf_val));         // 偽 (0)

printf("%f is finite: %d\n", nan_val, isfinite(nan_val));     // 偽 (0)
printf("%f is infinite: %d\n", nan_val, isinf(nan_val));     // 偽 (0)
printf("%f is NaN: %d\n", nan_val, isnan(nan_val));         // 真 (1)

// 安全な比較例
printf("1.0 > 0.0 : %d\n", isgreater(1.0, 0.0)); // 真 (1)
printf("1.0 < 0.0 : %d\n", isless(1.0, 0.0));    // 偽 (0)
printf("NaN > 0.0 : %d\n", isgreater(nan_val, 0.0)); // 偽 (0)
printf("NaN < 0.0 : %d\n", isless(nan_val, 0.0));    // 偽 (0)
printf("NaN == NaN : %d\n", (nan_val == nan_val)); // 偽 (0) - 通常の比較演算子
printf("NaN is unordered with NaN : %d\n", isunordered(nan_val, nan_val)); // 真 (1)

return 0;

}
“`

3.10. 特殊関数 (Special Functions) (C99以降)

より専門的な数学分野で用いられる関数の一部が <math.h> に含まれています。

• erf, erfc: 誤差関数とその補関数。確率・統計学で正規分布に関連して現れます。
• tgamma: ガンマ関数 (Γ(x))。階乗の一般化など。
• lgamma: ガンマ関数の絶対値の自然対数 (log(|Γ(x)|))。Γ(x) が非常に大きな値になる場合にオーバーフローを避けるために使われます。

これらの関数は特定の分野で重要ですが、一般的なプログラミングで日常的に使用されることは少ないかもしれません。

4. 数学定数 (Mathematical Constants)

<math.h> は、円周率 π や自然対数の底 e といったよく使われる数学定数を定義していることがあります。ただし、標準Cではこれらの定数の定義は保証されていません。多くの処理系では慣習として M_PI や M_E といったマクロ名で定義されていますが、これらはPOSIX標準やGNU拡張などに由来するものであり、厳密な標準C準拠の環境では利用できない場合があります。

移植性を最優先する場合は、これらの定数を自分で定義するか、計算によって得る必要があります。例えば、円周率 π は atan(1.0) * 4.0 で求めることができます。

使用例 (多くの環境で動作しますが非標準):
“`c

include // M_PI が定義されている場合

include

int main() {
// M_PI が定義されているか確認 (多くのシステムでは STDC_WANT_LIB_EXT1 または _USE_MATH_DEFINES が必要かもしれない)

ifdef M_PI

printf("円周率 M_PI = %f\n", M_PI);

else

// M_PI が定義されていない場合の代替手段
printf("M_PI is not defined. Calculating pi: %f\n", atan(1.0) * 4.0);

endif

return 0;

}
“`

5. 浮動小数点環境とエラー処理

数学関数を使用する上で、エラー処理は重要な側面です。<math.h> の関数は、引数が定義域外であったり（例: sqrt(-1.0)）、結果が表現可能な範囲を超えたり（例: exp(1000.0)）した場合にエラーを報告します。

C言語におけるエラー報告メカニズムには、主に二つの方法があります。

1. errno マクロ (<errno.h>): 関数がエラーを検出した場合、グローバル変数 errno に特定の値（<errno.h> で定義される EDOM や ERANGE など）を設定します。プログラマは関数の呼び出し後に errno の値をチェックしてエラーを検出します。ただし、errno は他の標準ライブラリ関数でも使用されるため、数学関数呼び出しの前に errno をゼロにクリアしておくのが一般的です。

   • EDOM (Domain Error): 引数が関数の定義域外である場合に発生します。
   • ERANGE (Range Error): 計算結果が表現可能な範囲外（オーバーフローやアンダーフロー）である場合に発生します。

   使用例:
   “`c

   include

   include

   include // errno を使用するために必要

   int main() {
   errno = 0; // errno をクリア
   double result = sqrt(-1.0);

   if (errno == EDOM) {
       perror("sqrt(-1.0)"); // エラーメッセージを出力
   } else {
       printf("sqrt(-1.0) = %f\n", result); // NaN が出力されることが多い
   }
   
   errno = 0; // errno をクリア
   double big_result = exp(1000.0);
   
   if (errno == ERANGE) {
       perror("exp(1000.0)"); // エラーメッセージを出力
   } else {
       printf("exp(1000.0) = %f\n", big_result); // INF が出力されることが多い
   }
   
   return 0;
   

   }
   “`

2. 浮動小数点例外 (<fenv.h> と C99以降の math_errhandling): C99以降、浮動小数点環境に関する標準ヘッダー <fenv.h> が導入されました。これにより、丸めモードの制御や、浮動小数点演算で発生する例外（不正確な演算 FE_INEXACT, ゼロ除算 FE_DIVBYZERO, オーバーフロー FE_OVERFLOW, アンダーフロー FE_UNDERFLOW, 無効な演算 FE_INVALID など）のフラグ操作が可能になりました。

   <math.h> の関数は、これらの例外フラグを設定することでエラーを報告することもあります。C99標準では、マクロ math_errhandling が定義されており、処理系が errno を使用するか (MATH_ERRNO)、浮動小数点例外を使用するか (MATH_ERREXCEPT)、またはその両方を使用するかを示します。

   使用例 (C99以降, <fenv.h> 使用):
   “`c

   include

   include

   include // 浮動小数点環境を操作するために必要

   pragma STDC FENV_ACCESS ON // 浮動小数点環境へのアクセスを許可

   int main() {
   // FE_INVALID 例外フラグをクリア
   feclearexcept(FE_INVALID);

   double result = sqrt(-1.0);
   
   // FE_INVALID フラグが設定されたかチェック
   if (fetestexcept(FE_INVALID)) {
       printf("sqrt(-1.0) resulted in FE_INVALID exception.\n");
   } else {
       printf("sqrt(-1.0) = %f\n", result);
   }
   
   // FE_OVERFLOW 例外フラグをクリア
   feclearexcept(FE_OVERFLOW);
   
   double big_result = exp(1000.0);
   
   // FE_OVERFLOW フラグが設定されたかチェック
   if (fetestexcept(FE_OVERFLOW)) {
        printf("exp(1000.0) resulted in FE_OVERFLOW exception.\n");
   } else {
        printf("exp(1000.0) = %f\n", big_result);
   }
   
   // math_errhandling の値を確認
   printf("math_errhandling: ");
   if (math_errhandling & MATH_ERRNO) {
       printf("MATH_ERRNO ");
   }
   if (math_errhandling & MATH_ERREXCEPT) {
       printf("MATH_ERREXCEPT");
   }
   printf("\n");
   
   return 0;
   

   }
   ``
**注意点:**の機能は処理系のサポートに依存します。また、#pragma STDC FENV_ACCESS ON` を使用して、コンパイラの最適化が浮動小数点環境の副作用（例外フラグの設定など）を無視しないように指示する必要があります。

どちらの方法が使用されるかは処理系に依存しますが、現代的な環境では MATH_ERREXCEPT をサポートしていることが多いです。信頼性の高い数値計算を行うには、これらのエラー報告メカニズムを理解し、適切にエラーをチェックまたは処理することが重要です。また、NaNやINFといった特殊な戻り値のチェックには、C99以降の分類関数 (isnan, isinf など) が役立ちます。

6. データ型と型総称マクロ (tgmath.h)

前述のように、C99標準では float, double, long double それぞれに対応する数学関数が追加されました (sinf, sin, sinl など)。これにより、プログラマは使用する浮動小数点数の型に合わせて関数を選ぶことができます。

しかし、同じ計算（例えばサイン）を異なる型の変数に対して行いたい場合に、変数の型に応じて sinf, sin, sinl を使い分けるのは煩雑です。この問題を解決するため、C99では <tgmath.h> ヘッダーファイルが導入されました。

<tgmath.h> をインクルードすると、多くの <math.h> および <complex.h>（複素数関数）の関数名と同じ名前の「型総称マクロ」が利用可能になります。これらのマクロは、引数の型を調べ、自動的に対応する float, double, long double または複素数型の関数を選択して呼び出します。

例えば、<tgmath.h> をインクルードした状態で sin(x) を呼び出すと、以下のようになります。

• x が float 型の場合、 sinf(x) が呼び出される
• x が double 型の場合、 sin(x) (元の double 版関数) が呼び出される
• x が long double 型の場合、 sinl(x) が呼び出される

これにより、プログラマは型の違いを意識することなく、同じマクロ名で数学関数を利用できます。

使用例 (C99以降, <tgmath.h>):
“`c

include // 型総称マクロを使用

include

int main() {
float f_val = 1.0f;
double d_val = 1.0;
long double ld_val = 1.0L;

// tgmath.h により、引数の型に応じて適切な sin 関数が選択される
printf("sin(float 1.0) = %f\n", sin(f_val));
printf("sin(double 1.0) = %lf\n", sin(d_val));
printf("sin(long double 1.0) = %Lf\n", sin(ld_val));

// sqrt も同様
printf("sqrt(float 4.0) = %f\n", sqrt(4.0f));
printf("sqrt(double 9.0) = %lf\n", sqrt(9.0));
printf("sqrt(long double 16.0) = %Lf\n", sqrt(16.0L));

return 0;

}
``` はコードの簡潔さと可読性を向上させ、型の変更に対するプログラムの堅牢性を高めます。ただし、これはC99以降の標準であるため、古い環境では利用できない場合があります。

7. コンパイルとリンク

<math.h> で宣言されている関数は、通常、C標準ライブラリの一部として提供されますが、多くの場合、独立した数学ライブラリ (libm) に実装されています。そのため、プログラムをコンパイルする際に、この数学ライブラリをリンクする必要があります。

GCCなどの多くのコンパイラでは、数学関数を使用するソースファイルをコンパイルする際に、-lm オプションをリンカーに渡す必要があります。

例:
bash
gcc your_program.c -o your_program -lm
-lm は “link math library” を意味します。これにより、コンパイラは数学関数に必要な定義を libm から探し出して実行ファイルに含める（または動的リンクする）ようになります。

ただし、環境によっては -lm オプションが不要な場合もあります（例: Windows上のMinGWなど、一部の環境では標準Cライブラリに数学関数が統合されている）。しかし、多くのUnix系システムでは -lm が必要ですので、習慣として付けておくのが安全です。

8. 応用分野

<math.h> の関数は、非常に多くの分野で利用されています。

• 物理学・工学: 力学計算、波動解析、信号処理、制御理論など、数値モデルやシミュレーションに不可欠です。
• 統計学: 確率分布関数の計算、データ解析など。
• コンピュータグラフィックス: 3D座標変換、照明計算、テクスチャマッピングなど、ベクトルの操作や角度計算に三角関数やユークリッド距離などが使われます。
• ゲーム開発: 物理エンジン、AIの経路探索、アニメーション計算など、様々な数値計算に利用されます。
• 金融工学: オプション価格計算など、数学的なモデルに基づいた計算。
• 最適化問題: 数学的関数の最小値・最大値を求めるアルゴリズムなどで使用されます。

<math.h> はこれらの応用分野で必要となる基本的な数学的演算を提供し、より複雑なアルゴリズムを構築するための土台となります。

9. 精度の問題と注意点

浮動小数点数の計算は、常に精度に関する注意が必要です。コンピュータが実数を有限のビット数で表現するため、多くの実数は正確に表現できず、丸め誤差が発生します。この誤差は計算を重ねるごとに蓄積される可能性があります。

<math.h> の関数も浮動小数点演算に基づいており、入力値の精度や内部アルゴリズムによって誤差を含む可能性があります。特に注意すべき点として以下があります。

• 丸め誤差: 計算の途中で発生する小数点以下の切り捨てや丸めによる誤差。
• 桁落ち: 非常に近い値同士の引き算によって、有効数字が大幅に失われる現象。
• 情報落ち: 桁数の大きく異なる値同士の加算や引き算で、絶対値の小さい数の情報が失われる現象。

これらの問題は <math.h> の関数自体の精度にも影響を与えます。例えば、sin(M_PI) は数学的には正確に0ですが、浮動小数点数の M_PI はπの近似値であるため、計算結果は厳密な0ではなく、非常に小さい浮動小数点数になるのが一般的です。

また、関数の定義域外の入力に対しては、前述のエラー（EDOM, ERANGE）が発生したり、NaNやINFが生成されたりします。これらの特殊な値が計算中に伝播していくと、最終的な結果がNaNになったり意図しない値になったりすることがあります。プログラム中でこれらの可能性を考慮し、必要に応じて分類関数でチェックしたり、エラー処理を実装したりすることが重要です。

10. まとめ

C言語の <math.h> ライブラリは、数値計算を行う上で非常に強力で便利なツール群を提供します。三角関数、指数・対数、べき乗、平方根といった基本的な関数から、浮動小数点数の詳細な操作や特殊な関数の計算まで、幅広い機能が利用可能です。

記事を通して、主要な関数の使い方、C99以降で追加された float/long double 版や型総称マクロ <tgmath.h>、エラー処理の仕組み、そして浮動小数点演算における精度に関する注意点などを解説しました。

<math.h> の関数を効果的に活用するためには、各関数の数学的な定義、引数の単位（ラジアンなど）、定義域と値域、そしてエラー発生時の挙動を理解することが重要です。また、浮動小数点数の性質上、厳密な等号比較は避け、許容誤差の範囲内での比較を行うなどの配慮も必要になります。

この記事が、C言語での数値計算における <math.h> の理解を深め、より正確で堅牢なプログラムを開発するための一助となれば幸いです。数学関数は多くの科学技術分野の基礎となるため、これらの機能を習得することは、C言語プログラマとしてのスキルを大きく向上させるでしょう。

---