Matlabのfsolveの代替としてcsolveがおすすめなようです。 http://sims.princeton.edu/yftp/optimize/mfiles/ ただし引数は同じではないです。csolveの引数は…

九州大学の滑降シンプレックス法のアルゴリズムをフローチャートにしてみました。

演算Octave 和A+B 差A-B 行列の乗算A*B 各要素間の積 A.*B 各要素間の除算A./B 各要素間の指数算A.^B 転置*A’ Aの(m,n)要素A(m,n) AのN番目の列を選択A(N,:) AのN番目の…

fortranとは直接の関係はないですが、数値計算をする上での乱数の記事でとても参考になったのは、 良い乱数・悪い乱数です。 fortranからCの関数を呼ぶサンプルのプログラムは Physics 410: Comput…

Fortran90は固定形式では1行72文字(桁)まで、自由形式では1行132文字(桁)までの制限があります。 これはコンパイル時のオプションで無効にできます。 具体的にはコンパイル時のオプションで -ffree-lin…

Fortran90でコマンドラインから引数を読み込むやり方を調べてみました。 C-Fortran Interfaceの”Command Line Arguments”をfortran90化にしたも…