function dx=ode1(t,x) %主函數(shù)
global m n a b c d1 d2;
% m=5;%“魚”的數(shù)量
% n=2;%維度-平面
% a=1;
% b=20;
% c=0.2;
% d1=0.1;
% d2=1;
dx=zeros((m+1)*n,1);%產(chǎn)生微分函數(shù)矩陣，整個微分函數(shù)放到一列。
%zeros是產(chǎn)生一個全部為零的矩陣，這里是12*1的矩陣，意思是首先全部清零
for i=1:m   %循環(huán)開始
    for j=1:n %第i個魚的第j維坐標(biāo)
       for k=1:m
            g=zeros(n,1);%產(chǎn)生
            g=x(1+(i-1)*n : n+(i-1)*n)- x(1+(k-1)*n : n+(k-1)*n);%產(chǎn)生向量的差，xi與xk的差。
            dx(j+(i-1)*n)=dx(j+(i-1)*n)+(x(j+(i-1)*n)-x(j+(k-1)*n))*(-a+b*exp(-(g'*g)/c));%產(chǎn)生"魚"的運(yùn)動方程中g(shù)
       end
       dx(j+(i-1)*n)=dx(j+(i-1)*n)+d1*(x(j+(i-1)*n)-x(j+m*n)); %產(chǎn)生"魚"的運(yùn)動方程f 
       dx(j+m*n)=dx(j+m*n)-d2*(x(j+m*n)-x(j+(i-1)*n)); %敵人的運(yùn)動方程         
    end
end