for p in range(10):
for q in range(10):
print("p=%s, q=%s, f=%s", p/10+0.1, q/10+0.1, f(p/10+0.1,q/10+0.1))
p=%s, q=%s, f=%s 0.1 0.1 0.13123614957214594
p=%s, q=%s, f=%s 0.1 0.2 0.17222926954522666
p=%s, q=%s, f=%s 0.1 0.30000000000000004 0.22602706178738802
p=%s, q=%s, f=%s 0.1 0.4 0.2966292128808233
p=%s, q=%s, f=%s 0.1 0.5 0.38928475749095637
p=%s, q=%s, f=%s 0.1 0.6 0.510882326602397
p=%s, q=%s, f=%s 0.1 0.7 0.6704622942775809
p=%s, q=%s, f=%s 0.1 0.7999999999999999 0.8798888993429671
p=%s, q=%s, f=%s 0.1 0.9 1.1547323120104451
p=%s, q=%s, f=%s 0.1 1.0 1.5154262241479266
p=%s, q=%s, f=%s 0.2 0.1 0.26247229914429193
p=%s, q=%s, f=%s 0.2 0.2 0.3444585390904532
p=%s, q=%s, f=%s 0.2 0.30000000000000004 0.45205412357477603
p=%s, q=%s, f=%s 0.2 0.4 0.5932584257616464
p=%s, q=%s, f=%s 0.2 0.5 0.7785695149819126
p=%s, q=%s, f=%s 0.2 0.6 1.0217646532047937
p=%s, q=%s, f=%s 0.2 0.7 1.3409245885551615
p=%s, q=%s, f=%s 0.2 0.7999999999999999 1.759777798685934
p=%s, q=%s, f=%s 0.2 0.9 2.3094646240208903
p=%s, q=%s, f=%s 0.2 1.0 3.0308524482958528
p=%s, q=%s, f=%s 0.30000000000000004 0.1 0.3937084487164379
p=%s, q=%s, f=%s 0.30000000000000004 0.2 0.5166878086356799
p=%s, q=%s, f=%s 0.30000000000000004 0.30000000000000004 0.678081185362164
p=%s, q=%s, f=%s 0.30000000000000004 0.4 0.8898876386424697
p=%s, q=%s, f=%s 0.30000000000000004 0.5 1.1678542724728689
p=%s, q=%s, f=%s 0.30000000000000004 0.6 1.5326469798071907
p=%s, q=%s, f=%s 0.30000000000000004 0.7 2.0113868828327424
p=%s, q=%s, f=%s 0.30000000000000004 0.7999999999999999 2.639666698028901
p=%s, q=%s, f=%s 0.30000000000000004 0.9 3.4641969360313354
p=%s, q=%s, f=%s 0.30000000000000004 1.0 4.546278672443779
p=%s, q=%s, f=%s 0.4 0.1 0.5249445982885838
p=%s, q=%s, f=%s 0.4 0.2 0.6889170781809064
p=%s, q=%s, f=%s 0.4 0.30000000000000004 0.9041082471495521
p=%s, q=%s, f=%s 0.4 0.4 1.186516851523293
p=%s, q=%s, f=%s 0.4 0.5 1.557139029963825
p=%s, q=%s, f=%s 0.4 0.6 2.0435293064095874
p=%s, q=%s, f=%s 0.4 0.7 2.681849177110323
p=%s, q=%s, f=%s 0.4 0.7999999999999999 3.5195555973718675
p=%s, q=%s, f=%s 0.4 0.9 4.61892924804178
p=%s, q=%s, f=%s 0.4 1.0 6.061704896591705
p=%s, q=%s, f=%s 0.5 0.1 0.6561807478607297
p=%s, q=%s, f=%s 0.5 0.2 0.8611463477261331
p=%s, q=%s, f=%s 0.5 0.30000000000000004 1.13013530893694
p=%s, q=%s, f=%s 0.5 0.4 1.4831460644041161
p=%s, q=%s, f=%s 0.5 0.5 1.9464237874547814
p=%s, q=%s, f=%s 0.5 0.6 2.554411633011984
p=%s, q=%s, f=%s 0.5 0.7 3.352311471387903
p=%s, q=%s, f=%s 0.5 0.7999999999999999 4.399444496714834
p=%s, q=%s, f=%s 0.5 0.9 5.773661560052224
p=%s, q=%s, f=%s 0.5 1.0 7.57713112073963
p=%s, q=%s, f=%s 0.6 0.1 0.7874168974328756
p=%s, q=%s, f=%s 0.6 0.2 1.0333756172713595
p=%s, q=%s, f=%s 0.6 0.30000000000000004 1.3561623707243278
p=%s, q=%s, f=%s 0.6 0.4 1.7797752772849391
p=%s, q=%s, f=%s 0.6 0.5 2.3357085449457373
p=%s, q=%s, f=%s 0.6 0.6 3.065293959614381
p=%s, q=%s, f=%s 0.6 0.7 4.022773765665484
p=%s, q=%s, f=%s 0.6 0.7999999999999999 5.279333396057801
p=%s, q=%s, f=%s 0.6 0.9 6.92839387206267
p=%s, q=%s, f=%s 0.6 1.0 9.092557344887554
p=%s, q=%s, f=%s 0.7 0.1 0.9186530470050215
p=%s, q=%s, f=%s 0.7 0.2 1.205604886816586
p=%s, q=%s, f=%s 0.7 0.30000000000000004 1.5821894325117158
p=%s, q=%s, f=%s 0.7 0.4 2.0764044901657623
p=%s, q=%s, f=%s 0.7 0.5 2.7249933024366935
p=%s, q=%s, f=%s 0.7 0.6 3.5761762862167776
p=%s, q=%s, f=%s 0.7 0.7 4.6932360599430645
p=%s, q=%s, f=%s 0.7 0.7999999999999999 6.159222295400768
p=%s, q=%s, f=%s 0.7 0.9 8.083126184073114
p=%s, q=%s, f=%s 0.7 1.0 10.607983569035483
p=%s, q=%s, f=%s 0.7999999999999999 0.1 1.0498891965771675
p=%s, q=%s, f=%s 0.7999999999999999 0.2 1.3778341563618128
p=%s, q=%s, f=%s 0.7999999999999999 0.30000000000000004 1.8082164942991037
p=%s, q=%s, f=%s 0.7999999999999999 0.4 2.3730337030465853
p=%s, q=%s, f=%s 0.7999999999999999 0.5 3.1142780599276496
p=%s, q=%s, f=%s 0.7999999999999999 0.6 4.087058612819174
p=%s, q=%s, f=%s 0.7999999999999999 0.7 5.363698354220644
p=%s, q=%s, f=%s 0.7999999999999999 0.7999999999999999 7.039111194743734
p=%s, q=%s, f=%s 0.7999999999999999 0.9 9.23785849608356
p=%s, q=%s, f=%s 0.7999999999999999 1.0 12.123409793183411
p=%s, q=%s, f=%s 0.9 0.1 1.1811253461493134
p=%s, q=%s, f=%s 0.9 0.2 1.5500634259070394
p=%s, q=%s, f=%s 0.9 0.30000000000000004 2.034243556086492
p=%s, q=%s, f=%s 0.9 0.4 2.669662915927409
p=%s, q=%s, f=%s 0.9 0.5 3.503562817418606
p=%s, q=%s, f=%s 0.9 0.6 4.597940939421571
p=%s, q=%s, f=%s 0.9 0.7 6.034160648498226
p=%s, q=%s, f=%s 0.9 0.7999999999999999 7.919000094086702
p=%s, q=%s, f=%s 0.9 0.9 10.392590808094004
p=%s, q=%s, f=%s 0.9 1.0 13.638836017331336
p=%s, q=%s, f=%s 1.0 0.1 1.3123614957214593
p=%s, q=%s, f=%s 1.0 0.2 1.722292695452266
p=%s, q=%s, f=%s 1.0 0.30000000000000004 2.26027061787388
p=%s, q=%s, f=%s 1.0 0.4 2.966292128808232
p=%s, q=%s, f=%s 1.0 0.5 3.8928475749095623
p=%s, q=%s, f=%s 1.0 0.6 5.108823266023968
p=%s, q=%s, f=%s 1.0 0.7 6.704622942775806
p=%s, q=%s, f=%s 1.0 0.7999999999999999 8.798888993429669
p=%s, q=%s, f=%s 1.0 0.9 11.54732312010445
p=%s, q=%s, f=%s 1.0 1.0 15.154262241479262
strider05052016-05-26 20:08:05 +08:00#!/usr/bin/python3from math import log, exp## ln(f)=ln(p)+q1*f=ln(p)+q1*exp[ln(f)]def f(p,q):lp=log(p)epsilon=1e-5## initial guessy0=0y1=1while (y1-y0)>epsilon ...
在实际的数学建模应用中,我们会遇到很多约束条件是二次的,三次的或者是高次函数的情况,这样用 optimize.linprog()来
解
决就显得不适用了,因此我们使用scipy.optimize下得minimize函数来
解
决这个问题。
官方文档:SciPy API referencehttps://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.minimize.html#scipy.optimize.minimize
scipy.o
def func(y):
return eval("""(23.474/ (1+y)**0.024 + 22.316/(1+y)**1.024 + 21.158/(1+y)**2.024-64.896)""")
root = optimize.newton(func, x0=0, tol=0.0000000001)
print(root)
转自https://
python
-china.org/t/1152def solve(eq, var):
eq1 = eq.replace("=", "-(") + ")"
c = eval(eq1, {var: 1j})
return -c.real / c.imag
原创内容,转载注明出处!
(上篇用
python
拟合2019nCov感染人数的文章被不少博主转载了,发的比较早,不少博主在文章基础上添加新内容也新发了新的更新后的预测,或者加入一些新的模块。博文链接如下:)
python
实现logistic增长模型拟合2019-nCov确诊人数
python
实现logistic增长模型拟合20...
Python
一元
二次方程
求
根
1、任务简介
在之前的博客中我分享了使用Java进行
一元
二次方程
求
根的方法,在学习了
Python
之后我也想使用
Python
编写一个类似的程序,故在编写成功后将该任务分享出来。
2、任务代码
学习过Java和
Python
的人都知道,
Python
的语法比Java简洁得多,并且目前已经广泛应用于爬虫开发、web开发、人工智能和机器学习等主流方向,是一种面向对象的语言,我学...
4*x0*x0 - 2*sin(x1*x2) = 0
x1*x2 - 1.5 = 0from scipy.optimize import fsolve
from math import sin,cosdef f(x):
x0 = float(x[0])
x1 = float(x[1])
要在
Python
中
求解
二阶微分方程,可以使用数值
解法
或
解
析
解法
。使用数值
解法
时,可以使用差分法,而使用
解
析
解法
时,可以根据方程的形式进行分类
求解
。
对于数值
解法
,可以使用差分法来
求解
一般的二阶线性常微分方程的边值问题。差分法将微分方程转化为差分方程,通过逼近微分方程的导数来
求解
。具体的数值
解
算例可以参考引用中提到的文章。
对于
解
析
解法
,可以根据二阶微分方程的形式进行分类
求解
。常见的形式包括可分离变量的微分方程、一阶齐次(非齐次)线性微分方程、二阶常系数微分方程和高阶常系数微分方程。针对不同的形式,可以采用不同的方法来
求解
。具体的
解
析
解
方法可以参考引用中提到的文章。
此外,基于
Python
的微分方程数值
解
也是可行的。可以使用
Python
中的数值计算库来实现对微分方程的数值
求解
。引用中的文章提供了常微分
方程解
析
解法
和基于
Python
的微分方程数值
解
算例的实现,可以作为参考和指导。
综上所述,要在
Python
中
解
二阶微分方程,可以选择使用数值
解法
(如差分法)或
解
析
解法
(根据不同的形式分类
求解
),也可以结合使用
解
析
解法
和数值
解法
。具体的方法和实现可以参考引用的资料。<span class="em">1</span><span class="em">2</span><span class="em">3</span>
#### 引用[.reference_title]
- *1* *2* *3* [常微分方程的
解
析
解
(方法归纳)以及基于
Python
的二阶微分方程边值问题的数值算例实现](https://blog.csdn.net/qq_42818403/article/details/120613079)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 100%"]
[ .reference_list ]