p = boltzmann(voltage, a);
plot(voltage, p)
y = -1;