1、求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
print(filter(is_prime, range(1, 1000000)))
2、求任意整数的级数
Sympy支持级数和非级数求和(sum)运算,也可以直接用它的Sum类。Sum类的功能还是比较丰富的,支持很多级数运算,比如:“is_absolutely_convergent()”函数和“is_convergent()” 函数,可以直接判断一个级数是否收敛。
注:Sympy除了求和,类似的也提供了累乘计算类 —— Product。
from sympy import *
init_printing()
n = Symbol('n', integer=True)
expr = (2 * n - 1) / 2 ** n
Sum(expr, (n, 1, oo)), Sum(expr, (n, 1, oo)).doit()
(∑n=1∞2−n(2n−1),3)
本文转自xinrenbaodao 51CTO博客,原文链接:http://blog.51cto.com/11832904/2046858,如需转载请自行联系原作者