cl=linspace(-5,5,1000);
cd=0.021+0.0537*cl.^2;

clsper=[-0.203 -0.019 0.166 0.352 0.536 0.722 0.907 1.093 1.278 1.382 1.412 1.385];
cdsper=[0.023 0.021 0.024 0.029 0.037 0.048 0.064 0.083 0.108 0.138 0.176 0.182];


plot(cd,cl,'b',cdsper,clsper,'r');
title('polare velivolo completo')
grid on 
xlabel('Cd')
ylabel('Cl')
legend('analitica','sperimentale')
axis([0 0.3 0 2.5])