Python 学習:連立方程式を解く

2015-11-03 :  PCクリニック
本文の前に、
-・・・ -・-
現時点での blogramのランクインカテゴリは、
6、2、1、0、1、 0、0、0、0、0(41)で、換算ポイント 84pt 。
「Firefox」が5位にアップ。「化学業界」「硝子業界」、
「e-radio」「グルコサミン」「FM COCOLO」はbg値変動のみ。
「C言語」「Perl」「Python」「FM青森」は全く変化無し。
・-・ - -・

さて、本文。

本ブログ内で、“Python solve”と検索すると、・・・
--->「検索結果

連立方程式を解く方法は、・・・・・

1年半以上前(2014-03-18)の記事:
Python 学習:Solve(Findroot)
でしたか。
“ scipy.optimize.broyden1”


でも他にもあった。

やはり1年半以上前(2014-03-26)の記事:
Python 学習:SymPy の学習
の続き???


例えば、

  2 * x**2 + 3 * y**2 = 10
  y = 4 * x - 1

の様な2元2次連立方程式を解くとしたら、・・・・・

from sympy import *

x, y = symbols( 'x y' )

eq1 = 2*x**2 + 3*y**2 - 10
eq2 = y - 4*x - 1

w = solve( [ eq1, eq2 ] , [ x, y ] )
print w

とすると解ける。


但し、答えは、・・・・・

  [(-6/25 + sqrt(494)/50, 1/25 + 2*sqrt(494)/25), ・・・]

厳密解


数値化するには、

  print N(w[0][0]), N(w[0][1])

とかするのかな。

 N() 関数または、 .evalf() メソッドを使う。

その結果は、

  0.204522215417857 1.81808886167143

となる。


念のために、検算?

  2 * 0.2**2 + 3 * 1.818**2 -----> 9.995372
1.818 - 4 * 0.2 -----> 1.018

ですね。


本日はここまで。


見ていただいた序でとは厚かましい限りですが、
お帰りに投票して頂けるとなお嬉しいです。 ⇒ blogram投票ボタン


151023
関連記事
スポンサーサイト

コメントの投稿

管理者にだけ表示を許可する

おきてがみ/blogram
blogram投票ボタン



おきてがみ

最新記事
カレンダー
09 | 2017/10 | 11
1 2 3 4 5 6 7
8 9 10 11 12 13 14
15 16 17 18 19 20 21
22 23 24 25 26 27 28
29 30 31 - - - -
月別アーカイブ
カテゴリ
最新コメント
検索フォーム
リンク
プロフィール

<紙>

Author:<紙>
ようこそ。
「パソコンヲタクの雑記帳」
もろもろなことを綴っています。
パソコン ヲタクってねくら?
画像は kami でなく kani です。

カウンター(fc2、i2i) /Google Analytics


i2i(from 2010-08-24)
Total =
Today  =  
Yesterday=
アンチエイジング

Google Analytics
ブックマーク