沒仔細想修正方法。直接給你「比較快」的作法。
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 ...)