% calib
clear;
clc;
close all;
qrupt=[0;0;0];
for i=1:6
    nomfic=['fig',num2str(i),'.txt'];
    M=load(nomfic);
    hold on
    if i<=3
        qrupt(i)=max(M(:,2));
        figure(1)
        plot(M(:,1),M(:,2),'-*');
    else
        figure(2)
        plot(M(:,1),M(:,2),'-*');
    end
    
end
figure(1)
grid on
legend('p_0=50 kP_a','p_0=100 kP_a','p_0=200 kP_a','Location','northwest')
xlabel('\epsilon_1 (%)')
ylabel('q (kP_a)')

figure(2)
grid on
xlabel('\epsilon_1 (%)')
ylabel('\epsilon_v (%)')


%%%%%%%%%%%%%%%%%%%%%%
% C (kPa) et Phi (degre)
% regression lineaire 
s3rupt = [50;100;200];
% moindres carres donne :



%%%%%%%%%%%%%%%%%%%%%%
% E50 (en Mpa)

%%%%%%%%%%%%%%%%%%%%%%
% nu tangent (ur ref)


%%%%%%%%%%%%%%%%%%%%%
% psi (en degre)


%%%%%%%%%%%%%%%%%%%%%%
% puissance m (pour pref= 100 kPa)