We propose a variational method for determining the surface patterns of cylindrical gels for both swelling and shrinking. Exact solutions are calculated for the initial stages of such peristaltic instabilities. The morphology and the formation mechanisms depend on a competition between bulk elastic energy and surface tension.