function rho = L2_x(rho)
global delx;
rho = sum(rho.^2)*delx;