Code to calculate correspondece between fractional anisotropy and Lambda_th