How do I plot P( U, V | C = c) in Python?

Question

My ultimate goal is to get the values that will then be used to visualize "the maps" of

• P(U, V | C = c)
• P(U | C = c) * P(V | C = c)
• |P(U, V | C = c) - P(U | C = c) * P(V | C = c)|

Here is some information about the data that I have.

Both U and V are continuous RVs. C I believe is ordered discreet as it consists only of the following values and the order of those values matters, as well as 5 > 2 == True.

C = {2, 3, 4, 5}

Otherwise, structurally all 3 variables have the same metadata

class:  ndarray
shape:  (90,)

I know how to plot

plt.hist(U, bins = 5, density = True)

and

plt.hist(V, bins = 5, density = True)

But how do I plot P(U, V | C = c)? Or P(U | C = c) * P(V | C = c)? Or |P(U, V | C = c) - P(U | C = c) * P(V | C = c)|? How do I get at those probabilities for each of the datapoints?