% local blade detailsprompt = ' please enter the details of the blade radius in metres : ';radius_r= input(prompt);prompt1 = ' Please enter the chord lenght for your blade element in metres: ';local_chord_C= input(prompt1);prompt3= ' Please enter the wind speed in metres/second : ';wind_speed= input(prompt3); % constant variables for the bladegamma= 92.6;frequency= 2;omega= 2*3.14*frequency;Radius_R=0.22;Area_A= 3.14*Radius_R*Radius_R;density_air= 1.225;efficiency_gearbox= 0.9;efficiecny_generator= 0.8;number_blades=3;F_tip=1; % formulastip_speed_ratio= Radius_R *omega/wind_speed;local_tip_speed_ratio= tip_speed_ratio*radius_r / Radius_R;local_solidity= number_blades*local_chord_C/2*3.14 * radius_r;Si_radians= 2/3*atan(1/local_tip_speed_ratio); % Si in radiansSi_degrees = Si_radians * 57.3;alpha= gamma + Si_degrees -90;coefficient_Lift = 2* 3.14 * sin(alpha);angular_induction_factor = 1/(1 + (4*F_tip* sin(alpha)*sin(alpha))/ (local_solidity* coefficient_Lift*cos(alpha)));axial_induction_factor = (1- 3*angular_induction_factor)/(4*angular_induction_factor-1);Coefficience_performance = (8 * F_tip * local_tip_speed_ratio*local_tip_speed_ratio*local_tip_speed_ratio*axial_induction_factor*(1-axial_induction_factor))/(tip_speed_ratio * number_blades);Work_extracted = (Coefficience_performance * efficiency_gearbox * efficiecny_generator * density_air * Area_A * wind_speed *wind_speed*wind_speed*0.5);fprintf(' %i the work extracted from the turbine is \n',Work_extracted);