code: def fundamental_arguments(t):
  
      
      a = fa4 * t
      a += fa3
      a *= t
      a += fa2
      a *= t
      a += fa1
      a *= t
      a += fa0
      fmod(a, ASEC360, out=a)
      a *= ASEC2RAD
      if getattr(t, 'shape', ()):
          return a
      return a[:,0]
docs: Compute the fundamental arguments (mean elements) of Sun and Moon.
  
      `t` - TDB time in Julian centuries since J2000.0, as float or NumPy array
  
      Outputs fundamental arguments, in radians:
            a[0] = l (mean anomaly of the Moon)
            a[1] = l' (mean anomaly of the Sun)
            a[2] = F (mean argument of the latitude of the Moon)
            a[3] = D (mean elongation of the Moon from the Sun)
            a[4] = Omega (mean longitude of the Moon's ascending node);
                   from Simon section 3.4(b.3),
                   precession = 5028.8200 arcsec/cy)