p1=c(.3,.3,.1,.3)
p2=c(.3,.3,.1,.3)
x1=c(10,10,2,8)
x2=c(5,8,1,6)
N1=30
N2=20
dmultinom(x1,N1,p1)*dmultinom(x2,N2,p2)