octave-連立一次方程式を解くについて

octaveで連立一次方程式を解く

p.mを以下のように組んで

 A = [
 41,0,16,13,0,0,0,2,0,0,2,0,20,2,0,8,1,3,1,0,0,0,0,0,0,0,0,2,1,2,0,5,0,3,0,0;
 31,0,4,1,1,2,0,0,0,0,1,1,8,0,0,6,0,5,6,0,0,0,0,0,0,0,0,39,0,4,1,38,0,0,0,0;
 35,0,2,2,0,1,0,0,0,0,2,0,2,0,0,12,0,1,3,0,0,0,0,0,4,0,0,6,0,2,0,10,0,4,7,7;
 ];
 
 b = {
 95.15;
 71.07;
 68.52;
 ];
 
 A\b

実行時は

 octave p.m

で終わり。

ベクトルの入力評価を出力させないためには

A=[....]; って　]の後に ;をつけると結果が表示されない。

余談csvから行列をつくるには以下のようにする。（もちっと最適化できるけど着陸がメインだったので・・・）

 cat score.csv | ruby -pe '$_.chop!;$_+=";\n"' > A.csv
 echo "A = [" > p.m ; cat A.csv >>p.m; echo "]"  >>p.m; echo >>p.m ; echo "b = ["  >>p.m; cat B.csv  >>p.m; echo "]"  >>p.m ; echo "A\b" >>p.m