Notebook
これは日々の作業を通して学んだことや毎日の生活で気づいたことをを記録しておく備忘録である。
HTML ファイル生成日時: 2025/08/15 15:56:28.294 (台灣標準時)
sympy パッケージ
sympy パッケージというものがあって、方程式を解析的に解かせる、というよ
うなことができるようなので、試してみたでござる。
ガウス積分
以下のようにすればガウス積分を計算してもらうことができるようでござる。
#!/usr/pkg/bin/python3.9
# Time-stamp: <2022/04/18 21:01:15 (CST) daisuke>
# importing sympy module
import sympy
# variable
x = sympy.Symbol ('x')
# function f(x)
f = sympy.exp (-x**2)
# integration of f(x)
I = sympy.integrate (f, (x, -sympy.oo, sympy.oo))
# printing result
print (I)
実行結果は以下の通りにござる。
% ./test_00.py
sqrt(pi)
式の展開
式の展開の例。
#!/usr/pkg/bin/python3.9
# Time-stamp: <2022/04/18 21:01:20 (CST) daisuke>
# importing sympy module
import sympy
# variable
x = sympy.Symbol ('x')
# function f(x)
f = (x+1) * (x-1)
# expanding f(x)
f2 = sympy.expand (f)
# printing result
print (f, "=", f2)
実行結果は以下の通りにござる。
% ./test_01.py
(x - 1)*(x + 1) = x**2 - 1
極限
極限の求め方。
#!/usr/pkg/bin/python3.9
# Time-stamp: <2022/04/18 21:01:25 (CST) daisuke>
# importing sympy module
import sympy
# variable
x = sympy.Symbol ('x')
# function f(x)
f = sympy.sin (x) / x
# lim (x -> 0) sin(x)/x
lim_0 = sympy.limit (f, x, 0)
# printing result
print ("lim_0 sin(x)/x =", lim_0)
実行結果は以下の通り。
% ./test_02.py
lim_0 sin(x)/x = 1
代数方程式を解く
二次方程式を解く。
#!/usr/pkg/bin/python3.9
# Time-stamp: <2022/04/18 21:01:29 (CST) daisuke>
# importing sympy module
import sympy
# variable
x, a, b, c = sympy.symbols ('x a b c')
# function f(x)
f = a*x**2 + b*x + c
# solving f(x) = 0
sol = sympy.solve (f, x)
# printing result
print ("x =", sol)
実行結果は以下の通り。
% ./test_03.py
x = [(-b + sqrt(-4*a*c + b**2))/(2*a), -(b + sqrt(-4*a*c + b**2))/(2*a)]
微分方程式を解く
微分方程式を解く。
#!/usr/pkg/bin/python3.9
# Time-stamp: <2022/04/18 21:01:34 (CST) daisuke>
# importing sympy module
import sympy
# variable
x, y, k = sympy.symbols ('x y k')
# function y(x)
y = sympy.Function ('y')
# solving dy/dx = -ky
sol = sympy.dsolve (sympy.Eq (y(x).diff (x), -k*y(x)), y(x))
# printing result
print ("y =", sol)
実行結果は以下の通り。
% ./test_04.py
y = Eq(y(x), C1*exp(-k*x))
式の整理
式の整理の例。
#!/usr/pkg/bin/python3.9
# Time-stamp: <2022/04/18 21:01:38 (CST) daisuke>
# importing sympy module
import sympy
# variable
x, y = sympy.symbols ('x y')
# function f(x,y)
f = (x+y)**2 - (x-y)**2
# solving f(x) = 0
f2 = sympy.simplify (f)
# printing result
print (f, "=", f2)
実行結果。
% ./test_05.py
-(x - y)**2 + (x + y)**2 = 4*x*y
n=0 の Lane-Emden 方程式を解く
n=0 の Lane-Emden 方程式を解いてみる。
#!/usr/pkg/bin/python3.9
# Time-stamp: <2022/04/18 21:08:41 (CST) daisuke>
# importing sympy module
import sympy
# variable
x = sympy.Symbol ('x')
# function y(x)
y = sympy.Function ('y')
# differential equation
lane_emden_0 = sympy.Eq (x**-2 * (x**2 * y(x).diff (x)).diff (x), -1)
# solving lane-emden equation of n=0
sol = sympy.dsolve (lane_emden_0, y(x))
# printing result
print ("y =", sol)
実行結果は以下の通り。
% ./test_06.py
y = Eq(y(x), C1 + C2/x - x**2/6)
- About this article:
- author: daisuke
- file: 20220418_0.html
- category: Computer___Python
- title: sympy パッケージ
- mode: public
- last modified: 2022/04/18 21:28:31 (Taiwan Standard Time)
- html generated: 2025/08/15 15:56:28.294 (Taiwan Standard Time)
Frequently accessed files
- Misc___Taiwan/20240207_00.html
- 149 page views
- title: 台灣から台灣の外に EMS で荷物を発送する方法
- Misc___Taiwan/20240819_00.html
- 129 page views
- title: 住所から台灣の郵便番号を調べる方法
- Computer___NetBSD/20230119_00.html
- 123 page views
- title: NetBSD でバイナリーパッケージを利用する方法
- Computer___TeX/20231107_00.html
- 110 page views
- title: LaTeX での counter の利用方法について
- Computer___Network/20230516_00.html
- 104 page views
- title: OpenVPN 2.6 を使い VPN Gate に接続するときの注意点
- Misc___Taiwan/20240903_01.html
- 85 page views
- title: 台湾の郵便局で EMS を利用して荷物を海外に送る方法
- Computer___Network/20230726_00.html
- 78 page views
- title: git の SSL certificate problem の解決方法
- Computer___NetBSD/20250301_01.html
- 73 page views
- title: yt-dlp で YouTube の動画をダウンロードするときのこと
- Computer___Network/20241214_00.html
- 63 page views
- title: Google Colab で Python 3.12 を使う方法
- Computer___NetBSD/20240805_03.html
- 62 page views
- title: NetBSD 10 のインストール
- Computer___Network/20240130_00.html
- 59 page views
- title: Google Colaboratory で Python 3.12 を使う方法
- Computer___TeX/20240414_00.html
- 55 page views
- title: LuaTeX での fontspec を使ったフォントの設定の方法
- Food___Taiwan/20230609_06.html
- 55 page views
- title: 「越芳越南牛肉河粉」のフォー
- Computer___TeX/20230726_01.html
- 52 page views
- title: Beamer の Metropolis テーマで block の色を変える方法
- Misc___Taiwan/20250728_01.html
- 51 page views
- title: 「我那小小多山的國家」という言葉が流行っている
- Food___Taiwan/20230608_04.html
- 50 page views
- title: 「田園美食屋」の咖哩煎雞腿
- Computer___TeX/20240414_01.html
- 50 page views
- title: LuaTeX を使って日本語と中国語の漢字が混在した PDF ファイルを作成する方法
- Computer___Python/20240101_00.html
- 49 page views
- title: Matplotlib の 3D plot においての注意点
- Computer___FreeBSD/20220621_0.html
- 49 page views
- title: FreeBSD での X.org の設定の仕方
- Misc___Japan/20240610_00.html
- 47 page views
- title: NHK ラジオの「聞き逃し」の変更点 (2024 年 06 月)
- Computer___NetBSD/20250409_00.html
- 45 page views
- title: Raspberry Pi 5 に NetBSD-10.1 をインストールする手順
- Computer___NetBSD/20241102_00.html
- 44 page views
- title: ImageMagick の convert コマンドについて
- Science___Math/20220420_0.html
- 44 page views
- title: ラプラシアンの三次元極座標表示
- Computer___TeX/20240410_00.html
- 43 page views
- title: pdfTeX を使って PDF ファイルを作る方法
- Food___Taiwan/20230711_15.html
- 43 page views
- title: 「東方美早餐店」での朝ご飯 (2023 年 06 月中旬)
- Computer___TeX/20240411_00.html
- 41 page views
- title: LuaTeX を使って PDF ファイルを作成する方法
- Computer___NetBSD/20240810_00.html
- 40 page views
- title: Raspberry Pi 5 に NetBSD-10 をインストールしてみた
- Computer___NetBSD/20250728_00.html
- 37 page views
- title: pkgsrc の geography/gdal-lib について
- Computer___Network/20230508_00.html
- 36 page views
- title: git push するときにパスワードの入力を省略する方法
- Misc___Taiwan/20220831_4.html
- 36 page views
- title: 台南の民權派出所
HTML file generated by Kinoshita Daisuke.
