ix1=1:nx; ix2=(nx+1):(2*nx); ix3=(2*nx+1):(3*nx); ix4=(3*nx+1):(4*nx);
iy1=1:nx; iy2=(nx+1):(2*nx); iy3=(2*nx+1):(3*nx); irv=nx:-1:1;