[[練習問題]] ガウスの掃き出し法を使っています。 # coding: utf-8 a = raw_input('1次の連立方程式を、例えばax+by=c dx+ey=fであれば、a,b,c|d,e,fと入力してください。') a = a.split('|') #入力された文字列を式単位に分割してリスト型に変換します。 for h in range(len(a)): #式を係数単位で分割してリスト型に変換します。 a[h] = a[h].split(',') for g in range(len(a[h])): #係数をfloat型に変換します。 a[h][g] = float(a[h][g]) for i in range(len(a)): if a[i][i] == 0: #行列の(i,i)成分が0である場合、i以降の行で、i成分が0でないj行とi行を入れ替えます。 for j in range(i,len(a)): if a[j][i] != 0: a[i],a[j] = a[j],a[i] break else: pass tpv = float(a[i][i]) #i行の対角成分をtpvとします。 for k in range(len(a[i])): #i行を対角成分で割ります。 a[i][k] = a[i][k] / tpv for l in range(i): #第1行からi-1行までのループ処理です。 rat = a[l][i] #比率を定めます。 for m in range(len(a[l])): #i行を(l,i)倍して、l行から引きます。 a[l][m] = a[l][m] - rat * a[i][m] for n in range(i+1,len(a)): #第i+1行から最後の行までのループ処理です。 rat = a[n][i] for o in range(len(a[n])): a[n][o] = a[n][o] - rat * a[i][o] print(a)