The most important sources for laser-interferometric gravitational-wave detectors like LIGO or VIRGO are catastrophic events such as coalescence of a neutron-star binary. The final phase, or the last three milliseconds, of coalescence is considered. We describe results of numerical simulations of coalescing binary neutron stars using Newtonian and post-Newtonian hydrodynamics code and then discuss recent development of our 3D GR code.