function xf = interpolate(x)
  %Interpolación lineal de centroides a caras para elementos equidistantes. Elementos de frontera, interpolación de orden 0. 
  xf = [x;x(end)];
  xf(2:end-1) = 0.5*(x(1:end-1) + x(2:end));
endfunction
