'1-SI 2-NO'
Q1=input('calibrazione già effettuata?: ')

if Q1==2

file_in1=input('Nome file calibrazione 1?   ','s');
file_in2=input('Nome file calibrazione 2?   ','s');
file_in3=input('Nome file calibrazione 3?   ','s');

file_in4=input('Nome file misura ?   ','s');


% PRIMO LIQUIDO salto il testo iniziale
 n=fopen(file_in1,'r');
 app=fgets(n);
 %i=1;
 % while  i<6
 %     app=fgets(n);
 %     i=i+1;
 % end
 C1=fscanf(n,'%f',[4 401]);
 C1=C1';
 
% SECONDO LIQUIDO salto il testo iniziale
 n1=fopen(file_in2,'r');
 app=fgets(n1);
 %i=1;
 % while  i<6
 %  app=fgets(n1);
 %  i=i+1;
 % end
 C2=fscanf(n1,'%f',[4 401]);
 C2=C2';

 % TERZO LIQUIDO salto il testo iniziale
 n2=fopen(file_in3,'r');
 app=fgets(n2);
 %i=1;
 %while  i<6
 %     app=fgets(n2);
 %     i=i+1;
 %end
 C3=fscanf(n2,'%f',[4 401]);
 C3=C3';



f1=C1(:,1);
Z0=50;
w1=2*pi*f1;

% LIQUIDO SOTTO MISURA 
 n3=fopen(file_in4,'r');
 app=fgets(n3);
 %i=1;
 %while  i<6
 %     app=fgets(n3);
 %     i=i+1;
 %end
 M=fscanf(n3,'%f',[4 401]);
 M=M';

f1=C1(:,1);
Z0=50;
w1=2*pi*f1;


e1=1;  % il primo liquido di calibrazione è l'aria

'1-acetone 2-etanolo 3-isopropil 4-Toluene 5-Iodiometano 6-cloroformio 7-metanolo 8-acqua distillata  '
M1=input('secondo liquido di calibrazione: ')

T1=input('temperatura liquido (°C): ')

e2=autoE(M1,T1,f1);

'1-acetone 2-etanolo 3-isopropil 4-Toluene 5-Iodiometano 6-cloroformio 7-metanolo 8-acqua distillata '

M2=input('terzo liquido di calibrazione: ')
T2=input('temperatura liquido (°C): ')

e3=autoE(M2,T2,f1);


else
file_in4=input('Nome file misura ?   ','s');

% LIQUIDO SOTTO MISURA 
 n3=fopen(file_in4,'r');
 app=fgets(n3);
% i=1;
% while  i<6
%      app=fgets(n3);
%      i=i+1;
% end
 M=fscanf(n3,'%f',[4 401]);
 M=M';
 
end

'1-acetone 2-etanolo 3-isopropil 4-Toluene 5-Iodiometano 6-cloroformio 7-metanolo 8-acqua 9-aria '
D1=input('liquido misurato è?: ')
if D1==9
    et=1;
else
    
T3=input('temperatura liquido (°C): ')
et=autoE(D1,T3,f1)
end

 
Epst1=real(et)
Epst2=-imag(et)

p1=complex(C1(:,3),C1(:,4));
p2=complex(C2(:,3),C2(:,4));
p3=complex(C3(:,3),C3(:,4));

pm=complex(M(:,3),M(:,4))

aln=(pm-p1).*(p2-p3).*(e3-e1)
ald=(pm-p2).*(p3-p1).*(e2-e3)
al=aln./ald

Eps=(e1-e2.*al)./(1-al)

Eps1=real(Eps)

Eps2=-imag(Eps)
figure;
plot(f1,Eps1)
hold on
plot(f1,Epst1,'r')

xlabel('f Hz')
ylabel('Eps1')




figure;

plot(f1,Eps2)
hold on

plot(f1,Epst2,'r')
xlabel('f Hz')
ylabel('Eps2')


%%%%%%%%%%%%%%%%%%%fitting%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%





% epssmin=input('epss min?   ');
% epssmax=input('epss max?   ');
% epsinfmin=input('epsinf min?   ');
% epsinfmax=input('epsinf max?   ');
% frmin=input('fr min (Hz)?   ');
% frmax=input('fr max (Hz)?   ');
% sigmasmin=0;
% sigmasmax=1e-9;
% 
% 
% [N_freq dummy]=size(C1)
% 
% freq=f1;
% epsrc_r=Eps1;
% epsrc_i=Eps2;
% epsrc_c=epsrc_r-1i*epsrc_i;
% 
% foptions = fitoptions('Method','NonLinearLeastSquares','Lower',[epssmin epsinfmin frmin 0 sigmasmin],'Upper',[epssmax epsinfmax frmax 1 sigmasmax],'Algorithm','Trust-Region','MaxFunEvals',1e4);
% type = fittype('abs(epsinf+(epss-epsinf)/(1+(1i*freq/fr)^(1-beta))+sigmas/(1i*2*pi*freq*8.8541878176e-12)-epsrc_c)','indep','freq','problem','epsrc_c','coeff',{'epss','epsinf','fr','beta','sigmas'},'options',foptions);
% fitted_abs = fit(freq,freq*0,type,'problem',epsrc_c);
% param_abs=coeffvalues(fitted_abs);
% disp(param_abs(1));
% disp(param_abs(2));
% disp(param_abs(3));
% disp(param_abs(4));
% disp(param_abs(5));
% 
% 
% epsrc_CC=f1;
% for n=1:1:N_freq
%     [epsrc_CC(n,2) epsrc_CC(n,3)] = Cole_Cole(freq(n),param_abs(1),param_abs(2),param_abs(3),param_abs(4),param_abs(5));
% end
% 
% 
% figure;
% semilogx(freq,Eps1,'k');
% grid on
% hold on
% semilogx(freq,epsrc_CC(:,2),'r');
% figure;
% semilogx(freq,Eps2,'k');
% grid on
% hold on
% semilogx(freq,epsrc_CC(:,3),'r');
