function [a,b,c]=reeval(delta,a,b,c,N)
for i=1:N
   a(i)=a(i)+delta(1,i);
   b(i)=b(i)+delta(2,i);
   c(i)=c(i)+delta(3,i);
end