The aim of this modelling work is to assess shape changes of cometary nuclei caused by sublimation of ices. The simplest possible model is assumed with the nucleus being initially spherical and its thermal conductivity being neglected. We have calculated the time-dependent sublimation flux versus cometographic latitude. If the rotation axis of the comet is inclined to the orbital plane, then sublimation leads to non-symmetrical changes of the nucleus shape. Calculations were performed for the nuclei of comets Hale–Bopp and Wirtanen.