Polyolefins made with Ziegler-Natta catalysts have non-uniform distributions of molecular weight (MWD) and chemical composition (CCD). The MWD is usually measured by high-temperature gel permeation chromatography (GPC) and the CCD by either temperature rising elution fractionation (TREF) or crystallization analysis fractionation (CRYSTAF). A mathematical model is needed to quantify the information provided by these analytical techniques and to relate it to the presence of multiple site types on Ziegler-Natta catalysts. We developed a robust computer algorithm to deconvolute the MWD and CCD of polyolefins simultaneously using Flory's most probable distribution and the cumulative CCD component of Stockmayer's distribution, which includes the soluble fraction commonly present in linear low-density polyethylene (LLDPE) resins and have applied this procedure for the first time to several industrial LLDPE resins. The deconvolution results are reproducible and consistent with theoretical expectations.