A method for accurate calculation of the harmonics generated bg a bipolar HVdc link is described. The method illustrates the importance of including detailed representation of the mutual coupling effects of de transmission lines, even when smoothing reactors are included. The overall solution is achieved bg means of a unified Newton algorithm in the harmonic domain.