function xdot=P6_1fun(t,x);

gm=.05;cm=1;gc=.1;I1=0;I2=0;

A=[-(gm+gc)/cm gc/cm;
    gc/cm  -(gm+gc)/cm];
F=[I1/cm;I2/cm];

xdot=A*x+F;