We outline polarization performance calculations and predictions for the Daniel K. Inouye Solar Telescope (DKIST) optics and show Mueller matrices for two of the first light instruments. Telescope polarization is due to polarization dependent mirror reflectivity and rotations between groups of mirrors as the telescope moves in altitude and azimuth. The Zemax optical modeling software has polarization ray-trace capabilities and predicts system performance given a coating prescription. We develop a model coating formula that approximates measured witness sample polarization properties. Estimates show the DKIST telescope Mueller matrix as functions of wavelength, azimuth, elevation, and field angle for the Cryogenic Near Infra-Red Spectro-Polarimeter (CryoNIRSP) and for the Visible SpectroPolarimeter (ViSP). Footprint variation is substantial. We estimate 2% variation of some Mueller matrix elements over the 5 arc minute CryoNIRSP field. We validate the Zemax model by show limiting cases for flat mirrors in collimated and powered designs that compare well with theoretical approximations and are testable with lab ellipsometers.