/* pure additive sequence */
/* 'seqa' in GAUSS has some error if its increment is a decimal. */
/* This comes from 'for' command. Try avoid 'for' loop in GAUSS. */
new; cls;
x=pureseqa(0,0.001,1001);
y=seqa(0,0.001,1001);
x[1001]==1;
y[1001]<=1;




proc pureseqa(start,inc,n);
   local x,e;
   e=1e-16;
   x=start|(start+inc*seqa(1,1,n-1));
   x=x.*(abs(x-round(x))>=e)+round(x).*(abs(x-round(x))<e);
   retp(start|(start+inc*seqa(1,1,n-1)));
endp;