% Simulation of a circuit with two capacitors % Feb 2016, J. Gaspar R= 1e3; C= 1e-3; A= [-1 1; 1 -2]/(R*C); B= [1/C 0]'; C= [0 1]; % note this overwrites C D= 0; step(ss(A,B,C,D))