function daoX=daoshu(X,V,d,N)
global A
global B
global C
global D
  for x=1:N 
     for i=1:N 
          m1=0;  m2=0;  m3=0;  m4=0; 
                for j=1:N 
                    if j~=i 
                        m1=m1+V(x,j); 
                    end                    
                end 
                  
                for y=1:N 
                    if y~=x 
                        m2=m2+V(y,i); 
                    end 
                end 
               
                for y=1:N 
                    for j=1:N 
                        m3=m3+V(y,j); 
                    end 
                end 
        
                for y=1:N 
                    if y~=x 
                        if i==1 
                            m4=m4+d(x,y)*(V(y,i+1)+V(y,N)); 
                        elseif i==N 
                            m4=m4+d(x,y)*(V(y,i-1)+V(y,1)); 
                        else 
                            m4=m4+d(x,y)*(V(y,i+1)+V(y,i-1)); 
                        end 
                    end   
                end 
               
                daoX(x,i)=-X(x,i)-A*m1-B*m2-C*(m3-N)-D*m4; 
            end 
     end