Устойчивость

> if steady(0,0,2,4.5,1,1,3,2)=1 then print(`(0,0)-Steady`):else print(`(0,0)-UnSteady`):fi;

> if steady(0,3,2,4.5,1,1,3,2)=1 then print(`(0,3)-Steady`):else print(`(0,3)-UnSteady`):fi;

> if steady(4.5,0,2,4.5,1,1,3,2)=1 then print(`(4.5,0)-Steady`):else print(`(4.5,0)-UnSteady`):fi;

> dfieldplot([D(x)(t)=f_11(x,y,t),D(y)(t)=f_22(x,y,t)],[x(t),y(t)],t=0..15,x=0..5,y=0..5,title=`{dx/dt=f_11(x,y),dy/dt=f_22(x,y)}`);

> F:=dsolve({D(x)(t)=f_11(x,y,t),D(y)(t)=f_22(x,y,t),x(0)=6,y(0)=5},{x(t),y(t)},type=numeric):

> g_1:=plots[odeplot](F,[t,x(t)],0..27,labels=[t,x],color=red):

> g_2:=plots[odeplot](F,[t,y(t)],0..27,labels=[t,y],color=blue):

> display({g_1,g_2},title=`First kind-red, second kind-blue`);