n=100;
alpha=0.1;


[Cp_k, Cl_k] = ktreff(n, 0, 0.08, 0.066, 0.15, "ktreff.xy"); 
points = load("ktreff.xy"); 
[cp, cl] = dpan(n, alpha, points(:,1), points(:,2)); 
hold on
plot(points(:, 1), points(:, 2)); 
plot(points(:,1), -(cp(:)));
%plot(temp(:,1), -(Cp_k(:)));

legend("airfoil", "dpan", "ktreff")

axis equal
cl
hold off