clc
clear

global m n a b c d1 d2 ;

m=5;%“魚”的數量
n=3;
a=1;
b=20;
c=0.2;
d1=0.1;
d2=1;
S=50;%活動范圍
x0=[S*(rand(m*n,1));0;-200;-200];%產生一個隨機的矩陣12*2，所以每次畫出來的圖是不一樣的
[T,X]=ode45(@ode1,[0 15],x0);%這就是求解這個微分方程，用的是ode45函數,dx定義用的是ode1函數,0-100是時間范圍
for i=1:3:3*(m)
    plot3(X(:,i),X(:,i+1),X(:,i+2));
    hold on
end;
    plot3(X(:,3*m+1),X(:,3*m+2),X(:,3*m+3),'g');

for i=1:3:3*(m+1)
   comet3(X(:,i),X(:,i+1),X(:,i+2),0.01);
hold on
end;