function qdot = twotankelectrical(t,q)

R1 = 215;
R2 = 1000;
R3 = 9800;

C1 = 220*10^-6;
C2 = 100*10^-6;

Vin = 8;

qdot1 = (Vin-(q(1)/C1))/R1 - (q(1)/C1-q(2)/C2)/R2;
qdot2 = (q(1)/C1-q(2)/C2)/R2 - (q(2)/C2)/R3;

qdot = [qdot1;qdot2];