One surface, the so-called Cayley Cubic (with four ordinary nodes as singularities): cayley_simple.sux

The Cayley Cubic together with its Hessian surface (i.e. the determinant of the matrix of its second partial derivatives): cayley_hesse.sux