function ut=utilda(u,i,j);
ut=ones(u)*nan;
ut(i,j)=(u(i-1,j)+u(i,j)+u(i+1,j))/3;
