We present a non-iterative algorithm to reconstruct the isotropic acoustic wave speed from the measurement of the Neumann-to-Dirichlet map. The algorithm is designed based on the boundary control method and involves only computations that are stable. We prove the convergence of the algorithm and present its numerical implementation. The effectiveness of the algorithm is validated on both constant speed and variable speed, with full and partial boundary measurement as well as different levels of noise.