« Python+numbaでz^2+cのジュリア集合を色んなcで描く。 | トップページ | 5/21のGoogleトップページはアフリカの楽器ムビラ!自由に弾けるのでやってみたが、、、学校のチャイムしか思いつかなかった!→というのをやってみたが、Windows10の画面キャプチャ、めちゃくちゃノイジーで音が割れる(GameDVR) »

2020年5月21日 (木)

Python+Sympyで全部係数を変数にした連立3元一次方程式を計算。init_printingを使うと表示が綺麗なのと、解が辞書型なのに[’x’]とかで取り出せない。[x]じゃないと。

手でもできる計算ではありますが、とてもめんどくさい

a-b*x1*y1-c*x1+y1=0
a-b*x2*y2-c*x2+y2=0
a-b*x3*y3-c*x3+y3=0

という計算をする必要が出てきた。a,b,cが未知数で、xi,yiは与えられてるデータ。

ここはSympyで計算してみよう。var('.....')で記号と定義して、solveで解く。一瞬、、、だけど表示がテキストでイマイチ。単なるprintを使ってるから、、、

Sympy001

pprint(pretty print)を使うといいと聞いたことがあるのでやってみると、

Sympy002

さっきよりはいいけどアスキーアート感が半端ない。

ここはinit_printing()をつかうとlatex表示するらしいということでやってみた。

表示はdisplayを使う。

Sympy003

おお、キレイにでた。が、長すぎて切れてる。。。

未知数ごとに取り出そうとしてs['a']とかするとエラー。

なんで?辞書型なのに?と思ったら、

s[a]とかしないとダメらしい。やってみると

 

Sympy004

うまくいった。

もっと簡単にならないかとsimplifyを試すと?

5sympy004

これ以上は単純にはならないか。

 

« Python+numbaでz^2+cのジュリア集合を色んなcで描く。 | トップページ | 5/21のGoogleトップページはアフリカの楽器ムビラ!自由に弾けるのでやってみたが、、、学校のチャイムしか思いつかなかった!→というのをやってみたが、Windows10の画面キャプチャ、めちゃくちゃノイジーで音が割れる(GameDVR) »

日記・コラム・つぶやき」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

« Python+numbaでz^2+cのジュリア集合を色んなcで描く。 | トップページ | 5/21のGoogleトップページはアフリカの楽器ムビラ!自由に弾けるのでやってみたが、、、学校のチャイムしか思いつかなかった!→というのをやってみたが、Windows10の画面キャプチャ、めちゃくちゃノイジーで音が割れる(GameDVR) »

最近の記事

2021年2月
  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            
フォト
無料ブログはココログ