In this paper, we present a probabilistic algorithm to compute the coefficients of modular forms of level one. Focusing on the Ramanujan's tau function, we give the explicit complexity of the algorithm. From a practical viewpoint, the algorithm is particularly well suited for implementations.