Competitive Species With a Parameter
We will be performing some matrix calculations. I am using the Linear Algebra, rather than the linalg, package.
with(LinearAlgebra):
Differential equations.
xdot := x*(2-x-a*y);
ydot := y*(3-y-2*a*x);
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USV4ZG90RigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2JS1GLTY5USJ4RihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8tRl9vNjNRMSZJbnZpc2libGVUaW1lcztGKEZib0Zlb0Znby9Gam9RJDBlbUYoL0ZdcEZfckZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYlLUZfbzYzUSIoRigvRmNvUSdwcmVmaXhGKC9GZm9GO0Znby9Gam9RLnRoaW5tYXRoc3BhY2VGKC9GXXBGanIvRl9wRjtGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2Jy1JI21uR0YlNjlRIjJGKEYwRjNGNi9GOkY4RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbi9Gam5RJ25vcm1hbEYoRlxvLUZfbzYzUSgmbWludXM7RihGYm9GZW9GZ28vRmpvUTBtZWRpdW1tYXRoc3BhY2VGKC9GXXBGanNGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHFGaHFGZnMtRiQ2JS1GLTY5USJhRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GW3ItRi02OVEieUYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvLUZfbzYzUSIpRigvRmNvUShwb3N0Zml4RihGaHJGZ29GaXIvRl1wUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYoRlxzRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxNyMtX0YpSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSV4ZG90R0YoKiZJInhHRigiIiIsKCIiI0ZldUZkdSEiIiomSSJhR0YoRmV1SSJ5R0YoRmV1Rmh1RmV1NyNGY3U=
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USV5ZG90RigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2JS1GLTY5USJ5RihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8tRl9vNjNRMSZJbnZpc2libGVUaW1lcztGKEZib0Zlb0Znby9Gam9RJDBlbUYoL0ZdcEZfckZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYlLUZfbzYzUSIoRigvRmNvUSdwcmVmaXhGKC9GZm9GO0Znby9Gam9RLnRoaW5tYXRoc3BhY2VGKC9GXXBGanIvRl9wRjtGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2Jy1JI21uR0YlNjlRIjNGKEYwRjNGNi9GOkY4RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbi9Gam5RJ25vcm1hbEYoRlxvLUZfbzYzUSgmbWludXM7RihGYm9GZW9GZ28vRmpvUTBtZWRpdW1tYXRoc3BhY2VGKC9GXXBGanNGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHFGaHFGZnMtRiQ2Jy1GYHM2OVEiMkYoRjBGM0Y2RmNzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZkc0Zcb0Zbci1GLTY5USJhRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG9GW3ItRi02OVEieEYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvLUZfbzYzUSIpRigvRmNvUShwb3N0Zml4RihGaHJGZ29GaXIvRl1wUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYoRlxzRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxNyMtX0YpSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSV5ZG90R0YoKiZJInlHRigiIiIsKCIiJEZodUZndSEiIiooIiIjRmh1SSJhR0YoRmh1SSJ4R0YoRmh1Rlt2Rmh1NyNGZnU=
Find equilibria. Note that three of the equilibria are independent of a, but the fourth equilibrium is in the first quadrant if and only if a < 2/3 or 3/4 < a.
solve({xdot=0, ydot=0}, {x, y});
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiktRiQ2JS1JI21vR0YlNjNRInxmckYoLyUlZm9ybUdRJ3ByZWZpeEYoLyUmZmVuY2VHUSV0cnVlRigvJSpzZXBhcmF0b3JHUSZmYWxzZUYoLyUnbHNwYWNlR1EudGhpbm1hdGhzcGFjZUYoLyUncnNwYWNlR0Y9LyUpc3RyZXRjaHlHRjcvJSpzeW1tZXRyaWNHRjovJShtYXhzaXplR1EpaW5maW5pdHlGKC8lKG1pbnNpemVHUSIxRigvJShsYXJnZW9wR0Y6LyUubW92YWJsZWxpbWl0c0dGOi8lJ2FjY2VudEdGOi8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUlc2l6ZUdRIzEyRigvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLUYkNiUtRiQ2JS1JI21pR0YlNjlRInhGKC8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdGVS8lJWJvbGRHRjovJSdpdGFsaWNHRjcvJSp1bmRlcmxpbmVHRjovJSpzdWJzY3JpcHRHRjovJSxzdXBlcnNjcmlwdEdGOi8lK2ZvcmVncm91bmRHRlgvJStiYWNrZ3JvdW5kR0Zlbi8lJ29wYXF1ZUdGOi8lK2V4ZWN1dGFibGVHRjovJSlyZWFkb25seUdGNy8lKWNvbXBvc2VkR0Y6LyUqY29udmVydGVkR0Y6LyUraW1zZWxlY3RlZEdGOi8lLHBsYWNlaG9sZGVyR0Y6LyUwZm9udF9zdHlsZV9uYW1lR0ZSLyUqbWF0aGNvbG9yR0ZYLyUvbWF0aGJhY2tncm91bmRHRmVuLyUrZm9udGZhbWlseUdGYG8vJSxtYXRodmFyaWFudEdRJ2l0YWxpY0YoLyUpbWF0aHNpemVHRlUtRi82M1EiPUYoL0YzUSZpbmZpeEYoL0Y2RjpGOC9GPFEvdGhpY2ttYXRoc3BhY2VGKC9GP0Zjci9GQUY6RkJGREZHRkpGTEZORlBGU0ZWRlktSSNtbkdGJTY5USIwRihGXm9GYW9GY28vRmZvRjpGZ29GaW9GW3BGXXBGX3BGYXBGY3BGZXBGZ3BGaXBGW3FGXXFGX3FGYXFGY3FGZXEvRmhxUSdub3JtYWxGKEZqcS1GLzYzUSIsRihGX3JGYXIvRjlGNy9GPFEkMGVtRigvRj9RM3Zlcnl0aGlja21hdGhzcGFjZUYoRmVyRkJGREZHRkpGTEZORlBGU0ZWRlktRiQ2JS1GW282OVEieUYoRl5vRmFvRmNvRmVvRmdvRmlvRltwRl1wRl9wRmFwRmNwRmVwRmdwRmlwRltxRl1xRl9xRmFxRmNxRmVxRmdxRmpxRlxyRmZyLUYvNjNRInxockYoL0YzUShwb3N0Zml4RihGNUY4RjsvRj9RMnZlcnl0aGlubWF0aHNwYWNlRihGQEZCRkRGR0ZKRkxGTkZQRlNGVkZZRl1zLUYkNiVGLi1GJDYlRmhuRl1zLUYkNiVGZ3NGXHItRmdyNjlRIjNGKEZeb0Zhb0Zjb0ZqckZnb0Zpb0ZbcEZdcEZfcEZhcEZjcEZlcEZncEZpcEZbcUZdcUZfcUZhcUZjcUZlcUZbc0ZqcUZqc0Zdcy1GJDYlRi4tRiQ2JUZlc0Zdcy1GJDYlRmpuRlxyLUZncjY5USIyRihGXm9GYW9GY29GanJGZ29GaW9GW3BGXXBGX3BGYXBGY3BGZXBGZ3BGaXBGW3FGXXFGX3FGYXFGY3FGZXFGW3NGanFGanNGXXMtRiQ2JUYuLUYkNiUtRiQ2JUZnc0Zcci1JJm1mcmFjR0YlNiotRiQ2Iy1GJDYmLUYvNjNRKiZ1bWludXMwO0YoRjJGYXJGOEZhcy9GP0Zic0ZlckZCRkRGR0ZKRkxGTkZQRlNGVkZZRmd0LUYvNjNRIitGKEZfckZhckY4L0Y8UTBtZWRpdW1tYXRoc3BhY2VGKC9GP0ZodkZlckZCRkRGR0ZKRkxGTkZQRlNGVkZZLUYkNiUtRmdyNjlRIjRGKEZeb0Zhb0Zjb0ZqckZnb0Zpb0ZbcEZdcEZfcEZhcEZjcEZlcEZncEZpcEZbcUZdcUZfcUZhcUZjcUZlcUZbc0ZqcS1GLzYzUTEmSW52aXNpYmxlVGltZXM7RihGX3JGYXJGOEZhc0ZjdkZlckZCRkRGR0ZKRkxGTkZQRlNGVkZZLUZbbzY5USJhRihGXm9GYW9GY29GZW9GZ29GaW9GW3BGXXBGX3BGYXBGY3BGZXBGZ3BGaXBGW3FGXXFGX3FGYXFGY3FGZXFGZ3FGanEtRiQ2Iy1GJDYmRmB2LUZncjY5RklGXm9GYW9GY29GanJGZ29GaW9GW3BGXXBGX3BGYXBGY3BGZXBGZ3BGaXBGW3FGXXFGX3FGYXFGY3FGZXFGW3NGanFGZHYtRiQ2JUZgdUZfdy1JJW1zdXBHRiU2JUZid0ZgdS8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRigvJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRmh4LyUpYmV2ZWxsZWRHRjpGVkZZRl1zLUYkNiVGam5GXHItRmp1NiotRiQ2Iy1GJDYmRmB2RmB1RmR2LUYkNiVGZ3RGX3dGYndGZXdGY3hGZnhGaXhGW3lGVkZZRmpzNyM2JjwkL0kieEdGKCIiIS9JInlHRihGXHo8JEZqeS9GXnoiIiQ8JEZdei9GW3oiIiM8JC9GXnoqJiwmRmF6ISIiKiYiIiUiIiJJImFHRihGXFtsRlxbbEZcW2wsJkZcW2xGaXoqJkZkekZcW2wpRl1bbEZkekZcW2xGXFtsRml6L0ZbeiomLCZGZHpGaXoqJkZhekZcW2xGXVtsRlxbbEZcW2xGXFtsRl5bbEZpeg==
fsolve(-2+3*a=0,a);
fsolve(-3+4*a=0,a);
fsolve(-1+2*a^2=0,a);
NiQtSSNtbkc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGKDY5US0wLjY2NjY2NjY2NjdGKC8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdRIzEyRigvJSVib2xkR1EmZmFsc2VGKC8lJ2l0YWxpY0dGNS8lKnVuZGVybGluZUdGNS8lKnN1YnNjcmlwdEdGNS8lLHN1cGVyc2NyaXB0R0Y1LyUrZm9yZWdyb3VuZEdRKlswLDAsMjU1XUYoLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV1GKC8lJ29wYXF1ZUdGNS8lK2V4ZWN1dGFibGVHRjUvJSlyZWFkb25seUdRJXRydWVGKC8lKWNvbXBvc2VkR0Y1LyUqY29udmVydGVkR0Y1LyUraW1zZWxlY3RlZEdGNS8lLHBsYWNlaG9sZGVyR0Y1LyUwZm9udF9zdHlsZV9uYW1lR1EqMkR+T3V0cHV0RigvJSptYXRoY29sb3JHRkAvJS9tYXRoYmFja2dyb3VuZEdGQy8lK2ZvbnRmYW1pbHlHRi8vJSxtYXRodmFyaWFudEdRJ25vcm1hbEYoLyUpbWF0aHNpemVHRjI3IyQiK25tbW1tISM1
NiQtSSNtbkc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGKDY5US0wLjc1MDAwMDAwMDBGKC8lJ2ZhbWlseUdRMFRpbWVzfk5ld35Sb21hbkYoLyUlc2l6ZUdRIzEyRigvJSVib2xkR1EmZmFsc2VGKC8lJ2l0YWxpY0dGNS8lKnVuZGVybGluZUdGNS8lKnN1YnNjcmlwdEdGNS8lLHN1cGVyc2NyaXB0R0Y1LyUrZm9yZWdyb3VuZEdRKlswLDAsMjU1XUYoLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV1GKC8lJ29wYXF1ZUdGNS8lK2V4ZWN1dGFibGVHRjUvJSlyZWFkb25seUdRJXRydWVGKC8lKWNvbXBvc2VkR0Y1LyUqY29udmVydGVkR0Y1LyUraW1zZWxlY3RlZEdGNS8lLHBsYWNlaG9sZGVyR0Y1LyUwZm9udF9zdHlsZV9uYW1lR1EqMkR+T3V0cHV0RigvJSptYXRoY29sb3JHRkAvJS9tYXRoYmFja2dyb3VuZEdGQy8lK2ZvbnRmYW1pbHlHRi8vJSxtYXRodmFyaWFudEdRJ25vcm1hbEYoLyUpbWF0aHNpemVHRjI3IyQiKysrKyt2ISM1
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtbkdGJTY5UTUmdW1pbnVzMDsuNzA3MTA2NzgxMkYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR0Y4LyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R1EldHJ1ZUYoLyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGQy8lL21hdGhiYWNrZ3JvdW5kR0ZGLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIixGKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRk0vJSdsc3BhY2VHUSQwZW1GKC8lJ3JzcGFjZUdRM3Zlcnl0aGlja21hdGhzcGFjZUYoLyUpc3RyZXRjaHlHRjgvJSpzeW1tZXRyaWNHRjgvJShtYXhzaXplR1EpaW5maW5pdHlGKC8lKG1pbnNpemVHUSIxRigvJShsYXJnZW9wR0Y4LyUubW92YWJsZWxpbWl0c0dGOC8lJ2FjY2VudEdGOC8lMGZvbnRfc3R5bGVfbmFtZUdGWC8lJXNpemVHRjUvJStmb3JlZ3JvdW5kR0ZDLyUrYmFja2dyb3VuZEdGRi1GLTY5US0wLjcwNzEwNjc4MTJGKEYwRjNGNkY5RjtGPUY/RkFGREZHRklGS0ZORlBGUkZURlZGWUZlbkZnbkZpbkZcbzcjNiQkISs3eTFycSEjNSQiKzd5MXJxRl5y
Find Jacobian.
J := <<diff(xdot,x) | diff(xdot,y)>, <diff(ydot,x) | diff(ydot,y)>>;
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USJKRigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2JS1GX282M1EiW0YoL0Zjb1EncHJlZml4RigvRmZvRjtGZ28vRmpvUS50aGlubWF0aHNwYWNlRigvRl1wRl9yL0ZfcEY7RmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiMtSSdtdGFibGVHRiU2JC1JJG10ckdGJTYkLUkkbXRkR0YlNiMtRiQ2Jy1JI21uR0YlNjlRIjJGKEYwRjNGNi9GOkY4RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbi9Gam5RJ25vcm1hbEYoRlxvLUZfbzYzUSgmbWludXM7RihGYm9GZW9GZ28vRmpvUTBtZWRpdW1tYXRoc3BhY2VGKC9GXXBGanNGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JUZfcy1GX282M1ExJkludmlzaWJsZVRpbWVzO0YoRmJvRmVvRmdvL0Zqb1EkMGVtRigvRl1wRmJ0Rl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYtNjlRInhGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0Zmcy1GJDYlLUYtNjlRImFGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcb0ZedC1GLTY5USJ5RihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8tRltzNiMtRiQ2JC1GX282M1EqJnVtaW51czA7RihGW3JGZW9GZ29GYXRGY3RGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JUZpdEZedEZkdC1GaHI2JC1GW3M2Iy1GJDYkRmN1LUYkNidGX3NGXnRGaXRGXnRGXHUtRltzNiMtRiQ2Jy1GYHM2OVEiM0YoRjBGM0Y2RmNzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZkc0Zcb0Zmcy1GJDYlRl9zRl50Rlx1RmZzLUYkNidGX3NGXnRGaXRGXnRGZHQtRl9vNjNRIl1GKC9GY29RKHBvc3RmaXhGKEZdckZnb0Zeci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGYXJGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHE3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JIkpHRigtSSdSVEFCTEVHRig2JSIqI0gxaDktSSdNQVRSSVhHRig2IzckNyQsKCIiIyIiIiomRmR4RmV4SSJ4R0YoRmV4ISIiKiZJImFHRihGZXhJInlHRihGZXhGaHgsJComRmp4RmV4Rmd4RmV4Rmh4NyQsJCooRmR4RmV4Rmp4RmV4Rlt5RmV4Rmh4LCgiIiRGZXgqJkZkeEZleEZbeUZleEZoeCooRmR4RmV4Rmp4RmV4Rmd4RmV4Rmh4SSdNYXRyaXhHNiQlKnByb3RlY3RlZEdGKjcjLUZleTYjL0kkJWlkR0YoRl14
Neither Survive Equilibrium. Source for all a > 0.
J0 := subs({y = 0, x = 0}, J);
DD := Determinant(J0); T := Trace(J0); T^2-4*DD;
Eigenvectors(J0);
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USNKMEYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiUtRl9vNjNRIltGKC9GY29RJ3ByZWZpeEYoL0Zmb0Y7RmdvL0Zqb1EudGhpbm1hdGhzcGFjZUYoL0ZdcEZfci9GX3BGO0ZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYjLUknbXRhYmxlR0YlNiQtSSRtdHJHRiU2JC1JJG10ZEdGJTYjLUkjbW5HRiU2OVEiMkYoRjBGM0Y2L0Y6RjhGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduL0ZqblEnbm9ybWFsRihGXG8tRltzNiMtRl5zNjlRIjBGKEYwRjNGNkZhc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GYnNGXG8tRmhyNiRGZHMtRltzNiMtRl5zNjlRIjNGKEYwRjNGNkZhc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GYnNGXG8tRl9vNjNRIl1GKC9GY29RKHBvc3RmaXhGKEZdckZnb0Zeci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGYXJGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHE3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JI0owR0YoLUknUlRBQkxFR0YoNiUiKltYOFkiLUknTUFUUklYR0YoNiM3JDckIiIjIiIhNyRGaXUiIiRJJ01hdHJpeEc2JCUqcHJvdGVjdGVkR0YqNyMtRlx2NiMvSSQlaWRHRihGYnU=
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USNEREYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUkjbW5HRiU2OVEiNkYoRjBGM0Y2L0Y6RjhGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduL0ZqblEnbm9ybWFsRihGXG83Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JI0RER0YoIiInNyNGZXI=
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USJURigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctSSNtbkdGJTY5USI1RihGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcbzcjLV9GKUksbXByaW50c2xhc2hHRig2JDcjPkkiVEdGKCIiJjcjRmVy
NiQtSSNtbkc2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGKDY5USIxRigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHRjUvJSp1bmRlcmxpbmVHRjUvJSpzdWJzY3JpcHRHRjUvJSxzdXBlcnNjcmlwdEdGNS8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjUvJStleGVjdXRhYmxlR0Y1LyUpcmVhZG9ubHlHUSV0cnVlRigvJSljb21wb3NlZEdGNS8lKmNvbnZlcnRlZEdGNS8lK2ltc2VsZWN0ZWRHRjUvJSxwbGFjZWhvbGRlckdGNS8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZALyUvbWF0aGJhY2tncm91bmRHRkMvJStmb250ZmFtaWx5R0YvLyUsbWF0aHZhcmlhbnRHUSdub3JtYWxGKC8lKW1hdGhzaXplR0YyNyMiIiI=
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtRiQ2JS1JI21vR0YlNjNRIltGKC8lJWZvcm1HUSdwcmVmaXhGKC8lJmZlbmNlR1EldHJ1ZUYoLyUqc2VwYXJhdG9yR1EmZmFsc2VGKC8lJ2xzcGFjZUdRLnRoaW5tYXRoc3BhY2VGKC8lJ3JzcGFjZUdGPS8lKXN0cmV0Y2h5R0Y3LyUqc3ltbWV0cmljR0Y6LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOi8lLm1vdmFibGVsaW1pdHNHRjovJSdhY2NlbnRHRjovJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lJXNpemVHUSMxMkYoLyUrZm9yZWdyb3VuZEdRKlswLDAsMjU1XUYoLyUrYmFja2dyb3VuZEdRLlsyNTUsMjU1LDI1NV1GKC1GJDYjLUknbXRhYmxlR0YlNiQtSSRtdHJHRiU2Iy1JJG10ZEdGJTYjLUkjbW5HRiU2OVEiMkYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR0ZVLyUlYm9sZEdGOi8lJ2l0YWxpY0dGOi8lKnVuZGVybGluZUdGOi8lKnN1YnNjcmlwdEdGOi8lLHN1cGVyc2NyaXB0R0Y6LyUrZm9yZWdyb3VuZEdGWC8lK2JhY2tncm91bmRHRmVuLyUnb3BhcXVlR0Y6LyUrZXhlY3V0YWJsZUdGOi8lKXJlYWRvbmx5R0Y3LyUpY29tcG9zZWRHRjovJSpjb252ZXJ0ZWRHRjovJStpbXNlbGVjdGVkR0Y6LyUscGxhY2Vob2xkZXJHRjovJTBmb250X3N0eWxlX25hbWVHRlIvJSptYXRoY29sb3JHRlgvJS9tYXRoYmFja2dyb3VuZEdGZW4vJStmb250ZmFtaWx5R0Znby8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRigvJSltYXRoc2l6ZUdGVS1GXG82Iy1GX282Iy1GYm82OVEiM0YoRmVvRmhvRmpvRlxwRl5wRmBwRmJwRmRwRmZwRmhwRmpwRlxxRl5xRmBxRmJxRmRxRmZxRmhxRmpxRlxyRl5yRmFyLUYvNjNRIl1GKC9GM1EocG9zdGZpeEYoRjVGOEY7L0Y/UTJ2ZXJ5dGhpbm1hdGhzcGFjZUYoRkBGQkZERkdGSkZMRk5GUEZTRlZGWS1GLzYzUSIsRigvRjNRJmluZml4RigvRjZGOi9GOUY3L0Y8USQwZW1GKC9GP1EzdmVyeXRoaWNrbWF0aHNwYWNlRigvRkFGOkZCRkRGR0ZKRkxGTkZQRlNGVkZZLUYkNiVGLi1GJDYjLUZpbjYkLUZcbzYkLUZfbzYjLUZibzY5RklGZW9GaG9Gam9GXHBGXnBGYHBGYnBGZHBGZnBGaHBGanBGXHFGXnFGYHFGYnFGZHFGZnFGaHFGanFGXHJGXnJGYXItRl9vNiMtRmJvNjlRIjBGKEZlb0Zob0Zqb0ZccEZecEZgcEZicEZkcEZmcEZocEZqcEZccUZecUZgcUZicUZkcUZmcUZocUZqcUZcckZeckZhci1GXG82JEZpdEZldEZqcjcjLV9GKUksbXByaW50c2xhc2hHRig2JDckLUknUlRBQkxFR0YoNiUiKndxWVkiLUknTUFUUklYR0YoNiM3JDcjIiIjNyMiIiQmSSdWZWN0b3JHNiQlKnByb3RlY3RlZEdGKjYjSSdjb2x1bW5HRigtRmd1NiUiKisnPnE5LUZbdjYjNyQ3JCIiIiIiITckRmB3Rl93SSdNYXRyaXhHRmR2NyQtRmJ2NiMvSSQlaWRHRihGaXUtRmJ3NiMvRmd3Rmp2
Only Species Y Survives Equilibrium. Saddle if a < 2/3; sink if 2/3 < a; never a spiral sink although at a = 5/3 the sink has only one straight-line solution..
J1 := subs({x = 0, y = 3}, J);
DD := Determinant(J1); T := Trace(J1); disc := simplify(T^2-4*DD);
E1 := Eigenvectors(J1);
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USNKMUYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiUtRl9vNjNRIltGKC9GY29RJ3ByZWZpeEYoL0Zmb0Y7RmdvL0Zqb1EudGhpbm1hdGhzcGFjZUYoL0ZdcEZfci9GX3BGO0ZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYjLUknbXRhYmxlR0YlNiQtSSRtdHJHRiU2JC1JJG10ZEdGJTYjLUYkNiUtSSNtbkdGJTY5USIyRihGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcby1GX282M1EoJm1pbnVzO0YoRmJvRmVvRmdvL0Zqb1EwbWVkaXVtbWF0aHNwYWNlRigvRl1wRmpzRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRmBzNjlRIjNGKEYwRjNGNkZjc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GZHNGXG8tRl9vNjNRMSZJbnZpc2libGVUaW1lcztGKEZib0Zlb0Znby9Gam9RJDBlbUYoL0ZdcEZldEZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GLTY5USJhRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8tRltzNiMtRmBzNjlRIjBGKEYwRjNGNkZjc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GZHNGXG8tRmhyNiQtRltzNiMtRiQ2JC1GX282M1EqJnVtaW51czA7RihGW3JGZW9GZ29GZHRGZnRGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JS1GYHM2OVEiNkYoRjBGM0Y2RmNzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZkc0Zcb0ZhdEZndC1GW3M2Iy1GYHM2OVEpJm1pbnVzOzNGKEYwRjNGNkZjc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GZHNGXG8tRl9vNjNRIl1GKC9GY29RKHBvc3RmaXhGKEZdckZnb0Zeci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGYXJGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHE3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JI0oxR0YoLUknUlRBQkxFR0YoNiUiKj95JXo5LUknTUFUUklYR0YoNiM3JDckLCYiIiMiIiIqJiIiJEZceEkiYUdGKEZceCEiIiIiITckLCQqJiIiJ0ZceEZfeEZceEZgeCEiJEknTWF0cml4RzYkJSpwcm90ZWN0ZWRHRio3Iy1GZ3g2Iy9JJCVpZEdGKEZkdw==
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USNEREYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiYtRl9vNjNRKiZ1bWludXMwO0YoL0Zjb1EncHJlZml4RihGZW9GZ28vRmpvUSQwZW1GKC9GXXBGXnJGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtSSNtbkdGJTY5USI2RihGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcby1GX282M1EiK0YoRmJvRmVvRmdvL0Zqb1EwbWVkaXVtbWF0aHNwYWNlRigvRl1wRltzRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRmFyNjlRIjlGKEYwRjNGNkZkckY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GZXJGXG8tRl9vNjNRMSZJbnZpc2libGVUaW1lcztGKEZib0Zlb0Znb0ZdckZfckZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GLTY5USJhRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG83Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JI0RER0YoLCYiIichIiIqJiIiKiIiIkkiYUdGKEZldEZldDcjRmB0
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USJURigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2Ji1GX282M1EqJnVtaW51czA7RigvRmNvUSdwcmVmaXhGKEZlb0Znby9Gam9RJDBlbUYoL0ZdcEZeckZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1JI21uR0YlNjlGZ3BGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcby1GX282M1EoJm1pbnVzO0YoRmJvRmVvRmdvL0Zqb1EwbWVkaXVtbWF0aHNwYWNlRigvRl1wRmpyRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRmFyNjlRIjNGKEYwRjNGNkZjckY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GZHJGXG8tRl9vNjNRMSZJbnZpc2libGVUaW1lcztGKEZib0Zlb0Znb0ZdckZfckZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GLTY5USJhRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG83Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JIlRHRigsJiIiIiEiIiomIiIkRmB0SSJhR0YoRmB0RmF0NyNGX3Q=
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USVkaXNjRigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2Jy1JI21uR0YlNjlRIzI1RihGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcby1GX282M1EoJm1pbnVzO0YoRmJvRmVvRmdvL0Zqb1EwbWVkaXVtbWF0aHNwYWNlRigvRl1wRmNyRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRmlxNjlRIzMwRihGMEYzRjZGXHJGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRl1yRlxvLUZfbzYzUTEmSW52aXNpYmxlVGltZXM7RihGYm9GZW9GZ28vRmpvUSQwZW1GKC9GXXBGXnNGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRi02OVEiYUYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvLUZfbzYzUSIrRihGYm9GZW9GZ29GYnJGZHJGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JS1GaXE2OVEiOUYoRjBGM0Y2RlxyRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZdckZcb0Zqci1JJW1zdXBHRiU2JUZgcy1GaXE2OVEiMkYoRjBGM0Y2RlxyRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZdckZcby8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRig3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JJWRpc2NHRigsKCIjRCIiIiomIiNJRl51SSJhR0YoRl51ISIiKiYiIipGXnUpRmF1IiIjRl51Rl51NyNGXHU=
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USNFMUYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiUtRiQ2JS1GX282M1EiW0YoL0Zjb1EncHJlZml4RigvRmZvRjtGZ28vRmpvUS50aGlubWF0aHNwYWNlRigvRl1wRmFyL0ZfcEY7RmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiMtSSdtdGFibGVHRiU2JC1JJG10ckdGJTYjLUkkbXRkR0YlNiMtSSNtbkdGJTY5USkmbWludXM7M0YoRjBGM0Y2L0Y6RjhGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduL0ZqblEnbm9ybWFsRihGXG8tRmpyNiMtRl1zNiMtRiQ2JS1GYHM2OVEiMkYoRjBGM0Y2RmNzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZkc0Zcby1GX282M1EoJm1pbnVzO0YoRmJvRmVvRmdvL0Zqb1EwbWVkaXVtbWF0aHNwYWNlRigvRl1wRmN0Rl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRmBzNjlRIjNGKEYwRjNGNkZjc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GZHNGXG8tRl9vNjNRMSZJbnZpc2libGVUaW1lcztGKEZib0Zlb0Znby9Gam9RJDBlbUYoL0ZdcEZedUZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GLTY5USJhRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8tRl9vNjNRIl1GKC9GY29RKHBvc3RmaXhGKEZfckZnb0Zgci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGY3JGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRl9vNjNRIixGKEZib0Zlby9GaG9GO0ZddS9GXXBRM3Zlcnl0aGlja21hdGhzcGFjZUYoRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiVGanEtRiQ2Iy1GZ3I2JC1GanI2JC1GXXM2Iy1GYHM2OVEiMEYoRjBGM0Y2RmNzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZkc0Zcby1GXXM2Iy1GJDYlLUkmbWZyYWNHRiU2Ki1GYHM2OUZncEYwRjNGNkZjc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GZHNGXG8tRmBzNjlRIjZGKEYwRjNGNkZjc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GZHNGXG8vJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRl54LyUpYmV2ZWxsZWRHRjhGYnFGZHFGanQtRmJ3NiotRiQ2Iy1GJDYmLUZfbzYzUSomdW1pbnVzMDtGKEZdckZlb0Znb0ZddUZfdUZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GYHM2OVEiNUYoRjBGM0Y2RmNzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZkc0Zcby1GX282M1EiK0YoRmJvRmVvRmdvRmJ0RmR0Rl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxRmV0LUYkNiNGYHVGaXdGXHhGX3hGYXhGYnFGZHEtRmpyNiQtRl1zNiNGZHdGZnlGY3U3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JI0UxR0YoNiQtSSdSVEFCTEVHRig2JSIqd3MoejktSSdNQVRSSVhHRig2IzckNyMhIiQ3IywmIiIjIiIiKiYiIiRGXltsSSJhR0YoRl5bbCEiIiZJJ1ZlY3Rvckc2JCUqcHJvdGVjdGVkR0YqNiNJJ2NvbHVtbkdGKC1GYno2JSIqK0FYWyItRmZ6NiM3JDckIiIhLCQqKCNGXltsIiInRl5bbCwmIiImRmJbbEZfW2xGXltsRl5bbEZhW2xGYltsRl5bbDckRl5bbEZeW2xJJ01hdHJpeEdGZVtsNyQtRmNbbDYjL0kkJWlkR0YoRmR6LUZoXGw2Iy9GXV1sRltcbA==
solve(DD=0);
solve(disc = 0);
NiQtSSZtZnJhY0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGKDYqLUkjbW5HRiU2OVEiMkYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR0Y4LyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R1EldHJ1ZUYoLyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGQy8lL21hdGhiYWNrZ3JvdW5kR0ZGLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRigvJSltYXRoc2l6ZUdGNS1GLTY5USIzRihGMEYzRjZGOUY7Rj1GP0ZBRkRGR0ZJRktGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8vJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRmZvLyUpYmV2ZWxsZWRHRjgvJStmb3JlZ3JvdW5kR0ZDLyUrYmFja2dyb3VuZEdGRjcjIyIiIyIiJA==
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSZtZnJhY0dGJTYqLUkjbW5HRiU2OVEiNUYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR0Y7LyUqdW5kZXJsaW5lR0Y7LyUqc3Vic2NyaXB0R0Y7LyUsc3VwZXJzY3JpcHRHRjsvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y7LyUrZXhlY3V0YWJsZUdGOy8lKXJlYWRvbmx5R1EldHJ1ZUYoLyUpY29tcG9zZWRHRjsvJSpjb252ZXJ0ZWRHRjsvJStpbXNlbGVjdGVkR0Y7LyUscGxhY2Vob2xkZXJHRjsvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRi8lL21hdGhiYWNrZ3JvdW5kR0ZJLyUrZm9udGZhbWlseUdGNS8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRigvJSltYXRoc2l6ZUdGOC1GMDY5USIzRihGM0Y2RjlGPEY+RkBGQkZERkdGSkZMRk5GUUZTRlVGV0ZZRmZuRmhuRmpuRlxvRl9vLyUubGluZXRoaWNrbmVzc0dRIjFGKC8lK2Rlbm9tYWxpZ25HUSdjZW50ZXJGKC8lKW51bWFsaWduR0Zpby8lKWJldmVsbGVkR0Y7LyUrZm9yZWdyb3VuZEdGRi8lK2JhY2tncm91bmRHRkktSSNtb0dGJTYzUSIsRigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y7LyUqc2VwYXJhdG9yR0ZQLyUnbHNwYWNlR1EkMGVtRigvJSdyc3BhY2VHUTN2ZXJ5dGhpY2ttYXRoc3BhY2VGKC8lKXN0cmV0Y2h5R0Y7LyUqc3ltbWV0cmljR0Y7LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOy8lLm1vdmFibGVsaW1pdHNHRjsvJSdhY2NlbnRHRjsvJTBmb250X3N0eWxlX25hbWVHRmVuLyUlc2l6ZUdGOEZecEZgcEYsNyM2JCMiIiYiIiRGaXI=
Only Species X Survives Equilibrium. Saddle if a < 3/4; sink if 3/4 < a; never a spiral sink although a sink with a single straight line solution at a = 5/4..
J2 := subs({y = 0, x = 2}, J);
DD := Determinant(J2); T := Trace(J2); disc := simplify(T^2-4*DD);
E2 := Eigenvectors(J2);
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USNKMkYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiUtRl9vNjNRIltGKC9GY29RJ3ByZWZpeEYoL0Zmb0Y7RmdvL0Zqb1EudGhpbm1hdGhzcGFjZUYoL0ZdcEZfci9GX3BGO0ZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYjLUknbXRhYmxlR0YlNiQtSSRtdHJHRiU2JC1JJG10ZEdGJTYjLUkjbW5HRiU2OVEpJm1pbnVzOzJGKEYwRjNGNi9GOkY4RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbi9Gam5RJ25vcm1hbEYoRlxvLUZbczYjLUYkNiQtRl9vNjNRKiZ1bWludXMwO0YoRltyRmVvRmdvL0Zqb1EkMGVtRigvRl1wRlx0Rl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRl5zNjlRIjJGKEYwRjNGNkZhc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GYnNGXG8tRl9vNjNRMSZJbnZpc2libGVUaW1lcztGKEZib0Zlb0Znb0ZbdEZddEZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GLTY5USJhRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8tRmhyNiQtRltzNiMtRl5zNjlRIjBGKEYwRjNGNkZhc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GYnNGXG8tRltzNiMtRiQ2JS1GXnM2OVEiM0YoRjBGM0Y2RmFzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZic0Zcby1GX282M1EoJm1pbnVzO0YoRmJvRmVvRmdvL0Zqb1EwbWVkaXVtbWF0aHNwYWNlRigvRl1wRlt2Rl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRl5zNjlRIjRGKEYwRjNGNkZhc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GYnNGXG9GY3RGZnQtRl9vNjNRIl1GKC9GY29RKHBvc3RmaXhGKEZdckZnb0Zeci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGYXJGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHE3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JI0oyR0YoLUknUlRBQkxFR0YoNiUiKk9XX1siLUknTUFUUklYR0YoNiM3JDckISIjLCQqJiIiIyIiIkkiYUdGKEZeeCEiIjckIiIhLCYiIiRGXngqJiIiJUZeeEZfeEZeeEZgeEknTWF0cml4RzYkJSpwcm90ZWN0ZWRHRio3Iy1GZ3g2Iy9JJCVpZEdGKEZkdw==
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USNEREYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiYtRl9vNjNRKiZ1bWludXMwO0YoL0Zjb1EncHJlZml4RihGZW9GZ28vRmpvUSQwZW1GKC9GXXBGXnJGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtSSNtbkdGJTY5USI2RihGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcby1GX282M1EiK0YoRmJvRmVvRmdvL0Zqb1EwbWVkaXVtbWF0aHNwYWNlRigvRl1wRltzRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRmFyNjlRIjhGKEYwRjNGNkZkckY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GZXJGXG8tRl9vNjNRMSZJbnZpc2libGVUaW1lcztGKEZib0Zlb0Znb0ZdckZfckZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GLTY5USJhRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG83Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JI0RER0YoLCYiIichIiIqJiIiKSIiIkkiYUdGKEZldEZldDcjRmB0
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USJURigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2JS1JI21uR0YlNjlGZ3BGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcby1GX282M1EoJm1pbnVzO0YoRmJvRmVvRmdvL0Zqb1EwbWVkaXVtbWF0aHNwYWNlRigvRl1wRmJyRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRmlxNjlRIjRGKEYwRjNGNkZbckY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GXHJGXG8tRl9vNjNRMSZJbnZpc2libGVUaW1lcztGKEZib0Zlb0Znby9Gam9RJDBlbUYoL0ZdcEZdc0ZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GLTY5USJhRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG83Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JIlRHRigsJiIiIkZbdComIiIlRlt0SSJhR0YoRlt0ISIiNyNGanM=
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USVkaXNjRigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2Jy1JI21uR0YlNjlRIzI1RihGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcby1GX282M1EoJm1pbnVzO0YoRmJvRmVvRmdvL0Zqb1EwbWVkaXVtbWF0aHNwYWNlRigvRl1wRmNyRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRmlxNjlRIzQwRihGMEYzRjZGXHJGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRl1yRlxvLUZfbzYzUTEmSW52aXNpYmxlVGltZXM7RihGYm9GZW9GZ28vRmpvUSQwZW1GKC9GXXBGXnNGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRi02OVEiYUYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvLUZfbzYzUSIrRihGYm9GZW9GZ29GYnJGZHJGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JS1GaXE2OVEjMTZGKEYwRjNGNkZcckY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GXXJGXG9GanItSSVtc3VwR0YlNiVGYHMtRmlxNjlRIjJGKEYwRjNGNkZcckY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GXXJGXG8vJTFzdXBlcnNjcmlwdHNoaWZ0R1EiMEYoNyMtX0YpSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSVkaXNjR0YoLCgiI0QiIiIqJiIjU0ZedUkiYUdGKEZedSEiIiomIiM7Rl51KUZhdSIiI0ZedUZedTcjRlx1
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USNFMkYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiUtRiQ2JS1GX282M1EiW0YoL0Zjb1EncHJlZml4RigvRmZvRjtGZ28vRmpvUS50aGlubWF0aHNwYWNlRigvRl1wRmFyL0ZfcEY7RmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiMtSSdtdGFibGVHRiU2JC1JJG10ckdGJTYjLUkkbXRkR0YlNiMtRiQ2JS1JI21uR0YlNjlRIjNGKEYwRjNGNi9GOkY4RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbi9Gam5RJ25vcm1hbEYoRlxvLUZfbzYzUSgmbWludXM7RihGYm9GZW9GZ28vRmpvUTBtZWRpdW1tYXRoc3BhY2VGKC9GXXBGXHRGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JS1GYnM2OVEiNEYoRjBGM0Y2RmVzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZmc0Zcby1GX282M1ExJkludmlzaWJsZVRpbWVzO0YoRmJvRmVvRmdvL0Zqb1EkMGVtRigvRl1wRmd0Rl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYtNjlRImFGKEYwRjNGNkY5RjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZpbkZcby1GanI2Iy1GXXM2Iy1GYnM2OVEpJm1pbnVzOzJGKEYwRjNGNkZlc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GZnNGXG8tRl9vNjNRIl1GKC9GY29RKHBvc3RmaXhGKEZfckZnb0Zgci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGY3JGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRl9vNjNRIixGKEZib0Zlby9GaG9GO0ZmdC9GXXBRM3Zlcnl0aGlja21hdGhzcGFjZUYoRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiVGanEtRiQ2Iy1GZ3I2JC1GanI2JC1GXXM2Iy1JJm1mcmFjR0YlNiotRiQ2JS1GYnM2OVEiMkYoRjBGM0Y2RmVzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZmc0Zcb0ZjdEZpdC1GJDYjLUYkNiYtRl9vNjNRKiZ1bWludXMwO0YoRl1yRmVvRmdvRmZ0Rmh0Rl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUZiczY5USI1RihGMEYzRjZGZXNGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmZzRlxvLUZfbzYzUSIrRihGYm9GZW9GZ29GW3RGXXRGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHFGXnQvJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRmR4LyUpYmV2ZWxsZWRHRjhGYnFGZHEtRl1zNiMtRmJzNjlGZ3BGMEYzRjZGZXNGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmZzRlxvLUZqcjYkRml4LUZdczYjLUZiczY5USIwRihGMEYzRjZGZXNGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmZzRlxvRmN1NyMtX0YpSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSNFMkdGKDYkLUknUlRBQkxFR0YoNiUiKl9xRFgiLUknTUFUUklYR0YoNiM3JDcjLCYiIiQiIiIqJiIiJUZoekkiYUdGKEZoeiEiIjcjISIjJkknVmVjdG9yRzYkJSpwcm90ZWN0ZWRHRio2I0knY29sdW1uR0YoLUZeejYlIipvJipbWCItRmJ6NiM3JDckLCQqKCIiI0ZoekZbW2xGaHosJiIiJkZcW2xGaXpGaHpGXFtsRmh6Rmh6NyRGaHoiIiFJJ01hdHJpeEdGYVtsNyQtRl9bbDYjL0kkJWlkR0YoRmB6LUZjXGw2Iy9GaFxsRmdbbA==
solve(DD=0);
solve(T=0);
solve(disc = 0);
NiQtSSZtZnJhY0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGKDYqLUkjbW5HRiU2OVEiM0YoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR0Y4LyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R1EldHJ1ZUYoLyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGQy8lL21hdGhiYWNrZ3JvdW5kR0ZGLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRigvJSltYXRoc2l6ZUdGNS1GLTY5USI0RihGMEYzRjZGOUY7Rj1GP0ZBRkRGR0ZJRktGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8vJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRmZvLyUpYmV2ZWxsZWRHRjgvJStmb3JlZ3JvdW5kR0ZDLyUrYmFja2dyb3VuZEdGRjcjIyIiJCIiJQ==
NiQtSSZtZnJhY0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGKDYqLUkjbW5HRiU2OVEiMUYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR0Y4LyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R1EldHJ1ZUYoLyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGQy8lL21hdGhiYWNrZ3JvdW5kR0ZGLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRigvJSltYXRoc2l6ZUdGNS1GLTY5USI0RihGMEYzRjZGOUY7Rj1GP0ZBRkRGR0ZJRktGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8vJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRmZvLyUpYmV2ZWxsZWRHRjgvJStmb3JlZ3JvdW5kR0ZDLyUrYmFja2dyb3VuZEdGRjcjIyIiIiIiJQ==
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSZtZnJhY0dGJTYqLUkjbW5HRiU2OVEiNUYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR0Y7LyUqdW5kZXJsaW5lR0Y7LyUqc3Vic2NyaXB0R0Y7LyUsc3VwZXJzY3JpcHRHRjsvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y7LyUrZXhlY3V0YWJsZUdGOy8lKXJlYWRvbmx5R1EldHJ1ZUYoLyUpY29tcG9zZWRHRjsvJSpjb252ZXJ0ZWRHRjsvJStpbXNlbGVjdGVkR0Y7LyUscGxhY2Vob2xkZXJHRjsvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRi8lL21hdGhiYWNrZ3JvdW5kR0ZJLyUrZm9udGZhbWlseUdGNS8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRigvJSltYXRoc2l6ZUdGOC1GMDY5USI0RihGM0Y2RjlGPEY+RkBGQkZERkdGSkZMRk5GUUZTRlVGV0ZZRmZuRmhuRmpuRlxvRl9vLyUubGluZXRoaWNrbmVzc0dRIjFGKC8lK2Rlbm9tYWxpZ25HUSdjZW50ZXJGKC8lKW51bWFsaWduR0Zpby8lKWJldmVsbGVkR0Y7LyUrZm9yZWdyb3VuZEdGRi8lK2JhY2tncm91bmRHRkktSSNtb0dGJTYzUSIsRigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y7LyUqc2VwYXJhdG9yR0ZQLyUnbHNwYWNlR1EkMGVtRigvJSdyc3BhY2VHUTN2ZXJ5dGhpY2ttYXRoc3BhY2VGKC8lKXN0cmV0Y2h5R0Y7LyUqc3ltbWV0cmljR0Y7LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOy8lLm1vdmFibGVsaW1pdHNHRjsvJSdhY2NlbnRHRjsvJTBmb250X3N0eWxlX25hbWVHRmVuLyUlc2l6ZUdGOEZecEZgcEYsNyM2JCMiIiYiIiVGaXI=
Both Species Survive Equilibrium. Both coordinates are positive if and only if a < 2/3 (sink) or 3/4 < a (saddle).
eq := {x = (-2+3*a)/(-1+2*a^2), y = (-3+4*a)/(-1+2*a^2)}:
J3 := subs(eq, J);
DD := Determinant(J3); T := simplify(Trace(J3)); disc := simplify(T^2-4*DD);
E3 := simplify([Eigenvectors(J3)]):
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USNKM0YoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiUtRl9vNjNRIltGKC9GY29RJ3ByZWZpeEYoL0Zmb0Y7RmdvL0Zqb1EudGhpbm1hdGhzcGFjZUYoL0ZdcEZfci9GX3BGO0ZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYjLUknbXRhYmxlR0YlNiQtSSRtdHJHRiU2JC1JJG10ZEdGJTYjLUYkNictSSNtbkdGJTY5USIyRihGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcby1GX282M1EoJm1pbnVzO0YoRmJvRmVvRmdvL0Zqb1EwbWVkaXVtbWF0aHNwYWNlRigvRl1wRmpzRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUkmbWZyYWNHRiU2Ki1GJDYlRl9zLUZfbzYzUTEmSW52aXNpYmxlVGltZXM7RihGYm9GZW9GZ28vRmpvUSQwZW1GKC9GXXBGZXRGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JS1GX282M1EiKEYoRltyRl1yRmdvRl5yRmByRmFyRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiYtRl9vNjNRKiZ1bWludXMwO0YoRltyRmVvRmdvRmR0RmZ0Rl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxRl9zLUZfbzYzUSIrRihGYm9GZW9GZ29GaXNGW3RGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2JS1GYHM2OVEiM0YoRjBGM0Y2RmNzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZkc0Zcb0ZhdC1GLTY5USJhRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8tRl9vNjNRIilGKC9GY29RKHBvc3RmaXhGKEZdckZnb0Zeci9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGYXJGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRiQ2Iy1GJDYmRl51LUZgczY5RmdwRjBGM0Y2RmNzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZkc0Zcb0ZhdS1GJDYlRl9zRmF0LUklbXN1cEdGJTYlRml1Rl9zLyUxc3VwZXJzY3JpcHRzaGlmdEdRIjBGKC8lLmxpbmV0aGlja25lc3NHUSIxRigvJStkZW5vbWFsaWduR1EnY2VudGVyRigvJSludW1hbGlnbkdGZncvJSliZXZlbGxlZEdGOEZicUZkcUZmcy1GXXQ2Ki1GJDYlRml1RmF0LUYkNiVGaXQtRiQ2JkZedUZmdUZhdS1GJDYlLUZgczY5USI0RihGMEYzRjZGY3NGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmRzRlxvRmF0Rml1Rlx2RmN2RmF3RmR3Rmd3Rml3RmJxRmRxLUZbczYjLUYkNiRGXnUtRl10NiotRiQ2JUZpdUZhdEZndEZjdkZhd0Zkd0Znd0Zpd0ZicUZkcS1GaHI2JC1GW3M2Iy1GJDYkRl51LUZddDYqLUYkNidGX3NGYXRGaXVGYXRGX3hGY3ZGYXdGZHdGZ3dGaXdGYnFGZHEtRltzNiMtRiQ2J0ZmdUZmcy1GXXQ2Ki1GJDYlRl9zRmF0Rl94RmN2RmF3RmR3Rmd3Rml3RmJxRmRxRmZzLUZddDYqLUYkNidGX3NGYXRGaXVGYXRGZ3RGY3ZGYXdGZHdGZ3dGaXdGYnFGZHEtRl9vNjNRIl1GKEZfdkZdckZnb0ZeckZhdkZhckZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcTcjLV9GKUksbXByaW50c2xhc2hHRig2JDcjPkkjSjNHRigtSSdSVEFCTEVHRig2JSIqX3I6WSItSSdNQVRSSVhHRig2IzckNyQsKCIiIyIiIiooRltcbEZcXGwsJkZbXGwhIiIqJiIiJEZcXGxJImFHRihGXFxsRlxcbEZcXGwsJkZcXGxGX1xsKiZGW1xsRlxcbClGYlxsRltcbEZcXGxGXFxsRl9cbEZfXGwqKEZiXGxGXFxsLCZGYVxsRl9cbComIiIlRlxcbEZiXGxGXFxsRlxcbEZcXGxGY1xsRl9cbEZfXGwsJCooRmJcbEZcXGxGXlxsRlxcbEZjXGxGX1xsRl9cbDckLCQqKkZbXGxGXFxsRmJcbEZcXGxGZ1xsRlxcbEZjXGxGX1xsRl9cbCwoRmFcbEZcXGwqKEZbXGxGXFxsRmdcbEZcXGxGY1xsRl9cbEZfXGwqKkZbXGxGXFxsRmJcbEZcXGxGXlxsRlxcbEZjXGxGX1xsRl9cbEknTWF0cml4RzYkJSpwcm90ZWN0ZWRHRio3Iy1GYl1sNiMvSSQlaWRHRihGZFts
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USNEREYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR1EldHJ1ZUYoLyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R0Y7LyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRC8lL21hdGhiYWNrZ3JvdW5kR0ZHLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1EnaXRhbGljRigvJSltYXRoc2l6ZUdGNS1JI21vR0YlNjNRIzo9RigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y4LyUqc2VwYXJhdG9yR0Y4LyUnbHNwYWNlR1EvdGhpY2ttYXRoc3BhY2VGKC8lJ3JzcGFjZUdGW3AvJSlzdHJldGNoeUdGOC8lKnN5bW1ldHJpY0dGOC8lKG1heHNpemVHUSlpbmZpbml0eUYoLyUobWluc2l6ZUdRIjFGKC8lKGxhcmdlb3BHRjgvJS5tb3ZhYmxlbGltaXRzR0Y4LyUnYWNjZW50R0Y4LyUwZm9udF9zdHlsZV9uYW1lR0ZYLyUlc2l6ZUdGNS8lK2ZvcmVncm91bmRHRkQvJStiYWNrZ3JvdW5kR0ZHLUYkNiQtRl9vNjNRKiZ1bWludXMwO0YoL0Zjb1EncHJlZml4RihGZW9GZ28vRmpvUSQwZW1GKC9GXXBGXnJGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtSSZtZnJhY0dGJTYqLUYkNiUtRiQ2JS1GX282M1EiKEYoRltyL0Zmb0Y7RmdvL0Zqb1EudGhpbm1hdGhzcGFjZUYoL0ZdcEZccy9GX3BGO0ZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYmRmhxLUkjbW5HRiU2OVEiMkYoRjBGM0Y2L0Y6RjhGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduL0ZqblEnbm9ybWFsRihGXG8tRl9vNjNRIitGKEZib0Zlb0Znby9Gam9RMG1lZGl1bW1hdGhzcGFjZUYoL0ZdcEZcdEZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYlLUZiczY5USIzRihGMEYzRjZGZXNGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmZzRlxvLUZfbzYzUTEmSW52aXNpYmxlVGltZXM7RihGYm9GZW9GZ29GXXJGX3JGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRi02OVEiYUYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvLUZfbzYzUSIpRigvRmNvUShwb3N0Zml4RihGanJGZ29GW3MvRl1wUTJ2ZXJ5dGhpbm1hdGhzcGFjZUYoRl5zRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxRmN0LUYkNiVGZ3ItRiQ2JkZocUZgdEZocy1GJDYlLUZiczY5USI0RihGMEYzRjZGZXNGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmZzRlxvRmN0RmZ0Rml0LUYkNiMtRiQ2JkZocS1GYnM2OUZncEYwRjNGNkZlc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GZnNGXG9GaHMtRiQ2JUZhc0ZjdC1JJW1zdXBHRiU2JUZmdEZhcy8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRigvJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRlx3LyUpYmV2ZWxsZWRHRjhGYnFGZHE3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JI0RER0YoLCQqKCwmIiIjISIiKiYiIiQiIiJJImFHRihGYHhGYHhGYHgsJkZfeEZdeComIiIlRmB4RmF4RmB4RmB4RmB4LCZGYHhGXXgqJkZceEZgeClGYXhGXHhGYHhGYHhGXXhGXXg3I0Zpdw==
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USJURigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctRiQ2JC1GX282M1EqJnVtaW51czA7RigvRmNvUSdwcmVmaXhGKEZlb0Znby9Gam9RJDBlbUYoL0ZdcEZeckZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1JJm1mcmFjR0YlNiotRiQ2Iy1GJDYmRmhxLUkjbW5HRiU2OVEiNUYoRjBGM0Y2L0Y6RjhGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduL0ZqblEnbm9ybWFsRihGXG8tRl9vNjNRIitGKEZib0Zlb0Znby9Gam9RMG1lZGl1bW1hdGhzcGFjZUYoL0ZdcEZic0ZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYlLUZocjY5USI3RihGMEYzRjZGW3NGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRlxzRlxvLUZfbzYzUTEmSW52aXNpYmxlVGltZXM7RihGYm9GZW9GZ29GXXJGX3JGXnBGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHEtRi02OVEiYUYoRjBGM0Y2RjlGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmluRlxvLUYkNiMtRiQ2JkZocS1GaHI2OUZncEYwRjNGNkZbc0Y8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GXHNGXG9GXnMtRiQ2JS1GaHI2OVEiMkYoRjBGM0Y2RltzRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZcc0Zcb0Zpcy1JJW1zdXBHRiU2JUZcdEZndC8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRigvJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRmV1LyUpYmV2ZWxsZWRHRjhGYnFGZHE3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JIlRHRigsJComLCYiIiYhIiIqJiIiKCIiIkkiYUdGKEZpdkZpdkZpdiwmRml2RmZ2KiYiIiNGaXYpRmp2Rl13Rml2Rml2RmZ2RmZ2NyNGYnY=
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSNtaUdGJTY5USVkaXNjRigvJSdmYW1pbHlHUTBUaW1lc35OZXd+Um9tYW5GKC8lJXNpemVHUSMxMkYoLyUlYm9sZEdRJmZhbHNlRigvJSdpdGFsaWNHUSV0cnVlRigvJSp1bmRlcmxpbmVHRjgvJSpzdWJzY3JpcHRHRjgvJSxzdXBlcnNjcmlwdEdGOC8lK2ZvcmVncm91bmRHUSpbMCwwLDI1NV1GKC8lK2JhY2tncm91bmRHUS5bMjU1LDI1NSwyNTVdRigvJSdvcGFxdWVHRjgvJStleGVjdXRhYmxlR0Y4LyUpcmVhZG9ubHlHRjsvJSljb21wb3NlZEdGOC8lKmNvbnZlcnRlZEdGOC8lK2ltc2VsZWN0ZWRHRjgvJSxwbGFjZWhvbGRlckdGOC8lMGZvbnRfc3R5bGVfbmFtZUdRKjJEfk91dHB1dEYoLyUqbWF0aGNvbG9yR0ZELyUvbWF0aGJhY2tncm91bmRHRkcvJStmb250ZmFtaWx5R0YyLyUsbWF0aHZhcmlhbnRHUSdpdGFsaWNGKC8lKW1hdGhzaXplR0Y1LUkjbW9HRiU2M1EjOj1GKC8lJWZvcm1HUSZpbmZpeEYoLyUmZmVuY2VHRjgvJSpzZXBhcmF0b3JHRjgvJSdsc3BhY2VHUS90aGlja21hdGhzcGFjZUYoLyUncnNwYWNlR0ZbcC8lKXN0cmV0Y2h5R0Y4LyUqc3ltbWV0cmljR0Y4LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOC8lLm1vdmFibGVsaW1pdHNHRjgvJSdhY2NlbnRHRjgvJTBmb250X3N0eWxlX25hbWVHRlgvJSVzaXplR0Y1LyUrZm9yZWdyb3VuZEdGRC8lK2JhY2tncm91bmRHRkctSSZtZnJhY0dGJTYqLUYkNiMtRiQ2Ky1JI21uR0YlNjlGZ3BGMEYzRjYvRjpGOEY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ24vRmpuUSdub3JtYWxGKEZcby1GX282M1EoJm1pbnVzO0YoRmJvRmVvRmdvL0Zqb1EwbWVkaXVtbWF0aHNwYWNlRigvRl1wRmdyRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiUtRl5yNjlRIjJGKEYwRjNGNkZgckY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GYXJGXG8tRl9vNjNRMSZJbnZpc2libGVUaW1lcztGKEZib0Zlb0Znby9Gam9RJDBlbUYoL0ZdcEZic0ZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GLTY5USJhRihGMEYzRjZGOUY8Rj5GQEZCRkVGSEZKRkxGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8tRl9vNjNRIitGKEZib0Zlb0Znb0ZmckZockZecEZgcEZicEZlcEZocEZqcEZccUZecUZgcUZicUZkcS1GJDYlLUZecjY5USM0OUYoRjBGM0Y2RmByRjxGPkZARkJGRUZIRkpGTEZORlBGUkZURlZGWUZlbkZnbkZhckZcb0Zecy1JJW1zdXBHRiU2JUZkc0Zbcy8lMXN1cGVyc2NyaXB0c2hpZnRHUSIwRihGY3ItRiQ2JS1GXnI2OVEkMTM2RihGMEYzRjZGYHJGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmFyRlxvRl5zLUZgdDYlRmRzLUZecjY5USIzRihGMEYzRjZGYHJGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmFyRlxvRmJ0RmdzLUYkNiUtRl5yNjlRIzk2RihGMEYzRjZGYHJGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmFyRlxvRl5zLUZgdDYlRmRzLUZecjY5USI0RihGMEYzRjZGYHJGPEY+RkBGQkZFRkhGSkZMRk5GUEZSRlRGVkZZRmVuRmduRmFyRlxvRmJ0LUYkNiMtRmB0NiUtRiQ2JS1GX282M1EiKEYoL0Zjb1EncHJlZml4RigvRmZvRjtGZ28vRmpvUS50aGlubWF0aHNwYWNlRigvRl1wRmZ2L0ZfcEY7RmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxLUYkNiYtRl9vNjNRKiZ1bWludXMwO0YoRmJ2RmVvRmdvRmFzRmNzRl5wRmBwRmJwRmVwRmhwRmpwRlxxRl5xRmBxRmJxRmRxRl1yRmdzLUYkNiVGW3NGXnNGX3QtRl9vNjNRIilGKC9GY29RKHBvc3RmaXhGKEZkdkZnb0Zldi9GXXBRMnZlcnl0aGlubWF0aHNwYWNlRihGaHZGYHBGYnBGZXBGaHBGanBGXHFGXnFGYHFGYnFGZHFGW3NGYnQvJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRlx4LyUpYmV2ZWxsZWRHRjhGYnFGZHE3Iy1fRilJLG1wcmludHNsYXNoR0YoNiQ3Iz5JJWRpc2NHRigqJiwsIiIiRlt5KiYiIiNGW3lJImFHRihGW3khIiIqJiIjXEZbeSlGXnlGXXlGW3lGW3kqJiIkTyJGW3kpRl55IiIkRlt5Rl95KiYiIycqRlt5KUZeeSIiJUZbeUZbeUZbeSksJkZbeUZfeSomRl15Rlt5RmJ5Rlt5Rlt5Rl15Rl95NyNGaXg=
solve(DD=0);
solve(T=0);
NiQtSSVtcm93RzYjL0krbW9kdWxlbmFtZUc2IkksVHlwZXNldHRpbmdHSShfc3lzbGliR0YoNiUtSSZtZnJhY0dGJTYqLUkjbW5HRiU2OVEiMkYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR0Y7LyUqdW5kZXJsaW5lR0Y7LyUqc3Vic2NyaXB0R0Y7LyUsc3VwZXJzY3JpcHRHRjsvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y7LyUrZXhlY3V0YWJsZUdGOy8lKXJlYWRvbmx5R1EldHJ1ZUYoLyUpY29tcG9zZWRHRjsvJSpjb252ZXJ0ZWRHRjsvJStpbXNlbGVjdGVkR0Y7LyUscGxhY2Vob2xkZXJHRjsvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGRi8lL21hdGhiYWNrZ3JvdW5kR0ZJLyUrZm9udGZhbWlseUdGNS8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRigvJSltYXRoc2l6ZUdGOC1GMDY5USIzRihGM0Y2RjlGPEY+RkBGQkZERkdGSkZMRk5GUUZTRlVGV0ZZRmZuRmhuRmpuRlxvRl9vLyUubGluZXRoaWNrbmVzc0dRIjFGKC8lK2Rlbm9tYWxpZ25HUSdjZW50ZXJGKC8lKW51bWFsaWduR0Zpby8lKWJldmVsbGVkR0Y7LyUrZm9yZWdyb3VuZEdGRi8lK2JhY2tncm91bmRHRkktSSNtb0dGJTYzUSIsRigvJSVmb3JtR1EmaW5maXhGKC8lJmZlbmNlR0Y7LyUqc2VwYXJhdG9yR0ZQLyUnbHNwYWNlR1EkMGVtRigvJSdyc3BhY2VHUTN2ZXJ5dGhpY2ttYXRoc3BhY2VGKC8lKXN0cmV0Y2h5R0Y7LyUqc3ltbWV0cmljR0Y7LyUobWF4c2l6ZUdRKWluZmluaXR5RigvJShtaW5zaXplR1EiMUYoLyUobGFyZ2VvcEdGOy8lLm1vdmFibGVsaW1pdHNHRjsvJSdhY2NlbnRHRjsvJTBmb250X3N0eWxlX25hbWVHRmVuLyUlc2l6ZUdGOEZecEZgcC1GLTYqRmFvLUYwNjlRIjRGKEYzRjZGOUY8Rj5GQEZCRkRGR0ZKRkxGTkZRRlNGVUZXRllGZm5GaG5Gam5GXG9GX29GZG9GZ29Gam9GXHBGXnBGYHA3IzYkIyIiIyIiJCNGYHMiIiU=
NiQtSSZtZnJhY0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGKDYqLUkjbW5HRiU2OVEiNUYoLyUnZmFtaWx5R1EwVGltZXN+TmV3flJvbWFuRigvJSVzaXplR1EjMTJGKC8lJWJvbGRHUSZmYWxzZUYoLyUnaXRhbGljR0Y4LyUqdW5kZXJsaW5lR0Y4LyUqc3Vic2NyaXB0R0Y4LyUsc3VwZXJzY3JpcHRHRjgvJStmb3JlZ3JvdW5kR1EqWzAsMCwyNTVdRigvJStiYWNrZ3JvdW5kR1EuWzI1NSwyNTUsMjU1XUYoLyUnb3BhcXVlR0Y4LyUrZXhlY3V0YWJsZUdGOC8lKXJlYWRvbmx5R1EldHJ1ZUYoLyUpY29tcG9zZWRHRjgvJSpjb252ZXJ0ZWRHRjgvJStpbXNlbGVjdGVkR0Y4LyUscGxhY2Vob2xkZXJHRjgvJTBmb250X3N0eWxlX25hbWVHUSoyRH5PdXRwdXRGKC8lKm1hdGhjb2xvckdGQy8lL21hdGhiYWNrZ3JvdW5kR0ZGLyUrZm9udGZhbWlseUdGMi8lLG1hdGh2YXJpYW50R1Enbm9ybWFsRigvJSltYXRoc2l6ZUdGNS1GLTY5USI3RihGMEYzRjZGOUY7Rj1GP0ZBRkRGR0ZJRktGTkZQRlJGVEZWRllGZW5GZ25GaW5GXG8vJS5saW5ldGhpY2tuZXNzR1EiMUYoLyUrZGVub21hbGlnbkdRJ2NlbnRlckYoLyUpbnVtYWxpZ25HRmZvLyUpYmV2ZWxsZWRHRjgvJStmb3JlZ3JvdW5kR0ZDLyUrYmFja2dyb3VuZEdGRjcjIyIiJiIiKA==
plot(disc,a=0..2/3);
plot(disc,a=3/4..1);
NiUtSSdDVVJWRVNHNiI2JDdVNyQkIiIhRiokIiIiRio3JCQiMys2NjY2TzlgOSEjPiQiM15iPUgvMys8KSohIz03JCQiM0tBQUEoNDt2ciNGMCQiMyc+bC87LHYwIykqRjM3JCQiM0FXV1clPk4lUlRGMCQiM01SdWUkeTNrKSoqRjM3JCQiM01XV1clcHoyZCZGMCQiMyRlb1V0J3oiMy4iISM8NyQkIjNxKSkpKSlRRUBgKnBGMCQiMyVvPmZGeSdldzVGQzckJCIzSnh4eF9PMDskKUYwJCIzWWddJnAjXCcqSDZGQzckJCIzM21tbVRiZiRvKkYwJCIzOXVUejxsOyY+IkZDNyQkIjNNV1clPiQ0ejQ2RjMkIjNxKipwYG5DNnM3RkM3JCQiMyUpKioqKlxQKG8yRCJGMyQiMytzQVA5bFNkOEZDNyQkIjMhKSkpKSkpKVE8eSZSIkYzJCIzPW0xaidSOEpYIkZDNyQkIjNQV1dXcCM0Tl8iRjMkIjNZJzNNOyZHUlY6RkM3JCQiM3NtbW07LUluO0YzJCIzUUxAMW5nK147RkM3JCQiM0hMTEwkZSJvNj1GMyQiM10nUksleWN1azxGQzckJCIzQUxMTCR5PjMmPkYzJCIzIUhTUkFOVyJ6PUZDNyQkIjMzNjZoKWZxcjIjRjMkIjMhR2puPE5rbSk+RkM3JCQiMyopKSkpKSkpUVJURkFGMyQiM2lGO0chNEgmPUBGQzckJCIzV2JiYmIoKm9hQkYzJCIzKm9DcW92JkdMQUZDNyQkIjNsKioqKlwoeVRGXSNGMyQiM3EyREIqUTAqcEJGQzckJCIzc0BBQUFIekxFRjMkIjNDOV8rOjxOJFwjRkM3JCQiMzkrK10oeXd2eCNGMyQiMzQvc0pGIj43aiNGQzckJCIzYG1tO0hJXDlIRjMkIjN3TStmPFtpa0ZGQzckJCIzOyopKSkpUXdddDAkRjMkIjMuJHknZiQ9NmUhSEZDNyQkIjMrQUFzZihRJik9JEYzJCIzZiQ9bmMpejRQSUZDNyQkIjNvNTY2T0YvSUxGMyQiMydcPixaIlFLIT0kRkM3JCQiMyU0NjYnKWZEcVokRjMkIjMyMU4iSDphMkwkRkM3JCQiMyUpKioqKlw3WShcZyRGMyQiM2ddX1o4YCxqTUZDNyQkIjNoeHh4X0k7VlBGMyQiM3BFbiN6Tig+Mk9GQzckJCIzdUtMTExfI2YpUUYzJCIzYS51KUhFNXd2JEZDNyQkIjNcKioqKipcPCFmRFNGMyQiMyVIJXBxWHpCMVJGQzckJCIzWyoqKipcKD1CMjslRjMkIjMlZm5SM3E+OjAlRkM3JCQiM3VLTExlYHc1VkYzJCIzIyplaz5UWHA5VUZDNyQkIjNlQEFBQVxlWFdGMyQiM0NePioqUnA+alZGQzckJCIzUUxMTCRlTCYqZSVGMyQiMyllZl9NJDQuQ1hGQzckJCIzW3h4RmxDKCo+WkYzJCIzUyZ5JFFoUTpzWUZDNyQkIjN1bW1tO21kaVtGMyQiMydHJlFFOUw/UFtGQzckJCIzTUBBczRndicqXEYzJCIzWyZILUJaM2gqXEZDNyQkIjMnKioqKipcUEAsUF5GMyQiM10oRzxNX1dvOyZGQzckJCIzZykpKSkpKVEiUVRGJkYzJCIzW1BgI0cqM19SYEZDNyQkIjNpbG07YWpwPGFGMyQiM0d1RCpbIls2R2JGQzckJCIzJXp4eHgoKmZmYiZGMyQiMyFSX08qKXk+Jz5kRkM3JCQiM0F4eHgtWE4ocCZGMyQiMyk0ZTQ+M1YpR2ZGQzckJCIzdkBBczQjeXYkZUYzJCIzUUhdYE5PaGFoRkM3JCQiM19MTExMI0hrJ2ZGMyQiMyF6MW5OOS1YUSdGQzckJCIzYUFBQVpxNTloRjMkIjM/SVVqUVVoJ28nRkM3JCQiMz9XV1dXPT5ZaUYzJCIzODwrLmF0MTVxRkM3JCQiMzsrK10oeUBxUSdGMyQiMydbelFdUiM9WXVGQzckJCIzOysrKyt6VGFrRjMkIjMmMz0iZVY3XzJ4RkM3JCQiMzkrK103UyI9XydGMyQiM1VrOzozUyozLSlGQzckJCIzbUtMZVIuQyVmJ0YzJCIzVj5KNEJ5QlUlKUZDNyQkIjNJbW1tbW1tbW1GMyQiM1lPMCsrKysrISpGQy1JJ0NPTE9VUkdGJTYmSSRSR0JHRiUkIiM1ISIiRilGKS1JK0FYRVNMQUJFTFNHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiU2JFEiYUYlUSFGJS1JJVZJRVdHRl1cbDYkO0YpJCIrbm1tbW0hIzVJKERFRkFVTFRHRiU=
NiUtSSdDVVJWRVNHNiI2JDdTNyQkIjMrKysrKysrK3YhIz0kIiIlIiIhNyQkIjMzbm1UJilHXGF2RiskIjM2ejgmUSgpSGtKJSEjPDckJCIzVExla2BvIT5nKEYrJCIzI1x6czpYMydbWEY0NyQkIjNDbTt6PilHX2woRiskIjMpcC1KJFwmZlB4JUY0NyQkIjM9bTthUVUhKjN4RiskIjMnKSpvRnFYbjIoXEY0NyQkIjNhTGVSWlhLaXhGKyQiM3AneWhebERROSZGNDckJCIzem0iej4sXz0ieUYrJCIzMyxzITRtIip5RyZGNDckJCIzWypcN0ckWzhqeUYrJCIzTTsmbygpSFFOVSZGNDckJCIzL24ieiUqZnJoInpGKyQiM1VWSjAlXGs9YiZGNDckJCIzMytEY3cjUSFwekYrJCIzPy1kJCpRJCpwcGNGNDckJCIza0xMM18iPU0tKUYrJCIzSTY9R3BQNCN5JkY0NyQkIjNTbTsvd2ZKciEpRiskIjM7SEdtNkR5dWVGNDckJCIzJioqKipcN2VQXzcpRiskIjMncCpbI1xeYkkoZkY0NyQkIjNzKioqXFBmIVF6IilGKyQiMzdxL2RuZTdtZ0Y0NyQkIjNlKioqXCg9dWJKIylGKyQiMzs8UClwdic+XmhGNDckJCIzUW0ielcoKlEqeSMpRiskIjN2K25BVUwqXEEnRjQ3JCQiM3NMTDNGLUdOJClGKyQiM1QkW0MlejAnKjNqRjQ3JCQiMylITEwkZSczSVEpRiskIjNJd1ZDMVJBeGpGNDckJCIzTitESnEiRyZRJSlGKyQiMz9XW0Z2Om9ga0Y0NyQkIjNpTEwkZU1zd1spRiskIjM0OSRcKm8jeSo9bEY0NyQkIjNoK0RKJkgiZlQmKUYrJCIzaUxdLCM0LiQpZSdGNDckJCIzJykqXFBmKVskSGYpRiskIjNBVV1VQlhGX21GNDckJCIzeUwkZWtgMWxrKUYrJCIzLVplYVNTNDxuRjQ3JCQiM2VMZSpbLi1kcClGKyQiM01VYDxJUS52bkY0NyQkIjNrbTsvRWd3WygpRiskIjM1Tj4oUU5yZiRvRjQ3JCQiMyJmO3olKmYlKVEhKSlGKyQiM19MLV0mPSZwKCpvRjQ3JCQiM2cqXChvemEnPSYpKUYrJCIzRUxfVUFFQV1wRjQ3JCQiMzRuO3pXaG8uKilGKyQiM0p4SSFSUip5MHFGNDckJCIzTyoqKipcaT5BZCopRiskIjNtcEMjPng6PzEoRjQ3JCQiMz4qKlxpOmpmNCEqRiskIjMqUS1yczVXZjYoRjQ3JCQiM0UqXDdgPnItMSpGKyQiMy1VTTgxaDxuckY0NyQkIjN5KipcUDRxYDsiKkYrJCIzJSp5VC4qPVVJQShGNDckJCIzO0xMJGVNJTRuIipGKyQiM3E+K0ViVFFzc0Y0NyQkIjNuKioqXFA0djVBKkYrJCIzWiJwSyxOR1VLKEY0NyQkIjNjbSJ6V24qKSpwIypGKyQiMy9mIUdacCwwUChGNDckJCIzVCoqKlw3Qm1NSypGKyQiM0RYT1QoKW9PP3VGNDckJCIzIUckZWsuTnl0JCpGKyQiMyFlbzFTbUttWShGNDckJCIzUytEY14mempVKkYrJCIzKHo5QlhFX1ZeKEY0NyQkIjNDTEwzLT0heVoqRiskIjNbS1A3LF1TZ3ZGNDckJCIzKikqXDdHOE87YCpGKyQiM2RbTDZEQCwzd0Y0NyQkIjNXbW1tIipcWyRlKkYrJCIzbUcrOkZkSWB3RjQ3JCQiM0ttO2FRel1PJypGKyQiMydbaVxBYiYzKnAoRjQ3JCQiM3VLZWtHPTQqbypGKyQiM3kmSHYsJD0oUnUoRjQ3JCQiM1EqKioqXGk0VFAoKkYrJCIzJSpHLHNCX3kleShGNDckJCIzL0wkM0Y5IXojeipGKyQiMykqPihwQSVSMkp5RjQ3JCQiM0ltbW1UPktVKSpGKyQiM1UxQSp5MV4/KHlGNDckJCIzMStESnFKOCYqKSpGKyQiM0hNOioqUXdKOnpGNDckJCIzMSt2b2Etb1gqKkYrJCIzcSJ5YShbXUxjekY0NyQkIiIiRi4kIiIpRi4tSSdDT0xPVVJHRiU2JkkkUkdCR0YlJCIjNSEiIiRGLkYuRmFbbC1JK0FYRVNMQUJFTFNHNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiU2JFEiYUYlUSFGJS1JJVZJRVdHRmRbbDYkOyQiKysrKyt2ISM1RmZ6SShERUZBVUxUR0Yl
plot(T,a=0..2/3);
plot(T,a=3/4..1);
subs(a=2/3,eq);
subs(a=3/4,eq);