function [Cl] = Cl(alpha)

alpha_data = 0:pi/20:pi/2; 
CL_data = [0 1 2 3 2 1 2 3 1 0 ];  
Cl = interp1(alpha_data,CL_data,alpha,'spline'); 
plot(x,y,'o',xi,yi)