We consider an inverse problem of recovering all spatial dependent coefficients in the time dependent Schrödinger equation defined on an open bounded domain in R n , n≥2, with smooth enough boundary. We show that by appropriately selecting a finite number of initial conditions and a fixed Dirichlet boundary condition, we may recover all the coefficients in a Lipschitz stable fashion from the corresponding finitely many boundary measurements made on a portion of the boundary. The proof is based on a direct approach, which was introduced in [11], to derive the stability estimate directly from the Carleman estimates without any cut-off procedure or compactness-uniqueness argument.