An algebraic, calculational, non‐fitting method of parametrization of the cubic and hexagonal crystal‐field (CF) Hamiltonians false(HCFfalse) with only two independent CF parameters B40 and B60 is proposed. It is based on the linear relation σ2false(Jfalse)=∑kAk2false(Jfalse)Sk2 between σ2false(Jfalse) the square of the second moment of the electron state false|J〉 CF splitting and the squares of the multipolar CF strengths Sk2=12k+1∑qfalse|Bkq|2, where AkJ=〈Jfalse|false|Cfalse(kfalse)false|false|J〉. For a pair of false|Ji〉 and false|Jj〉 electron states an algebraic solution for S42 and S62, and hence for the B40 and B60 CF parameters can be gained. The necessary conditions for physical correctness of the solutions false(Sk2≥0false) impose limitations on the observed σ2false(Jifalse)/σ2false(Jjfalse) ratios depending on the respective ratios Ak2false(Jifalse)/Ak2false(Jjfalse). These restrictions can be used to verify postulated CF splitting diagrams. In the high‐symmetry crystal fields, the electron states with J=2 or 5/2 undergo splitting only by the k=4 HCF multipole, and from their CF splitting second moments the B40 parameter can be directly obtained. For J=2, false|B40false|=6310Δfalse(Jfalse)false|A4false(Jfalse)false|, whereas for J=5/2, false|B40false|=7Δfalse(Jfalse)false|A4false(Jfalse)false|, where Δ is the energy gap between the two Stark's levels of the false|J〉 state. The proposed method is demonstrated and discussed for several systems of RE3+ ions in some high‐symmetry crystal matrices.