Calculation of zonal polynomials of 3×3 positive definite symmetric matrices