超雑

連立方程式を計算してみました。
function sim(e_a: Array): Array{ var i, j ,k, l, m; var e_c = new Array(); var e_s = new Array(); var sol = new Array(); var len = e_a[0].length - 1; for(i = 0; i < len; i ++){ e_c[i] = e_a[i].concat(); e_s[i] = e_a[i].concat(); } //0除算回避 var cha = false; var temp; var zer = new Array(); for(i = 0; i < len; i ++){ zer[i] = new Array(); for(j = 0; j < len; j ++){ if(e_c[i][j] == 0){ zer[i].push(j); } } if(zer[i].length == 0){ zer[i][0] = -1; } } for(i = 0; i < len; i ++){ for(j = 0; j < zer[i].length; j ++){ if(zer[i][j] == i){ cha = false; for(k = 0; k < len; k ++){ for(l = 0; l < zer[k].length; l ++){ for(m = 0; m < zer[i].length; m ++){ if(zer[k][l] == i || zer[i][m] == k){ break; }else if(l == zer[k].length - 1 && m == zer[i].length - 1){ temp = e_s[i].concat(); e_s[i] = e_s[k].concat(); e_s[k] = temp.concat(); temp = zer[i].concat(); zer[i] = zer[k].concat(); zer[k] = temp.concat(); cha = true; } } if(cha){ break; } } if(cha){ break; } } } } } for(i = 0; i < len; i ++){ e_c[i] = e_s[i].concat(); } //ガウス・ジョルダン法 for(i = 0; i < len; i ++){ for(j = 0; j < len + 1; j ++){ e_s[i][j] = e_c[i][j] / e_c[i][i]; } for(j = 0; j < len; j ++){ e_c[j] = e_s[j].concat(); } for(j = 0; j < len; j ++){ for(k = 0; k < len + 1; k ++){ if(j != i){ e_s[j][k] = e_c[j][k] - e_c[j][i] * e_c[i][k]; } } } for(j = 0; j < len; j ++){ e_c[j] = e_s[j].concat(); } } for(i = 0; i < len; i ++){ sol[i] = e_c[i][len]; } return sol; }
0除算防止の部分がかなり雑です。


使い方:
引数は係数を入れた配列を用意して、
それのさらに配列。
戻り値は配列。


例:
ax + by + cz = d, ex + fy + gz = h, ix + jy + kz = l の場合
引数は[[a, b, c, d], [e, f, g, h], [i, j, k, l]]
戻り値は[α, β, γ, δ]