iseed = mod(8121*iseed+28411,134456);
ran(ii,jj)= iseed/134456;
還有就是全局變量iseed要在調用該函數之前就有定義和賦值
我試下來就沒問題了