Summary  
This chapter covers how to perform numerical integration (definite integrals) and solve ordinary differential equations numerically using SciPy’s integrate module.

General domain of usage  
Scientific computing

数値積分は、曲線下の面積を計算したり、解析的な解が存在しない方程式を解いたりするための手法。科学技術計算では、定積分の評価や常微分方程式（ODE）の解法が必要となる場面が多く、正確な解が不明または非常に複雑な場合がある。SciPyの **scipy.integrate** モジュールは、これらのタスクに対して強力かつ使いやすいツールを提供し、わずか数行の `python` コードで数値積分やODEの数値解法を実現可能。

from scipy import integrate
import numpy as np

# Define the function to integrate
def f(x):
    return np.sin(x)

# Compute the definite integral of sin(x) from 0 to pi
result, error = integrate.quad(f, 0, np.pi)
print("Integral result:", result)
print("Estimated error:", error)

from scipy.integrate import solve_ivp
import numpy as np
import matplotlib.pyplot as plt

# Define the ODE: dy/dt = -2y
def dydt(t, y):
    return -2 * y

# Initial condition
y0 = [1]

# Time span for the solution
t_span = (0, 5)

# Solve the ODE
solution = solve_ivp(dydt, t_span, y0, t_eval=np.linspace(0, 5, 100))

# Plot the solution
plt.plot(solution.t, solution.y[0])
plt.xlabel("t")
plt.ylabel("y(t)")
plt.title("Solution of dy/dt = -2y with y(0) = 1")
plt.show()

`scipy.integrate.quad` を使用すると、関数は積分の計算値と誤差の推定値の両方を返します。上記の例では、`sin(x)` を 0 から π まで積分すると、2 に非常に近い結果が得られ、これは厳密な解析解と一致します。これにより、数値積分ルーチンの精度と信頼性の両方が示されます。

常微分方程式の場合、`scipy.integrate.solve_ivp` は指定された区間で解を計算します。ODE の例では、初期条件 `dy/dt = -2y` を持つ方程式 `y(0) = 1` は指数関数的減衰を表します。解は `y` が時間とともに滑らかに減少する様子を示し、簡単なプロットでこれを可視化できます。出力は期待される解析解 `y(t) = exp(-2t)` と一致します。

SciPy で定積分に使用される関数はどれですか？

`scipy.integrate.solve_ivp` は何を解きますか？

なぜ数値積分は科学技術計算において重要なのですか？

強力な数学的バックグラウンドと中級から上級のPythonスキルを持つ学生を対象に、SciPyを活用した科学計算、最適化、積分、高度な数学的操作に焦点を当てた包括的なコースです。

科学計算のためのSciPyの構造、主要モジュール、および基本機能を探求します。

SciPyの強力な線形代数機能を活用し、実世界の数学的問題を解決します。

科学および工学の問題に対して、SciPy を用いた最適化手法と根探索アルゴリズムを習得します。

SciPyを使用した積分、補間、基本的な信号処理を含む高度な数学的操作を学びます。

Scipy.integrateによる数値積分

1. SciPy で定積分に使用される関数はどれですか？

2. `scipy.integrate.solve_ivp` は何を解きますか？

3. なぜ数値積分は科学技術計算において重要なのですか？

Scipy.integrateによる数値積分

1. SciPy で定積分に使用される関数はどれですか？

2. scipy.integrate.solve_ivp は何を解きますか？

3. なぜ数値積分は科学技術計算において重要なのですか？

2. `scipy.integrate.solve_ivp` は何を解きますか？