function E = Emod(ex)
load mu_lamda_gafe16.ws
strainx = Sigma1/56e9;
strain_tot = strainx + Lambda_S;
modulus = Sigma1./strain_tot;
modulus(301) = (modulus(300)+modulus(302))/2;
E = interp1(strain_tot,modulus,ex,'linear');
end