0001 function dist = sigma_tanh(nlev,dl,du) 0002 0003 z=0.0; 0004 for k=1:nlev 0005 x1=dl+du; 0006 x1=x1*(nlev-1-k)/(nlev-1); 0007 x1=x1-dl; 0008 x1=tanh(x1); 0009 x2=tanh(dl); 0010 x3=x2+tanh(du); 0011 z(k+1)=(x1+x2)/x3-1.0; 0012 end; 0013