あぅ、訂正。
s/かなり小さい/から想定されるよりかなり小さい/

あと、初歩的な質問であれなのだけれど、
ifc で適当な範囲の乱数えるのってどうすればいいのかな？
g77 だと a(i,j)=mod(irand(),100) てなかんじで、
適当に型変換して double な a(i,j) に代入できてたんだけど、
ifc だと a(i,j)=drand() 以外実行してもアドレスエラーとかおこられし、
drand() にしても乱数の系列を回してくれなくて、
１つの種に対して１つの乱数しか生成してくれない。
普段は C のひとなのでこのへんの作法がよくわからないっす。