// q be the charge in coulombs and i be the current in Ampere
 
// i is given by i=(3*t^2-t)
 
//charge is to be found between t=1s and t=2s
 
integrate('3*t^2-t','t',1,2);
 
// the charge between 1s and 2s is 5.5 coulombs
 
q=ans;
 
disp("q=")
disp(q)
units='Coulombs C'
q=[string(q) units];
disp(q)
// in coulombs


// the charge is 5.5 coulombs