3、lve('Dy=y+2*x','y(0)=1','x');ezplot(y,[01])Result:1.2Code:clear[ty]=ode23('wf',[010],[01]);plot(t,y(:,1))M文件functionyd=wf(t,y)yd=[y(2);-y(1)*cos(t)]Result:以下为部分数值解和图像3．Rossler微分方程组：当固定参数b=2,c=4时，试讨论随参数a由小到大变化（如a∈(0,0.65))而方程解的变化情况，并且画出空间曲线图形，观察空间曲线是否形成混沌状？Code：clearglo
4、balafora=0:0.05:0.65;[tx]=ode45('P3',[010],[111]);plot(t,x)holdonendM文件：functionequ=fun(t,x)globalab=2;c=4;equ=[0-1-1;1a0;x(3)0-c]*x+[0;0;b];endResult：可见，图像呈混沌状。4.Apollo卫星的运动轨迹的绘制Code：clearglobaluu1u=1/82.45;u1=1-u;[tx]=ode45('P4',[010],[1.200-1.04935751]);plot(t,x)M文件
5、：functioneqt=P4(t,x)globaluu1r1=sqrt(((x(1)+u))^2+x(3)^2);r2=sqrt(((x(1)-u1)^2)+x(3)^2);eqt(1,1)=x(2);eqt(2,1)=2*x(4)+x(1)-u1*(x(1)+u)/r1^3-u*(x(1)-u1)/r2^3;eqt(3,1)=x(4);eqt(4,1)=-2*x(2)+x(3)-u1*x(3)/r1^3-u*x(3)/r2^3;Result：5．盐水的混合问题一个圆柱形的容器，内装350升的均匀混合的盐水溶液。如果纯水以每秒14
6、升的速度从容器顶部流入，同时，容器内的混合的盐水以每秒10.5升的速度从容器底部流出。开始时，容器内盐的含量为7千克。求经过时间t后容器内盐的含量。Code：clear[tx]=ode23('P5',0:0.1:10,[7350]);count=1;fort=0:0.1:10plot(t,x(count,1)/x(count,2),'--');xlabel('ʱ¼ä');ylabel('º¬ÑÎÁ¿£¨kg/L£©');holdoncount=count+1;endM文件：functionequ=P5(t,x)equ=[-10.5
7、*x(1)/x(2);3.5];result：9．River-bay系统水污染问题一条河流和河湾与大湖相连，位于湾上游的小河是造成湾污染的主要因素，另有一座铝厂恰好建在湾旁,也造成污染。当湾中污染物平均浓度达到1.6mg/l时，铝厂将被迫暂时关闭。假使该湾的容量为4,000,000公升，流入和流出河湾的水流速度均为40,000公升/天，若当前该河湾水中的污染物浓度为0.8mg/l，河水中污染物的浓度为0.5mg/l。要求1）建立湾中水污染状况的模型；2）计算30天后该河湾水的污染物浓度；3）该河湾水的污染物浓度是否能达到一个稳定值？
8、4）如将4,000,000mg污染物瞬间排入河水中，求铝厂必须关闭多长时间？5）列出并讨论影响河湾污染的模型中未考虑到的因素至少四种。（1）假设t时刻，河湾中污染物浓度为W，河水中污染物浓度为W1则浓度变化量等于流入的污染物减去流出的
