def set_p(p, Eg, W2, nmax, mmax):
p.ems = func_ems(Eg)
p.alpha = fetmodel.alpha_NP(Eg, p.ems)
p.W2 = W2
Cox = fetmodel.Cox_rect(epsOX, tOX, p.W1, W2)
Cc = fetmodel.Cc_rect(epsS, p.W1, W2)
p.Ceff = Cox * Cc / (Cox + Cc)
p.nmax = nmax
p.mmax = mmax
# Eg = 0.36
# epsOX = 8.5
# epsS = 8.9
# ems = 0.2
# tOX = 20e-9
Eg = 0.36
epsOX = 20
epsS = 15.15
ems = 0.023
tOX = 3e-9
temperature = 300
W1 = 30e-9 ## diameter
W2 = W1
alpha = fetmodel.alpha_NP(Eg, ems)
Cox = fetmodel.Cox_radial(epsOX, tOX, W1/2)
Cc = fetmodel.Cc_rect(epsS, W1, W1)
# zI = radius-Delta,
# Delta = \int rho(r) r^2 dr / \int rho(r) r dr
# Cc = fetmodel.Cc_radial(epsS, zI, W1-zI)
# See eq. in
# IEEE TRANSACTIONS ON ELECTRON DEVICES, VOL. 55, NO. 1, JANUARY 2008 411
# Modeling the Centroid and the Inversion Charge in
# Cylindrical Surrounding Gate MOSFETs,
# Including Quantum Effects
# J. B. Roldán, Andrés Godoy, Francisco Gámiz, Senior Member, IEEE, and M. Balaguer
nmax=5
mmax=5
print('Cox=', Cox,', Cc=', Cc)
p=fetmodel.parameters_ballistic(alpha=alpha,
Ceff=Cox*Cc/(Cox+Cc),
ems=ems,