A local Lagrange interpolation method based on $C^1$ cubic