--- a/mp_BBP.m
+++ b/mp_BBP.m
@@ -22,8 +22,7 @@
   i8 = 8 * mp (ii, precision);
   x = 4/(i8 + 1)-2/(i8 + 4)-1/(i8 + 5)-1/(i8 + 6);
   x = x * power (mp (16, precision), -ii);
-  terms = zeros([1 1+mp(order, precision)]);
-  [out, terms(end)] = deal (x);
+  [out, terms(1, 1+order)] = deal (x);
 
   %# compute the rest
   for ii = (order-1:-1:0)