A numerical procedure, which utilises polynomial powerseries expansions for the optimisation of multipanelwing structures in idealised critical flutterconditions, is introduced and developed. It arisesfrom the Rayleigh-Ritz method and employes trialpolynomial describing functions both for theflexural displacement and for the thicknessvariation over the multipanel surface. An idealisedstructural plate model, according to the Kirchhoff’stheory, together with a linearised supersonicaerodynamic approach, are supposed. The classicalEuler-Lagrange optimality criterion, based onvariational principles, has been utilised for theoptimisation operations, where by imposing thestationary conditions of the Lagrangian functionalexpression, a nonlinear algebraic equations systemis obtained, whose solution is found by anappropriate algorithm. By utilising an iterativeprocess it is possible to reach the referencestructure critical conditions, with an optimisedthickness distribution throughout the multipanelsurface. The final part of the work consists insearching the minimum weight of the multipanelplanform wing structure with optimised thicknessprofile vs the flutter frequency,considered as a variable imput parameter, for fixedflutter speed and equal to the critical one of thereference uniform structure.