一段Zemax程序，和Matlab/C+相似，好心人帮忙看看是干嘛的

做matlab程序的一个原则是在复杂度不变的情况下，尽量用矩阵运算代替循环运算以节省时间。这段程序中C12和C13的求法都可以转化为矩阵计算，但是C11由于对矩阵的局部做了循环相乘，这里很难改成对整个矩阵的运算。
所以最好是能修改一下算法，降低复杂度。或者避免矩阵的局部循环相乘。