An $(\epsilon,\delta)$-accurate level set estimation with a stopping criterion