Circular harmonics expansion method