M = 100;
nu = linspace(-.0001,.0001,M);
c = 2*nu./log(1+(2*nu./(1-nu));
cl = ones(M,1);
plot(nu,c, nu,cl)