// Calculates perihelion time of Earth from the sun. use sun.frink use planets.frink use binarySearch.frink mindist = 2 au maxdist = 0 au starttime = beginningOfYear[now[]] endtime = starttime + 365 days println["starttime is $starttime"] println["endtime is $endtime"] f = {|date| Planet.Earth.getCoordinates[date]@2} [mintime, mindist] = findALocalMinimum[f, starttime, endtime, 1 s] [maxtime, maxdist] = findALocalMaximum[f, starttime, endtime, 1 s] println["Minimum distance is " + (mindist -> "km") + ", " + (mindist -> "au") + ", " + (mindist-> "miles") + " at $mintime."] println["Maximum distance is " + (maxdist -> "km") + ", " + (maxdist -> "au") + ", " + (maxdist-> "miles") + " at $maxtime."]