计算各向异性介质的刚度矩阵，通过自洽理论来实现。-Computing the stiffness matrix and elements for TI by self-consistent method.The origin code is for VTI, so we extend the code to HTI and VTI.Interchanging the matrix and inclusion will give you stiff and soft estimations.