Gaussian kernel density estimation of the spin-axis orientations computed using the SciPy library (Virtanen et al. 2020) for a very prolate case (p = −6.7). Top panel: (ɛ_c, ϕ_c) map. The results in ecliptic coordinates, (λ_p, β_p), are shown in the middle panel and those in equatorial coordinates, (α_p, δ_p), are displayed in the bottom panel that includes the actual observations (in red) of interstellar object 1I/2017 U1 (‘Oumuamua) available from the MPC.