clc

ET=90//mJ/cm^2
EI=45//mJ/cm^2
gammma=1/[log(ET/EI)]
disp(gammma,"gamma is= ")
ET=7//mJ/cm^2
EI=12//mJ/cm^2
gammma=1/[log(EI/ET)]
disp(gammma,"gamma is= ")