Analytical indicial aerodynamic functions are calculated for several trapezoidal wings in subsonic flow, with a Mach number 0.3 ≤ Ma ≤ 0.7. The formulation herein proposed extends well-known aerodynamic theories, which are limited to thin aerofoils in incompressible flow, to generic trapezoidal wing planforms. Firstly, a thorough study is executed to assess the accuracy and limitation of analytical predictions, using unsteady results from two state-of-the-art computational fluid dynamics solvers as cross-validated benchmarks. Indicial functions are calculated for a step change in the angle of attack and for a sharp-edge gust, each for four wing configurations and three Mach numbers. Then, analytical and computational indicial responses are used to predict dynamic derivatives and the maximum lift coefficient following an encounter with a one-minus-cosine gust. It is found that the analytical results are in excellent agreement with the computational results for all test cases. In particular, the deviation of the analytical results from the computational results is within the scatter or uncertainty in the data arising from using two computational fluid dynamics solvers. This indicates the usefulness of the developed analytical theories. Keywords: Analytical approach, CFD, Compressible flow, Gust response, Indicial aerodynamics, Trapezoidal wing