一、求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