A simple algorithm and a computational program to numerically compute the electric field gradient and the concomitant quadrupolar nuclear splitting is developed for an arbitrary ionic crystal. The calculations are performed using a point charge model. The program provides three different ways for the data input: by Bravais lattices, by lattice parameters, or by introducing any spatial structure. The program calculates the components of the electric field gradient, the asymmetry parameter and the quadrupolar splitting for a given number of nearest neighbors with respect to the nuclear charge as origin. In addition, the program allows the use of different Sternheimer antishielding factors.A simple algorithm and a computational program to numerically compute the electric field gradient and the concomitant quadrupolar nuclear splitting is developed for an arbitrary ionic crystal. The calculations are performed using a point charge model. The program provides three different ways for the data input: by Bravais lattices, by lattice parameters, or by introducing any spatial structure. The program calculates the components of the electric field gradient, the asymmetry parameter and the quadrupolar splitting for a given number of nearest neighbors with respect to the nuclear charge as origin. In addition, the program allows the use of different Sternheimer antishielding factors.