function [d,v,ch]=charge(t,x)

global c Q num_ch m gap volt

q0 = -Q*num_ch/(2*pi);
a0 = num_ch*Q*volt/(2*m*gap);
d = gap+a0*t.^2;
v = 2*a0*t;
ch = q0*d./(d.^2.+x.^2).*(d./(d.^2.+x.^2).^(1/2)+v./c);
