This article presents a method for constructing a G 1 -smooth surface, composed of independently parametrized triangular polynomial Bézier patches, to fit scattered data points triangulated in R 3 with arbitrary topology. The method includes a variational technique to optimize the shape of the surface. A systematic development of the method is given, presenting general equations provided by the theory of manifolds, explaining the heuristic assumptions made to simplify calculations, and analyzing the numerical results obtained from fitting two test configurations of scattered data points. The goal of this work is to explore an alternative G 1 construction, inspired by the theory of manifolds, that is subject to fewer application constraints than approaches found in the technical literature; e.g., this approach imposes no artificial restrictions on the tangents of patch boundary curves at vertex points of a G 1 surface. The constructed surface shapes fit all test data surprisingly well for a noniterative method based on polynomial patches.