function dp=crb(t,p,c,e)
dp=c*p*(1-p)-e*p
end
c=0.05;
e=0.025;
C=0:0.02:0.9;
nC=length(C);
figure
for i=1:nC
[t,p]=ode45(@crb,[0,300],C(i),[],c,e);
plot(t,p)
hold on
end
axis([0.1 300 0 1])
三维：让e=0.03，0.04,0.05,0.06,0.07,0.08,0.09，0.1,0.2,0.3,0.4,0.5,0.6变化在一张图上效果是：
