function Vdot = twocan(t,V)

Kideal1 = .0015579419;
Kideal2 = .0005193107896;

V1 = V(1);
V2 = V(2);

if (V(1)<=0)
    V1=0;
end
    Vdot1 = -Kideal1*sqrt(V1);

if (V(2)<=0)
    V2=0;
end
    Vdot2 = Kideal1*sqrt(V1)-Kideal2*sqrt(V2);

Vdot = [Vdot1;Vdot2];