|
|

楼主 |
发表于 2019-6-11 17:53:15
|
显示全部楼层
import math
def Q(r):
return 1e14*math.cos(math.pi*r/(2*R))
D = 0.008
E = 8.5
R = 1
alpha = int(input("请输入(球2柱1板0):"))
N = 100
h = 0.01
AA = []
BB = []
CC = []
DD = []
for i in range(1,N-1):
AA.append(((i-1)*h+h/2)**alpha*D/h)
CC.append((i*h+h/2)**alpha*D/h)
BB.append((i*h)**alpha*E*h+ AA[i-1] + CC[i-1])
DD.append((i*h)**alpha/2*h*(Q((i-1)*0.01+h/2)+Q(i*0.01+h/2)))
U = []#贝塔
P = []#得它
U.append(BB[0])
P.append(DD[0])
for i in range(1,N-1):
U.append(BB[i] - CC[i-1]**2/U[i-1])
P.append(DD[i] - CC[i-1]*P[i-1]/U[i-1])
|
|