Abstract. We discuss a method for computing p≤x 1/p, using time about x 2/3 and space about x 1/3 . It is based on the Meissel-Lehmer algorithm for computing the prime-counting function π(x), which was adapted and improved by Lagarias, Miller, and Odlyzko. We used this algorithm to determine the first point at which the prime harmonic sum first crosses 4.