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

/* [wxMaxima: title   start ]
state equations
   [wxMaxima: title   end   ] */

/* [wxMaxima: comment start ]

1st step: form the equations, full system, 
reduced to algabraic (over Dx, y, x, u)

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
e1: i1+i2+i3-i4=0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e2: i4-i5=0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e3: i5+i6=0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e4: u1=v1;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e5: u2=v1;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e6: u3=v1;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e7: u4=v2-v1;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e8: u5=v3-v2;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e9: u6=v3;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e10: i1=C1*Du1;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e11: u2=R1*i2;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e12: i3=0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e13: u4=R2*i4;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e14: i5=C2*Du5;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e15: u6=A*u3;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
S: [e1, e2, e3, e4, e5, e6, e7, e8, e9, e10, e11, e12, e13, e14, e15];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]

2nd step: eliminate y vector

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
y: [i1, i2, i3, i4, i5, i6, u2, u3, u4, u6, v1, v2, v3];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
SE: eliminate(S, y);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]

3rd step: solve over Dx

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
SE: solve(SE, [Du1, Du5]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
SE1: ratsimp(SE[1][1], u1, u5);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
SE2: ratsimp(SE[1][2], u1, u5);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]

4th step: exclude memory elements 
from the equation system; solve over y

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
O: [e1, e2, e3, e4, e5, e6, e7, e8, e9, e11, e12, e13, e15];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
OE: solve(O, y);
/* [wxMaxima: input   end   ] */

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