Python中文
首页
教程
问答
标签
搜索
登录
注册
求数值积分的根
回答此问题可获得
20
贡献值,回答如果被采纳可获得
50
分。
<p>我正在尝试用Python复制这个Mathematica程序:</p> <p><img src="https://i.stack.imgur.com/0Km6p.png" alt="Mathematica program"/></p> <p>它找到数值积分的根,并形成这些数值的图。然而,我无法尝试逃跑。在</p> <p>当前尝试:</p> <p>从整合导入四边形 从scipy导入集成 从scipy.optimize公司导入F溶剂 pylab进口 将numpy作为np导入</p> <pre><code># Variables. boltzmann_const = 1.38e-23 planck_const = 6.62e-34 hbar = planck_const / ( 2 * np.pi ) transition_temp = 9.2 gap_energy_at_zero_kelvin = 3.528 / ( 2 * transition_temp * boltzmann_const ) debye_freq = ( 296 * boltzmann_const ) / hbar # For subtracting from root_of_integral a_const = np.log( ( 1.13 * hbar * debye_freq ) / ( boltzmann_const * transition_temp) ) # For simplifying function f. b_const = ( hbar * debye_freq ) / ( 2 * boltzmann_const) def f( coherence_length, temp ): # Defines the equation whose integral will have its roots found. Epsilon = coherence length. Delta = Gap energy. squareRoot = np.sqrt( coherence_length*coherence_length + gap_energy*gap_energy ) return np.tanh( ( ( b_const / temp ) * squareRoot ) / squareRoot ) def integrate( coherence_length, temp ): # Integrates equation f with respect to E, between 0 and 1. return integrate.quad( f, 0, 1, args = ( temp, ) )[0] def root_of_integral( temp ): # Finds the roots of the integral with a guess of 0.01. return fsolve( integrate, 0.01, args = ( temp, ) ) def gap_energy_values( temp ): # Subtracts a_const from each root found, to obtain the gap_energy_values. return root_of_integral( temp ) - a_const </code></pre>
0 条评论
分类:
Python问答
请先
登录
后评论
默认排序
时间排序
1 个回答
匿名
1天前
擅长:python、mysql、java
<p>这条线:</p> <pre><code>integral = (integrate.quad(lambda E: np.tanh(1477.92*np.sqrt(E**2+x**2))/np.sqrt(E**2+x**2), 0, 1) </code></pre> <p>括号不平衡:</p> ^{pr2}$ <p>如果你把它拆散了会容易得多</p> ^{3}$
请先
登录
后评论
针对此问题:
更多的回答
关注
89
关注
收藏
1
收藏,
216
浏览
网友 提问于 2天前
相关Python问题
无法使用Django/mongoengine连接到MongoDB(身份验证失败)
2 回答
无法使用Django\u mssql\u后端迁移到外部hos
4 回答
无法使用Django&Python3.4连接到MySql
5 回答
无法使用Django+nginx上载媒体文件
1 回答
无法使用Django1.6导入名称模式
1 回答
无法使用Django1.7和mongodb登录管理站点
5 回答
无法使用Djangoadmin创建项目,进程使用了错误的路径,因为我事先安装了错误的Python
8 回答
无法使用Djangockedi验证CBV中的字段
8 回答
无法使用Djangocketditor上载图像(错误400)
6 回答
无法使用Djangocron进行函数调用
10 回答
无法使用Djangofiler djang上载文件
1 回答
无法使用Djangokronos
2 回答
无法使用Djangomssql provid
1 回答
无法使用Djangomssql连接到带有Django 1.11的MS SQL Server 2016
7 回答
无法使用Djangomssq迁移Django数据库
8 回答
无法使用Djangonox创建用户
4 回答
无法使用Djangopyodb从Django查询SQL Server
8 回答
无法使用Djangopython3ldap连接到ldap
9 回答
无法使用Djangoredis连接到redis
2 回答
无法使用Django中的FK创建新表
5 回答