在Python中，可以使用scipy庫中的quad函數來進行數值驗證。quad函數用于對給定的函數進行數值積分，可以通過比較得到的數值積分結果與預期的解析解進行驗證。
以下是一個使用quad函數進行數值驗證的示例代碼：
from scipy.integrate import quad
# 定義要積分的函數
def integrand(x):
return x**2
# 使用quad函數進行數值積分
result, error = quad(integrand, 0, 1)
# 打印數值積分結果和誤差
print(f"數值積分結果: {result}")
print(f"誤差: {error}")
# 預期的解析解
expected_result = 1/3
# 比較數值積分結果與預期解析解
if abs(result - expected_result) < 1e-6:
print("數值驗證通過")
else:
print("數值驗證未通過")
在這個示例中，我們定義了一個要積分的函數integrand(x)，然后使用quad函數對該函數在區間[0, 1]上進行數值積分。我們將得到的數值積分結果與預期的解析解進行比較，如果它們的差異在一個很小的誤差范圍內，則認為數值驗證通過。
