|
极限问题
这是我的数值计算所用的代码。 python 是现在很流行的一种程序语言。
可以免费下载使用。代码无需编译就可执行。特此推荐
#!F:/Python/python
# Some study of |sin n|
import math
def asn(n):
return abs(math.sin(n+0.0))
def acs(n):
return abs(math.cos(n+0.0))
def pslt(n,f=';s';):
row = n/10 + 1 # number of lines
rcnt = 0 # number of current line items
rr = 0 # curent line number, start with 0
ps = 0 # position at current line, start with 0
while rr < row:
rcnt = max(10*min((n-10*rr)/10,1),n%10)
ps = 0
for ps in range(rcnt):
if f == ';s';:
print("%8.5f" % asn(10*rr+ps+1)),
elif f == ';c';:
print("%8.5f" % acn(10*rr+ps+1)),
print("")
rr += 1
def sta(n,f=';s';,sb = 20):
c = []
a = []
b = []
sbr = sb + 0.0 # make sure it';s float, not integer
for i in range(sb):
c.append(0)
a.append(i/sbr)
b.append((i+1)/sbr)
for k in range(n):
if f == ';s';:
x = asn(k)
elif f == ';c';:
x = acs(k)
for i in range(sb):
if a <= x and x < b:
c = c+1
break
print(c)
|
|