// e be the electron charge in coulombs (c)
e=-1.6*10^-19;
// n be the no. of electrons
n=4600;
//total charge q in coulombs (c)
q=n*e;
disp("q=")
disp(q)
units='coulombs C';
q1=[string(q) units];
disp(q1)
// In coulombs
// the total charge is -7.36*10^-16 coulombs