diff --git a/Principal_component_analysis/include/CGAL/linear_least_squares_fitting_segments_2.h b/Principal_component_analysis/include/CGAL/linear_least_squares_fitting_segments_2.h
index d7cf0448b93..8c0fbb8c7b3 100644
--- a/Principal_component_analysis/include/CGAL/linear_least_squares_fitting_segments_2.h
+++ b/Principal_component_analysis/include/CGAL/linear_least_squares_fitting_segments_2.h
@@ -104,7 +104,7 @@ linear_least_squares_fitting_2(InputIterator first,
 
   // Translate the 2nd order moment calculated about the origin to
   // the center of mass to get the covariance.
-  covariance[0] -= mass * CGAL::square(c.x()));
+  covariance[0] -= mass * CGAL::square(c.x());
   covariance[1] -= mass * (c.x() * c.y());
   covariance[2] -= mass * (CGAL::square(c.y()));