Allora, il codice che ho in questo momento è il seguente:
nu=0.25;
phi=0.10;
varepsilon=linspace(0,1.5,45);
if varepsilon<0.6
x=((1+nu)/2)+((sqrt(((1+nu^2)/2)-((1-nu)/2)))*(varepsilon/0.6))
else x=(sqrt((1+nu^2)/2))+(0.5*((sqrt(((1+nu^2)/2)-((1-nu)/2))))*((varepsilon-0.6)/0.6));
end
cl=1;
cd=0.012;
cds=(0.018*(cl^2));
sigh=1.2;
ar=1.7;
cda=((0.02)/(sigh*nu*ar));
cdtip=(0.029*0.015*(cl^(1.5)));
tanbeta1=1/(phi/(1-nu^2));
tanbeta2=(1-(varepsilon.*(phi/(1-nu^2))))/(phi/(1-nu^2));
tanbetam=0.5*(tanbeta1+tanbeta2);
a=(((8/pi)*(phi./((1-nu^2).*x))).^-1);
b=((8/pi)*((phi./((1-nu^2).*x))).*(1+((varepsilon.^2)/4)));
c=varepsilon;
d=(cl/(cd+cds+cds+cdtip));
e=tanbetam;
f=((2/pi)*(phi/((1-nu^2))));
g=(2./x);
h=((nu^2)./varepsilon)*(2-nu^2)*(0.4+(nu/2));
etaaer=1-((a+b-c)/(d-e))-(f*g.*(varepsilon+h));
figure
plot(varepsilon,etaaer,'b')
hold on
grid on
plot(varepsilon,etaaer, 'black', 'LineWidth', 2)
axis([0 1.5 0 1]);
set(gca, 'XTick',0 : 0.25 : 1.5);
xlabel('\epsilon_s')
ylabel('\eta_{aer}')
title('Grafico')
Non saprei come fare per il calcolo delle curve ad altri valori di nu nello stesso grafico.