  We describe a novel method for obtaining the gradient expansion for the free
energy of a clean BCS superconductor. We present explicit results up to fourth
order in the gradients of the order parameter.
