x_teorico = [0:.01:0.6];
L_teorico = 100*(x_teorico/0.6).^(0.625).*((1-x_teorico)/0.4).^(-1.625);
x_exp = xlsread('dati_toluene_II','a2:a31');
L_exp_I = xlsread('dati_toluene_II','b2:b31');
L_exp_II = xlsread('dati_toluene_II','c2:c31');
plot(L_teorico,x_teorico),grid,xlabel('L (mol)'),ylabel('x (mol B/mol)');
hold on;
plot(L_exp_I,x_exp,'r*'),grid,xlabel('L (mol)'),ylabel('x (mol B/mol)');
hold on;
plot(L_exp_II,x_exp,'g*'),grid,xlabel('L (mol)'),ylabel('x (mol B/mol)');
legend('curva teorica','I esperimento','II esperimento')