#N_POINTS=3
tscan(charge, N_POINTS, 1.0)