# Solve for the equil points
# x'=60x-3x^2-4xy, y'=42y-3y^2-2xyF1:=(x,y)->60*x-3*x^2-4*x*y; F2:=(x,y)->42*y-3*y^2-2*x*y; solve({F1(x,y)=0,F2(x,y)=0},{x,y});F:=unapply(<F1(x,y),F2(x,y)>,(x,y)):F(x,y);J:=unapply(<diff(F1(x,y),x),diff(F2(x,y),x)|diff(F1(x,y),y),diff(F2(x,y),y)>,(x,y)):J(x,y);with(LinearAlgebra):X:=[0,0,20,12]:Y:=[0,14,0,6]:for i from 1 to 4 do 'J'(X[i],Y[i]) =J(X[i],Y[i]); od;TTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODgzODY5MzgyNzI0NTM0WColKWFueXRoaW5nRzYiRiVbZ2whIyUhISEiIyIjeik2JCUieEclInlHW1tbW1tfW1tbW3MhIiRcW1tbW11bW1tbcyEiJVtbW1tbXVtbW1tvIiNneilGJ1xbW1tbXVtbW1tzISIjXVtbW1tbW1tbW3NGKlxbW1tbW1tbW1tvIiNVRiU=TTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODgzODY5MzgyNzI1MDE0WCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjJSIjIiMsKCUieEchIiclInlHISIlIiNnIiIiLCRGKSEiIywkRidGKiwoRidGLkYpRigiI1VGLEYlTTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODgzODY5Mzg0NjMzMDMwWCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjJSIjIiMiI2ciIiFGJyIjVUYlTTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODgzODY5Mzg0NjMzMTUwWCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjJSIjIiMiIiUhI0ciIiEhI1VGJQ==TTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODgzODY5Mzg0NjMzMjcwWCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjJSIjIiMhI2ciIiEhIyEpIiIjRiU=TTdSMApJQVJUQUJMRV9TQVZFLzE4NDQ2ODgzODY5Mzg0NjMzMzkwWCwlKWFueXRoaW5nRzYiRiVbZ2whIiUhISEjJSIjIiMhI08hIzchI1shIz1GJQ==