A method of characteristic surfaces: The approximate solution of three-dimensional supersonic flow fields, from the near field up to the sonic boom region