function q = prob3fun(x)
rho = 8960;
c = 389;
A = (-19.5*x.^2 + 15.5*x + 5)/1000;
T = 100(1-sinh(x));
q = rho.*c.*A.*T;