In a recent paper [PLB 837 (2023) 137669] we showed that all measured ($1s_\Lambda$, $1p_\Lambda$) pairs of $\Lambda$ binding energies in $\Lambda$-hypernuclei across the periodic table, $12\leq A \leq 208$, can be obtained from $\Lambda$-nucleus optical potential with only two adjustable $\Lambda N$ and NN$ parameters, associated leading linear quadratic terms nuclear density, derived by fitti...