I have already computed a 6-dimensional joint distribution function in Python using the pyvinecopulib library and numpy, and stored it in the variable cop. How can I calculate the probability on an interval such as (0.5 < F1(X1) < 0.6, 0.5 < F2(X2) < 0.6, 0.5 < F3(X3) < 0.6, 0.5 < F4(X4) < 0.6, 0.5 < F5(X5) < 0.6, 0.6 < F1(X6) < 0.6)?