Compact stars, e.g., black holes and neutron stars, are the most energetic objects in astrophysics. These objects are accompanied by extremely strong gravity and a high velocity, which approaches the speed of light. Therefore, compact objects should be dealt with in Einstein’s relativity. This article will briefly introduce a numerical method that will allow us to obtain general solutions in general relativity. Several applications using numerical relativistic simulations will also be presented.