In order to improve the effectiveness, increase the life cycle, and avoid the blade structural failure of wind turbines, the blades need to be perfectly designed. Knowing the flow angle and the geometric characteristics of the blade is necessary to calculate the values of the induction factors (axial and tangential), which are the basis of the Blade Element Momentum theory (BEM). The aforementioned equations form an implicit and nonlinear system. Consequently, a straightforward iterative solution process can be used to solve this problem. A theoretical study of the aerodynamic performance of a horizontal-axis wind turbine blade was introduced using the BEM. The main objective of the current work is to examine the wind turbine blade’s performance under specific initial and boundary conditions. In this study, NACA 4415 airfoils were selected to do this investigation. The turbine blades were divided into small segments to calculate the forces acting upon each segment to assess their impact on the final wind turbine blade design. An extensive and critical analysis of the chosen wind turbine was performed, including lift, drag, shear force, and bending moment calculations. Based on the computations, the values of total thrust force, torque, and power generation values for the optimal wind turbine were 3755.9 N, 1834 Nm, and 30.122 kW, respectively. Furthermore, a significant finding emerges from the analysis, indicating that the largest difference in power occurs at r/R=0.8, amounting to 5.5239 kW. The new Matlab code was validated. The key contribution of this study lies in enhancing turbine efficiency and reducing fatigue losses through optimizing wind turbine blade design to obtain the highest efficiency level. The outcomes demonstrate that the newly implemented MATLAB code exhibits exceptional accuracy in assessing aerodynamic performance, enabling efficient wind turbine blade design optimization.