我一直试图使用Sage Math绘制1
和一些非负数x
的算术几何平均值,但我很难让plot
函数接受一个Python函数作为输入:
plot(agm(1,x), (x,0,5))
产生以下错误:
TypeError: Cannot evaluate symbolic expression to a numeric value.
agm(x,y)
的定义如下:
def agm(x,y):
x,y = RDF(x),RDF(y)
a = (x+y)/2
g = sqrt(x*y)
diff = a-g
while diff > (1/10**11):
newa = (a + g)/2
newg = sqrt(a*g)
a,g = newa,newg
diff = a - g
return a
我在绘图时尝试用x
代替x.n()
,但这并没有改变任何内容。
如何绘制这个Python函数?