clear all, close all
[v_a,interp_power] = load_aircraft_eff()

v_w=40  ;       %windspeed in km/h
WD=60    ;       %winddrection in °


for TC=-180:180
    [v_best,res_eff_best(TC+181)] = eff_calc(interp_power,v_a, TC, v_w, WD );
end
legend('TC von -180° bis 180°' )

figure(2)
plot(-180:180,res_eff_best, 'linewidth', 3);
xlabel('windspeed/ km/h')
ylabel('efficiency / %')