k0 = 30.0297;

k = 30.096;



E = linspace(2,6);



singleE = k*E.^(-2);

doubleE = k0*(0.86*E.^(-2)+0.14*(E-0.043).^(-2));



plot(E,180*singleE);

hold on;

plot(E,180*doubleE,'r--');