function [i,j,x,n] = sparse_vec(i,j,x,n, ii,jj,xx)

    n = n + 1;
    L = length(i);
    if (n > L)
       i(2*L) = 0;
       j(2*L) = 0;
       x(2*L) = 0;
    end
    
    i(n) = ii;
    j(n) = jj;
    x(n) = xx;