function N=N_functie(L,A,r);

n=1;

N=A.*(1./r-1./(r+L));
