code: def julian_day(year, n):
      
      if n > 59 and isleap(year):
          n += 1
      return datetime(year, 1, 1) + timedelta(n - 1)
docs: Return the Julian day n(1 <= n <= 365).
  
      Leap days are not counted; that is, in all years -- including leap years --
      February 28 is day 59 and March 1 is day 60. It is impossible to explicitly
      refer to the occasional February 29.