consensus of multi agent system control

%==============
%2016.4.2
%===============
function consensus
clc
clear
%%begin
N=input（‘请输入N的值:‘）;
a=cell（NN）;%创建一个空的N*N矩阵
X=cell（1N）;
Y=cell（1N）;
d=cell（1N）;
B=round（rand（NN））;%由rand函数产生随机数并四舍五入
A=tril（B-1）+triu（B‘0）;%产生N*N对称矩阵
L=diag（sum（A2））-A;%L=D-A
eig（L）;%求L的全部特征值

for i=1:N
    X{1i}=i;
    Y{1i}=0;
end

for I=1:N
    d{1I}=[];
end

T0=0;
n=[];

for k=1:1:1200
    for j=1:N
        b=0;
        for m=1:N
            b=b+A（jm）*（Y{1m}-Y{1j}）;%公式（2）累加
        end