Rで自作関数:gravity_by_latitude

# 緯度(度)を入力すると標準重力加速度(m/s^2)を返す関数
# 1980年国際重力式(International Gravity Formula 1980)に基づく
gravity_by_latitude <- function(lat_deg) {
  phi <- lat_deg * pi / 180 # 度 → ラジアン変換

  g <- 9.780327 * (1 + 0.0053024 * sin(phi)^2 - 0.0000058 * sin(2 * phi)^2)

  return(g)
}

# --- サンプル計算 ---
# 札幌市: 北緯 43.0621度
# 那覇市: 北緯 26.2124度

sapporo_lat <- 43.0621
naha_lat <- 26.2124

g_sapporo <- gravity_by_latitude(sapporo_lat)
g_naha <- gravity_by_latitude(naha_lat)

cat(sprintf("札幌市(緯度 %.4f°)の重力加速度: %.6f m/s^2\n", sapporo_lat, g_sapporo))
cat(sprintf("那覇市(緯度 %.4f°)の重力加速度: %.6f m/s^2\n", naha_lat, g_naha))