for(k=i-m;k<i;k++)

                for(l=j-k+i;l<j;l++)

                {d11[1]+=(AA[k*width+l]-u0[1])*(AA[k*width+l]-u0[1]);}

                d11[1]=d11[1]/(2*m+1)/(2*n+1)*8;

            //////////////////////////////

            for(k=i-m;k<i;k++)

                for(l=j-k+i;l<j+m;l++)

                {u0[2]+=AA[k*width+l];}

                u0[2]=u0[2]/(2*m+1)/(2*n+1)*8;

            for(k=i-m;k<i;k++)

                for(l=j-k+i;l<j+m;l++)

                {d11[2]+=(AA[k*width+l]-u0[2])*(AA[k*width+l]-u0[2]);}

                d11[2]=d11[2]/(2*m+1)/(2*n+1)*8;

            //////////////////////////////

            for(k=i;k<i+m;k++)

                for(l=j+k-i;l<j+m;l++)

                {u0[3]+=AA[k*width+l];}

                u0[3]=u0[3]/(2*m+1)/(2*n+1)*8;

            for(k=i;k<i+m;k++)

                for(l=j+k-i;l<j+m;l++)

                {d11[3]+=(AA[k*width+l]-u0[3])*(AA[k*width+l]-u0[3]);}

                d11[3]=d11[3]/(2*m+1)/(2*n+1)*8;

            //////////////////////////////

            for(k=i;k<i+m;k++)

                for(l=j;l<j+k-i;l++)

                {u0[4]+=AA[k*width+l];}

                u0[4]=u0[4]/(2*m+1)/(2*n+1)*8;

            for(k=i;k<i+m;k++)

                for(l=j;l<j+k-i;l++)

                {d11[4]+=(AA[k*width+l]-u0[4])*(AA[k*width+l]-u0[4]);}

                d11[4]=d11[4]/(2*m+1)/(2*n+1)*8;

            //////////////////////////////

            for(k=i;k<i+m;k++)

                for(l=j-k+i;l<j+m;l++)

                {u0[5]+=AA[k*width+l];}

                u0[5]=u0[5]/(2*m+1)/(2*n+1)*8;

            for(k=i;k<i+m;k++)

                for(l=j-k+i;l<j+m;l++)

                {d11[5]+=(AA[k*width+l]-u0[5])*(AA[k*width+l]-u0[5]);}

                d11[5]=d11[5]/(2*m+1)/(2*n+1)*8;

            //////////////////////////////