Fast Computation of Orthogonal Systems with a Skew‐Symmetric Differentiation Matrix