Python 离散和#x2013的Sage实现;时间动力系统
因此,正如您所看到的,我有Python 离散和#x2013的Sage实现;时间动力系统,python,algorithm,sage,Python,Algorithm,Sage,因此,正如您所看到的,我有x1和x2的初始值,然后将其插入函数f1和f2,然后我调用了一个新变量k1,它就是f1+x1,我对f2做了同样的事情。然后,我得到的k1和k2成为x1和x2的新值,依此类推。我想这样做好几次,但我觉得没有for循环为N迭代次数这样做有点愚蠢。我正试图在圣人身上做到这一点。我在Sage中查找了一些编码,比如使用牛顿方法。但我看到他们在编程中都使用一个变量。有人能让我开始编写这个问题的代码吗。原则上,您只需重新命名即可。在这里,我将-*x1替换为-x1,尽管根据我得到的数字
x1
和x2
的初始值,然后将其插入函数f1
和f2
,然后我调用了一个新变量k1
,它就是f1+x1
,我对f2
做了同样的事情。然后,我得到的k1
和k2
成为x1
和x2
的新值,依此类推。我想这样做好几次,但我觉得没有for循环为N
迭代次数这样做有点愚蠢。我正试图在圣人身上做到这一点。我在Sage中查找了一些编码,比如使用牛顿方法。但我看到他们在编程中都使用一个变量。有人能让我开始编写这个问题的代码吗。原则上,您只需重新命名即可。在这里,我将-*x1
替换为-x1
,尽管根据我得到的数字,我认为这是不对的
x1=18;x2 = 10
f1_0=(40)*x1*(100000-x1-x2)-x1;k1_0=f1_0+x1;k1_0.n()
f2_0=(1/3)*x1;k2_0=f2_0+x2;k2_0.n()
######################################################
x1=k1_0.n();x2=k2_0.n()
f1_1=(40)*x1*(100000-x1-x2)-x1;k1_1=f1_1+x1;k1_1.n()
f2_1=(1/3)*x1;k2_1=f2_1+x2;k2_1.n()
######################################################
x1=k1_1.n();x2=k2_1.n()
f1_2=(40)*x1*(100000-x1-x2)-x1;k1_2=f1_2+x1;k1_2.n()
f2_2=(1/3)*x1;k2_2=f2_2+x2;k2_2.n()
什么是
-*x1
?你不能同时进行减法和乘法。好消息!如果这回答了您的问题,那么一定要单击复选标记,以便其他人知道这解决了它-并为自己获取Python教程,我认为这将真正为您带来回报!我做到了。你知道我如何将所有的x1都放在一个列表中,这样我就能够绘制它们和x2的相同内容。关于Stackoverflow,有很多问题可以帮助你做到这一点-看看如何创建Python列表并向其添加内容!正如我所说,一个好的Python教程将对您非常有帮助。提示:points(zip(list1,list2))
一旦您列出了x1和x2值。
x1=18;x2 = 10
for i in range(10):
f1=(40)*x1*(100000-x1-x2)-x1
k1=f1+x1
f2=(1/3)*x1
k2=f2+x2
print k1, k2
x1=k1.n()
x2=k2.n()