here is the matlab coding: % First, define symbolic variable syms h % Define bounds of integration (height of the structure) a = 0; b = 20; % Define the integrand components rho = 240000 / (log(21) * (h + 1)); % density A = 4 * (400 - h^2); % cross-sectional area g = 9.8; % gravitational constant integrand = rho * A * g * h; % total integrand for energy % Try symbolic antiderivative E_anti = int(integrand); % indefinite integral disp('Indefinite integral (antiderivative):') pretty(E_anti) % Compute definite integral symbolically E_def = int(integrand, h, a, b); disp('Exact total lifting energy (symbolic):') disp(E_def) % Approximate result numerically E_approx = vpa(E_def); disp('Numerical approximation of lifting energy:') disp(E_approx)