ni = max(2 * int((xmax - xmin) / hs_e) + 3, 5)