一、求100萬之內的全部素數
def is_prime(n):
if n == 1:
return False
for i in range(2, n//2 + 1):
if n % i == 0:
return False
return True ide
print(filter(is_prime, range(1, 1000000)))函數
二、求任意整數的級數
Sympy支持級數和非級數求和(sum)運算,也能夠直接用它的Sum類。Sum類的功能仍是比較豐富的,支持不少級數運算,好比:「is_absolutely_convergent()」函數和「is_convergent()」 函數,能夠直接判斷一個級數是否收斂。
注:Sympy除了求和,相似的也提供了累乘計算類 —— Product。
from sympy import *
init_printing()it
n = Symbol('n', integer=True)
expr = (2 * n - 1) / 2 ** n
Sum(expr, (n, 1, oo)), Sum(expr, (n, 1, oo)).doit()class
(∑n=1∞2−n(2n−1),3)import