A method is presented for calculating the shape sensitivity of a wing aeroelastic response with respect to changes in geometric shape. Yates’ modified strip method is used in conjunction with Giles' equivalent plate analysis to predict the flutter speed, frequency, and reduced frequency of the wing. Three methods are used to calculate the sensitivity of the eigenvalue. The first method is purely a finite difference calculation of the eigenvalue derivative directly from the solution of the flutter problem corresponding to the two different values of the shape parameters. The second method uses an analytic expression for the eigenvalue sensitivities of a general complex matrix, where the derivatives of the aerodynamic, mass, and stiffness matrices are computed using a finite difference approximation. The third method also uses an analytic expression for the eigenvalue sensitivities but the aerodynamic matrix is computed analytically. All three methods are found to be in good agreement with each other. The sensitivities of the eigenvalues were used to predict flutter speed, frequency , and reduced frequency. These approximations were found to be in good agreement with those obtained using a complete reanalysis. However, it is recommended that higher order terms be used in the calculations in order to assure greater accuracy. / Master of Science / incomplete_metadata