Compute the normal of a vertex by averaging the normals of all triangles adjacent to that vertex.
One final normalization of the sum of all normals merely ensures that the result has unit length.
With one normal for each vertex, the model looks smooth.