This thesis presents the development of AVCOPro (Aerial Vehicle Construction and Optimization Program), a module created in the Python programming language, that is currently only applicable for optimizing the design of aircraft with blended fuselage and wing configurations. An optimization problem is formulated that includes designing the shape of such aircraft while verifying longitudinal and lateral dynamic stability. Previously developed mathematical models are described for defining the geometry of aircraft with blended wings and fuselage, calculating inertial properties, determining aerodynamic characteristics and verifying dynamic characteristics of the vehicle. These models, implemented in Python modules, are compared with available results on several representative models to verify their accuracy. A Python module has been developed that integrates the previously developed modules with existing implementations of applicable optimization algorithms, using the PyTorch library for automatic differentiation via backpropagation. This enables efficient calculation of gradients needed for optimization using the Adam optimizer. The developed module was applied to optimize the shape of an unmanned aircraft for the delivery of organs for transplantation. Finally, the limitations of the current version of the program are highlighted, and an overview of guidelines for future development is provided.