The phase diagram of the Ising model in the presence of nearest- and next-nearest-neighbor interactions on a simple cubic lattice is studied within the framework of the differential operator technique. The Hamiltonian is solved by employing an effective-field theory with finite clusters consisting of a pair of spins. A functional form is also proposed for the free energy, similar to the Landau expansion, in order to obtain the phase diagram of the model. The transition from the ferromagnetic (or antiferromagnetic) phase to the disordered paramagnetic phase is of second order. On the other hand, a first-order transition is obtained from the lamellar phase to the paramagnetic phase, as well as from the lamellar phase to the ferromagnetic (or antiferromagnetic) phase, with the presence of a critical end point. An expected singular behavior of the first-order line at the critical end point is also obtained.