Extracted from http://prac.im.pwr.edu.pl/~bdyda/sem2019letni/autonomous-ode.wxm

autonomous-ode.wxm

/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 16.04.2 ] */

/* [wxMaxima: input   start ] */
plotdf([y,-sin(x) - b*y], [x,y],
        [parameters,"b=0.05"],
        [tstep,0.01],
        [x,-4,12], [y,-4,4], [direction,forward],
        [nsteps,1500], [sliders,"b=0:0.05"])$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
A:matrix([1,2],[3,2]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
eigenvalues(A);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
eivects(A);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
plotdf(transpose(A.[x,y])[1],[tstep,0.01],[nsteps,2000],[direction,forward]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
B:matrix([-1,2],[0,-3]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
eivects(B);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
plotdf(transpose(B.[x,y])[1],[tstep,0.01],[nsteps,2000],[direction,forward]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
C:matrix([0,1],[-1,0]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
eivects(C);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
plotdf(transpose(C.[x,y])[1],[tstep,0.01],[nsteps,2000],[direction,forward]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
D:matrix([5,1],[1,6]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
eivects(D);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
plotdf(transpose(D.[x,y])[1],[tstep,0.01],[nsteps,2000],[direction,forward]);
/* [wxMaxima: input   end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$