In the pioneering paper [13], the concept of Coherent Pair was introduced by Iserles et al. In particular, an algorithm to compute Fourier Coefficients in expansions of Sobolev orthogonal polynomials defined from coherent pairs of measures supported on an infinite subset of the real line is described. In this paper we extend such an algorithm in the framework of the so called Symmetric (1, 1)-Coherent Pairs presented in [8].