import spharpy

n_max = 2
sampling = spharpy.samplings.equal_area(0, n_points=500)

Y_nm = spharpy.SphericalHarmonics(
    n_max=n_max,
    coordinates=sampling,
    basis_type='complex',
    normalization='N3D',
    condon_shortley='auto'
).basis

axs, gs = plot_basis_functions(Y_nm, sampling)