The formulations of the finite-field approach to calculate the linear and non-linear optical coefficients (i, (ij, (ijk and (ijkl of a molecular system with different symmetries have been deduced and summarized. The possible choices of the energy sets of the 48 frequent point groups have been optimized and categorized into 11 classes. With the restriction of symmetry operators, a minimum of 9, no more than 21 energy points have to be calculated in order to determine the coefficients, except in the case of the first class to which C1 point group belongs and in which the 34 non-relative energy points selected in our uniform and general scheme are all needed. The symmetric operators that cause some of the tensor components to vanish have been demonstrated as well. The formulations of the finite-field approach to calculate the linear and non-linear optical coefficients (i, (ij, (ijk and (ijkl of a molecular system with different symmetries have been deduced and summarized. The possible choices of the energy sets of the 48 frequent point groups have been optimized and categorized into 11 classes. With the restriction of symmetry operators, a minimum of 9, no more than 21 energy points have to be calculated in order to determine the coefficients, except in the case of the first class to which C1 point group belongs and in which the 34 non-relative energy points selected in our uniform and general scheme are all needed. The symmetric operators that cause some of the tensor components to vanish have been demonstrated as well.