It is a numerical method for solving an inverse problem about the shape design of aerodynamic body in inviscid subsonic flows. This method transfers the original Euler equations into the stream-function plane, where it solves the Riemann problem across the streamline presenting the solid-wall and obtains the geometry of the solid-wall concurrently.