We present the Elliptica code for generation of initial data in simulated binary neutron star mergers. The relevant theoretical models of numerical relativity are described, and we demonstrate the generation of initial data for a simple binary neutron star system. Further, we analyze the accuracy of the numerical model it converges exponentially to the correct solution of Einstein's equations with respect to the computational grid size.