An efficient inverse method for 1.he supersonic/hypersonic axisymmetric body design is developed for the parabolized Navier-Stokes equations. The developed method is examined numerically for three extreme testcases in the supersonic(M_∞=3.0) and hypersonic(M_∞=6.28) speeds. The first one is a negative pressure distribution near a vacuum pressure and the second one is a positive pressure distribution over the whole region of the body. The last one is the case of abrupt change of pressure distribution to zero in the forward region of the body. These testcases show the robustness of the method. By introducing a regular-falsi method and by using a not-fully converged inverse solution, the convergence behavior was greatly improved.