Let L be an elliptic system of higher order homogeneous partial differ
ential operators. We establish in this article the equivalence in L-P
norm between the maximal function and the square function of solutions
to L in Lipschitz domains. Several applications of this result are di
scussed.