Modeling of the bivariate molecular weight distribution-copolymer composition distribution in RAFT copolymerization using probability generating functions