ダメ出し：妙な数値表記

湿度諸量の計算法（まとめ）　でみかけた

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 の方が多項式を計算するための定石手段）。