数値計算プログラム

連立方程式は行列の解になります。
一番オーソドックスなのが、「ガウスのはきだし法」です。

dim a(n,n+1)
i:=1
j:=1
repeat i
repeat j
input a(i,j)
until j>=n+1
until i>=n
k:=1
repeat k
w:=a(k,k)
j:=k
repeat j
a(k,j)*=a(k,j)/w
until j>=n+1
i:=1
repeat i
w:=a(i,k)
j:=k
repeat j
a(i,j):=a(i,j)-w*a(k,j)
until j>=n+1
until i>=n or i:=k
until k>=n or w

不定の処理は省いています