1 2 3 4 5 6 7 8 9 10 | from ptrap_gen import trap0,trap1
import math
def f_r_12_x(r_12_x):
v = trap1(1.00000000000000e-5,1.00000000000000,f_r_12_y,n,r_12_x)
return v
def f_r_12_y(r_12_x,r_12_y):
return math.exp(-4.0 * ((r_12_y**2 + r_12_x**2)**-6) + 4.0 * ((r_12_y**2 + r_12_x**2)**-3))
n = 10
v = trap0(1.00000000000000e-5,1.00000000000000,f_r_12_x,n)
print v
|