How can I simulate these joint random variables with the given PMF?

Question

If I have two Discrete Random variables X and Y which have the joint PMF:

Pxy=c(x2+y2), where x belongs to {1,2,4} and y belongs to {1,3}

I tried to make two d