diff --git a/lib/polyfit.f90 b/lib/polyfit.f90
index 779b99a9b..632223048 100644
--- a/lib/polyfit.f90
+++ b/lib/polyfit.f90
@@ -1,7 +1,7 @@
 subroutine polyfit(x,y,sigmay,npts,nterms,mode,a,chisqr)
   implicit real*8 (a-h,o-z)
   real*8 x(npts), y(npts), sigmay(npts), a(nterms)
-  real*8 sumx(19), sumy(10), array(10,10)
+  real*8 sumx(10), sumy(10), array(10,10)
 
 ! Accumulate weighted sums
   nmax = 2*nterms-1