Abstract This paper proposes a new algorithm for joint diagonalization of positive definite symmetric matrices and shows that it can be approximated by an independent component analysis problem. The good performance of the algorithm is illustrated with an example by compared with other approaches.