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

v_w=40;     %windspeed in km/h
WD=-180;    %winddirection in °
TC=0;       %true course in °

for v_w=1:100
    [v_best,res_eff_best(v_w)] = eff_calc(interp_power,v_a, TC, v_w, WD );
end
legend('v_w von 1 bis 100 km/h' )


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