function f = idos_n(E,n)
f = 

