Geometric piecewise cubic Bézier interpolating polynomial with $C^2$ continuity

Bézier curve is a parametric polynomial that is applied to produce good piecewise interpolation methods with more advantage over the other piecewise polynomials. It is, therefore, crucial to construct Bézier curves that are smooth and able to increase the accuracy of the solutions. Most of the known strategies for determining internal control points for piecewise Bezier curves achieve only partial smoothness, satisfying the first order of continuity. Some solutions allow you to construct interpolation polynomials with smoothness in width along the approximating curve. However, they are still unable to handle the locations of the inner control points. The partial smoothness and non-controlling locations of inner control points may affect the accuracy of the approximate curve of the dataset. In order to improve the smoothness and accuracy of the previous strategies, а new piecewise cubic Bézier polynomial with second-order of continuity $C^2$ is proposed in this study to estimate missing values. The proposed method employs geometric construction to find the inner control points for each adjacent subinterval of the given dataset. Not only the proposed method preserves stability and smoothness, the error analysis of numerical results also indicates that the resultant interpolating polynomial is more accurate than the ones produced by the existing methods.