alpha = 3*pi/180;
eps = 0.08;
beta = 0.06;
tau = 0.15;
filename="temp.xy";
cls = zeros(10, 1);
for n = 1:140
    [~, cl] = ktreff(n, alpha, eps, beta, tau, filename);
    points=load(filename);
    [cp, cl2] = dpan(n, alpha, points(:,1), points(:,2)); 
    cls = [cls; cl2];
end
hold on
plot(cls)
plot([0, 140], [cl, cl])
hold off
xlabel("number of panels")
ylabel("coefficient of lift")
legend("dpan.m", "transformed flow")
%figure
%plot(log(1./(1:360)),log(abs(cls - cl)))