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