%% Function (xsystem_case2.m) %%
function xdot = xsystem_case2(t,x)
global Q L V1o V2o V3o K1 K2 K3

V1 = V1o + K1*(x(2)); % voltage on sphere 1 [V]

V2 = V2o + K2*(x(2)); % voltage on sphere 2 [V]

V3 = V3o + K3*(x(2)); % voltage on sphere 3 [V]

if V1 > 30000
    V1 = 30000;
elseif V1 < -30000
    V1 = -30000;
end

if V3 > 30000
    V3 = 30000;
elseif V3 < -30000
    V3 = -30000;
end

%% Differential Equations %%

xdot(1) = x(2);
xdot(2) = V2*Q*((V1./(x(1)^2)) - (V3./(L -x(1))^2));

xdot = xdot';