没仔细想修正方法。直接给你“比较快”的作法。
matlab 要快要向量化写法，我在 octave + 弱笔电执行我的程式后的输出为：
answer same, for = 1.06e+01 s, vec = 1.71e-02 s
Vectorization 619.79 x faster then for loop.
Original outer number is 10000, after mod() is 0
这程式假设你ao,bo,co 都是常数、不随位置改变。
以下程式
(注：其中的 clear -x ... 在 matlab 要改用 clearvars -except ...)

说一下，for n ... 只是为了计时而重复算10次

感谢!! 我想可能我太死板在用循环了，没有想到用向量化去处理