diff --git a/fvnlib.f90 b/fvnlib.f90
index 460b9c8..512a22e 100644
--- a/fvnlib.f90
+++ b/fvnlib.f90
@@ -1626,7 +1626,7 @@ subroutine fvn_s_akima(n,x,y,br,co)
     ! br is just a copy of x
     br(:)=x(:)
     
-    allocate(var(n))
+    allocate(var(n+3))
     allocate(z(n))
     ! evaluate the variations
     do i=1, n-1
@@ -1688,7 +1688,7 @@ subroutine fvn_d_akima(n,x,y,br,co)
     ! br is just a copy of x
     br(:)=x(:)
 
-    allocate(var(n))
+    allocate(var(n+3))
     allocate(z(n))
     ! evaluate the variations
     do i=1, n-1