We develop a Monte Carlo algorithm named the Frozen Gaussian Sampling (FGS) to solve the semiclassical Schrodinger equation based on the frozen Gaussian approximation. Due to the highly oscillatory structure of the wave function,traditional mesh-based algorithms suffer from “the curse of dimensionality”, which gives rise to more severe computational burden when the semiclassical parameter issmall. The Frozen Gaussian sampling outperforms the existing algorithms in that it is mesh-free in computing the physical observables and is suitable for high dimensional problems. We also discussion the extension of the FGS approach to the mixed quantum-classical dynamical models.