def mCVTC(molFrac, vTC, vV, fW, eps=1.0):
    """
    mCVTC(molFrac, vTC, vV, fW, eps=1.0)
    
    multi-Component Vapor Thermal Conductivity in W/m/K

    Parameters:
        molFrac, list of mol fractions
        vTC, list of pure component vapor thermal conductivities (components 
             ordered as in molFrac, fW, and vV) 
        vV, list of pure component vapor viscosities (components ordered as
            in molFrac, fW, and vTC)
        fW, list of formula weights for each component in molFrac
        eps, epsilon - typical value is 1.0
    
    Returns:
        multi-component vapor thermal conductivity in W/m/K
    """
    tCm = 0
    Aij = []
    nMF = len(molFrac)
    
    for idx in range(nMF):
        Aij.append([])
        sumyjAij = 0
        for jdx in range(nMF):
            fWR=fW[idx]/fW[jdx]
            vR = vV[idx]/vV[jdx]
            if idx == jdx:
                sumyjAij += molFrac[jdx]
            elif idx < jdx:
                Aij[idx].append(eps*(1+vR**(1./2)*(1/fWR)**(1./4))**2/(8*(1+fWR))**(1./2))
                sumyjAij += molFrac[jdx]*Aij[idx][jdx-(1+idx)]
            else:
                sumyjAij += molFrac[jdx]*Aij[jdx][idx-(1+jdx)]*vR/fWR
        tCm += molFrac[idx]*vTC[idx]/sumyjAij
    
    return tCm