The first implementation of analytic gradients for the coupled-cluster singles, doubles, triples (CCSDT) model is described. The relevant theoretical expressions are given in a diagrammatic form together with the corresponding algebraic formulas. The computational requirements of CCSDT gradient calculations are discussed and their applicability demonstrated by performing benchmark calculations for molecular geometries with large correlation-consistent basis sets. A statistical analysis of the data reveals that CCSDT and CCSD(T) in most cases perform equally well. The CCSDT calculations thus provide further evidence for the high accuracy of the CCSD(T) approach.