docs: def right_ascension_at_time(self,t):
          '''RA of prime meridian'''
          α0 = self.right_ascension_at_epoch
          t0 = self.epoch
          ω = self.rotational_speed
          return (α0 + ω * (t - t0)) % (2*π)