2

我正在使用以下代码：

eps = MkInfinitesimal()
print(( (1 + eps)**2- 1**2)/eps < 2.00000000001)
print(( (1 + eps)**2- 1**2)/eps > 2)

输出是：

True
True

其他示例：证明当 x = pi 时 diff(x^2, x) = 2*pi 和当 x = e 时 diff(x^2,x) = 2*e

代码：

eps1 = MkInfinitesimal()
eps2 = MkInfinitesimal() # eps2 is infinitely smaller thant eps1
pi = Pi()
e = E()
print(( (pi + eps2)**2- pi**2)/eps2 < 2*pi + eps1)
print(( (pi + eps2)**2- pi**2)/eps2 > 2*pi)
print(( (e + eps2)**2- e**2)/eps2 < 2*e + eps1)
print(( (e + eps2)**2- e**2)/eps2 > 2*e)

输出：

True
True
True
True

其他示例：证明当 x = e 或 x = pi 时 diff(x^3, x) = 3*x^2。

代码：

eps1 = MkInfinitesimal()
eps2 = MkInfinitesimal() # eps2 is infinitely smaller thant eps1
pi = Pi()
e = E()
print(( (pi + eps2)**3- pi**3)/eps2 < 3*pi**2 + eps1)
print(( (pi + eps2)**3- pi**3)/eps2 > 3*pi**2)
print(( (e + eps2)**3- e**3)/eps2 < 3*e**2 + eps1)
print(( (e + eps2)**3- e**3)/eps2 > 3*e**2)

输出：

True
True
True
True

其他示例：

代码：

[x] = MkRoots([-1, -1, 0, 0, 0, 1])
[y] = MkRoots([-197, 3131, -31*x**2, 0, 0, 0, 0, x])
[z] = MkRoots([-735*x*y, 7*y**2, -1231*x**3, 0, 0, y])
print(x.decimal(10))
print(y.decimal(10))
print(z.decimal(10))
eps1 = MkInfinitesimal()
eps2 = MkInfinitesimal() # eps2 is infinitely smaller thant eps1
print(( (x + eps2)**2- x**2)/eps2 < 2*x + eps1)
print(( (x + eps2)**2- x**2)/eps2 > 2*x)
print(( (y + eps2)**2- y**2)/eps2 < 2*y + eps1)
print(( (y + eps2)**2- y**2)/eps2 > 2*y)
print(( (z + eps2)**2- z**2)/eps2 < 2*z + eps1)
print(( (z + eps2)**2- z**2)/eps2 > 2*z)

输出：

1.1673039782?
0.0629726948?
31.4453571397?
True
True
True
True
True
True

这个证明对吗？如果您知道更好的证明，请告诉我。非常感谢。

4

1 回答 1

2

这是一个很好的例子。Z3RCF API 将在 Z3Py 的下一个正式版本中提供。

顺便说一句，您可以在 Z3RCF 中创建许多无穷小。每一个都比之前创建的要小得多。这是相同的示例，但它2.0000000001通过使用两个不同的无穷小来避免 （也可以在此处获得）。

eps1 = MkInfinitesimal()
eps2 = MkInfinitesimal() # eps2 is infinitely smaller thant eps1
print(( (1 + eps2)**2- 1**2)/eps2 < 2 + eps1)
print(( (1 + eps2)**2- 1**2)/eps2 > 2)
于 2013-05-16T17:32:12.193 回答