Source

planets/jupiter/base.ts

import { JulianDay } from '@/types'
import { getDecimalYear } from '@/times'

// The value of K must be an integer
function getK (jd: JulianDay): number {
  const decimalYear = getDecimalYear(jd)
  const decimalK = 0.08430 * (decimalYear - 2011.20)
  return decimalK >= 0 ? Math.floor(decimalK) : Math.ceil(decimalK)
}

/**
 * Aphelion (time of passage at the closest point to the Sun)
 * @param {JulianDay} jd The julian day
 * @returns {JulianDay}
 * @memberof module:Jupiter
 */
export function getAphelion (jd: JulianDay): JulianDay {
  const kdash = getK(jd) + 0.5
  return 2455636.936 + 4332.897_065 * kdash
    + 0.000_1367 * Math.pow(kdash, 2)
}

/**
 * Perihelion (time of passage at the farthest point to the Sun)
 * @param {JulianDay} jd The julian day
 * @returns {JulianDay}
 * @memberof module:Jupiter
 */
export function getPerihelion (jd: JulianDay): JulianDay {
  const k = getK(jd)
  return 2455636.936 + 4332.897_065 * k + 0.000_1367 * Math.pow(k, 2)
}