expdecaydata;
close all;

p0 = [0.5, 6.0];
pest = lsqcurvefit(@(p, x) exp2func(x, p), p0, time, voltage);
fprintf('c=%4.1fmV, tau=%4.1fms\n', pest(1), pest(2));

hold on;
plot(time, voltage)
plot(time, exp2func(time, pest), 'LineWidth', 2)
hold off;
xlabel('time [ms]')
ylabel('voltage [mv]')
savefigpdf(gcf, 'exp2fit.pdf', 15, 8.5);