In this paper we present a novel scheme for surface generation. Every surface patch satisfy a 8th-order partial differential equation(PDE)with a shape parameter $c$ and eight boundary curves, which make the combined surface be a $C^1$ continuous surface. Moreover, when the shape parameter $c$ is given a low value, the proposed PDE yields a class of surfaces including bicubic surface; when $c$ is given a large value, the surface is approximated to the one generated by the Euler-Lagrange PDE $S_{xxyy}=0$, which include bilinear surfaces. Finally, our computational examples demonstrate the validity of the algorithm.