Adding an extra point in a grid in finite element code

Dears,

I am using this code that solves the 1D BVP -au'' + bu' + cu = f in (x0, x1) at x0 to x1.

I am diving the intervals to n=16.

I need to add an extra node (point) between the last two nodes so it will be between 15/16 and 16/16.

Any Idea how to do so?

Here is the part of the code concerning the intervals and the matrix forming:

n = 16;

x = linspace(leftbdr.x, rightbdr.x, n+1);

% assemble matrix A

A = sparse(n+1,n+1);

for i=1:n % on each sub-interval [x(i), x(i+1)]

localmat = assembleLocalMatrix(x(i),x(i+1),coefficients);

A(i:i+1,i:i+1) = A(i:i+1,i:i+1) + localmat;

end

% assemble right-hand side vector rhs = (f,v)

rhs = zeros(n+1,1);

for i=1:n % on each sub-interval [x(i), x(i+1)]

localrhs = assembleLocalRhs(x(i),x(i+1),f);

rhs(i:i+1) = rhs(i:i+1) + localrhs;

end

Thanks