2021-07-09から1日間の記事一覧

Python scipyで任意の方程式の解を求める

from scipy.optimize import fsolve import numpy as np # 変数 beta = 0.6955076793404303 ionpair = 4886 # 方程式 func = lambda delE : np.sqrt(delE/(np.log(ionpair*beta**2)-np.log(1-beta**2)-beta**2))*beta-92 # 初期値 initial_guess = 1000 # 解…