1 function [DLL] = RadialDiffCoeff(Kp,L)
2 %RadialDiffCoeff - This function returns the radial diffusion coefficients
3 %for L between 2 and 8 and for a provided value of Kp.
4 % This function enables one to reproduce the first row of figures in this paper
6 DLL = 10^(0.506*Kp - 9.325)*(L.^10);
8 %To see plot, uncomment following lines.
11 %plot(L,1./DLL(2,L),
'.-k');
12 %set(gca,
'yscale',
'log',
'ylim',[1e-2 1e6],
'ylabel',
'1/Dll',
'xlabel,'L
');