f = zeros(Ntime, dtype=float)