Add here a paragraph of the function description.
Te = 1 Te_rec = 1 Tfin = 100 Nbit = 16 a=0.5 b=-1 ci=0.675 Intmax = 2^(Nbit-1); ci1=-int(ci*Intmax) ci2=int(ci*Intmax) p_source=2*%pi*0.01 a_source=2^(Nbit-1)-1