We calculate the two-point correlation function <x(t2)x(t1)> for a subdiffusive continuous time random walk in a parabolic potential, generalizing well-known results for the single-time statistics to two times. A closed analytical expression is found for initial equilibrium, revealing a clear deviation from a Mittag-Leffler decay.