function xdot=P6_2fun(t,x,flag,alpha,I1,I2);

gm=.05;cm=1;gc=.1;

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

xdot=A*x+F;