function Out=usMatEc(dx,n,mc,indsEc)
%        Out=usMatEc(dx,n,mc,indsEc)
Out=spalloc(n,n,2*mc);
Out(indsEc)=dx/sqrt(2);
Out=Out+Out';
