求助：矩阵乘法中参数太小导致乘出来的结果系统默认为0怎么办？

本人是学光学工程的，激光原理中的矩阵乘法需要大量矩阵连乘，而且这些矩阵都是秩为1的矩阵，但是一般光学参数都是特别小的，与1相比非常非常小，几近为0。这样乘出来的结果直接按单位矩阵给出的结果，这不是我想要的，我也查了format long，但这种形式也解决不了我的问题，求大神指点。大神可以看下，完整代码已经贴出了，很明显fenzi那一项按照原公式就不应该为无穷大，但因为ABCD1矩阵最后是个单位矩阵了（被matlab约等于的太多了）所以就发生错误了。