練習問題/解答例/連立方程式を解く/Python
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
[[練習問題]]
ガウスの掃き出し法を使っています。
# coding: utf-8
a = raw_input('1次の連立方程式を、例えばax+by=c dx+ey=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以降の...
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)
終了行:
[[練習問題]]
ガウスの掃き出し法を使っています。
# coding: utf-8
a = raw_input('1次の連立方程式を、例えばax+by=c dx+ey=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以降の...
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)
ページ名:
