default(realprecision,38);
nf = nfinit('y^3 - 2);
nfeltembed(nf, y)
nfeltembed(nf, y, [2])
nfeltembed(nf, y, [1,2])
nfeltembed(nf, 2, [1..3])

nf=nfinit(x^4-3);
[nfeltembed(nf,x,i) | i <- [1..3]]

\\Errors
nfeltembed(nf, y, [1,2,4])

