function callPhiII

T0 = 1000.0;
n0 = 5e10;
V0 = -5000.0;
arrT0 = [1000.0 2000.0 3000.0 4000.0 5000.0];
arrN0 = [1e10 5e10 1e11 5e11];
arrV0 = [-1000.0 -2000.0 -3000.0 -4000.0 -5000.0 -10000.0];
nN0 = length(arrN0);
nV0 = length(arrV0);
nT0 = length(arrT0);

RN = 0.004;
kB = 1.3806503e-23; % m2 kg s-2 K-1
e = 1.60217646e-19; % C
eps = 8.85418782e-12; % m-3 kg-1 s4 A2
iN0 = 1;
iT0 = 1;
iV0 = 1;
for iN0 = 1:nN0
    for iT0 = 1:nT0 
        for iV0 = 1:nV0
            T0 = arrT0(iT0);
            N0 = arrN0(iN0);
            V0 = arrV0(iV0);
            [RS,pIIRS,rII,pII] = phiII(T0,N0,V0);
            %[RS,pIIRS,rII,pII] = phiII(T0,N0,V0);
            x = sqrt(eps*kB*T0/(e*e*N0));
            fprintf(1,'%g,%.0f,%.0f,%.2f,%f,%.4f\n',N0,T0,V0,RS,x/RS,pIIRS);
        end
        end
        end
%plot(rII,pII)


