要统计的是每个周期内达到某一位置处的电子流
下面是大致思路：
只需要设置变量fe， 统计每个步长dt到达这一位置的电荷ne 之后，在外部fe=fe+ne
注意 每次统计ne和nion之前，要将ne=0,nion=0.即ne和nion只是统计一个步长的通量之后在每个周期尾部：fe=0和 fion=0