This paper deals with time-harmonic problems with complex symmetric matrices in the electromagnetic field analysis. To solve such problems more stably and quickly, the Product-type Krylov Subspace Methods based on the COCG and COCR methods are proposed, and verified in finite element analyses of the time-harmonic eddy current and high-frequency electromagnetic field problems. As a result, in some cases, the number of iterations and the computational time have been reduced.