foil = load("ktreff.xy"); 
alpha_low = -10;
alpha_high = 20;

cls = zeros(alpha_high - alpha_low, 1);

for a = alpha_low:alpha_high
    [cp, cl] = dpan(100, a*pi/180, foil(:,1), foil(:,2));
    cls(1+a-alpha_low) = cl;
end