湿度諸量の計算法（まとめ） でみかけた
Okada <- function(t) {
## 入力：気温(℃)
## 出力：飽和水蒸気圧(hPa)
exp(1.809378 +
0.07266115 * t +
(-3.003879 * 10^-4) * t^2 +
(1.181765 * 10^-6) * t^3 +
(-3.863083 * 10^-9) * t^4)
}
Okada2 <- function(t) {
## 入力：気温(℃)
## 出力：飽和水蒸気圧(hPa)
exp(1.809378 +
0.07266115 * t -
3.003879e-4 * t^2 +
1.181765e-6 * t^3 -
3.863083e-9 * t^4)
}
Okada3 <- function(t) {
## 入力：気温(℃)
## 出力：飽和水蒸気圧(hPa)
exp((((-3.863083e-9 * t + 1.181765e-6) * t -3.003879e-4) * t + 0.07266115) * t + 1.809378)
}
> options(digits=16)
> Okada(30)
[1] 42.42093319229595
> Okada2(30)
[1] 42.42093319229595
> Okada3(30)
[1] 42.42093319229597
Okada2 と Okada3 では，計算量が違う（Okada3 の方が多項式を計算するための定石手段）。
