function xarea = tank_cross_section(h)
% Returns the cross-sectional area of a tank
xarea = pi*(h.^2+(1/3)*h.^sqrt(3)).^2;