仮想計算機構

IT業界と無縁な派遣社員のブログ

2021-06-01から1ヶ月間の記事一覧

【Python】ニュートンフラクタルの可視化

wikiの疑似コードを参考に実装します。 プログラム import numpy as np import cv2 Px,Py = 360,420 img = np.zeros((Px, Py,3), np.uint8) roots = [1+0j,-0.5+np.sqrt(3)/2j, -0.5-np.sqrt(3)/2j] r_hsv = cv2.cvtColor(np.uint8([[[255,0,0]]]),cv2.COLOR…

【Python】Lyapunovフラクタル

wikiの例を再現します。 繰り返し列:AB プログラム import numpy as np import cv2 b_hsv = cv2.cvtColor(np.uint8([[[255,0,0]]]),cv2.COLOR_RGB2HSV)[0][0] Px,Py = 360,360 s = [0,1] # AB S = [] # 周期列 M = int(Px*Py/len(s))+1 for i in range(M): S…

【Python】反復関数系の例

Sierpinskiの三角形 プログラム 実行結果 Barnsleyのシダ プログラム 実行結果 Koch曲線 プログラム 実行結果 Levy C曲線 プログラム 実行結果 参考 Sierpinskiの三角形 シェルピンスキーのギャスケット - Wikipedia プログラム import numpy as np import m…

Twitterでみかけた素数の一般項を計算してみる

「素数に一般項はない!」と主張する人をまれに見かけますが、それは嘘です。画像の式に好きな整数nを代入すると、n番目の素数が得られます。 pic.twitter.com/GF6hCJDVuF— 数学を愛する会 (@mathlava) June 6, 2021 面白そうだったので計算してみました。 …

PythonによるN体シミュレーション

はじめに N体シミュレーションを行いました。 初期値等は以下の情報を参考にしています。 Create Your Own N-body Simulation (With Python) | by Philip Mocz | The Startup | Medium nbody-python/nbody.py at master · pmocz/nbody-python · GitHub 上記…