Quad Normals

How do I calculate Quad normals?

Using perpnorm from glAux vect3d.h