mathematica的数值计算精度

txrtxs

最近在做一个有限元分析，发现用mathematica编的程序生成的单元质量和刚度矩阵没有对称，具体如下：
单元刚度矩阵eleKK=Transpose[bK1]. bK2. bK1
其中bK1为8乘240阶，bK2为8乘8阶矩阵。
bK2经验证为对称矩阵，从单元刚度矩阵表达式分析，应该为对称，但eleKK==Transpose[eleKK]为False，很奇怪，大家有没有遇到这种情况，急问！！！！！
附件后缀名改为.nb即可运行。