xdel(winsid())
clear;

t = [0:0.01:10];
a = 0.5;
b = 0.6;
k = 1;

y = k*(exp(-a*t) - exp(-b*t) );
ymax = max(y)*ones(1,length(y));

xpts = [1 1];
ypts = [0 max(y)];
plot(xpts*0.07, ypts, '--');
plot(xpts*1.1, ypts, '--');
plot(xpts*8.96, ypts, '--');

plot(t,y, 'black');
plot(t,ymax.*0.1, '--');
plot(t,ymax.*0.9, '--');

// titletxt = ['$y = k(e^{-\alpha t} - e^{-\beta t})$' ; strcat(['$k=', string(k), ', \alpha=', string(a), ', \beta=', string(b), '$']) ];

titletxt = '$u(t) = k(e^{-\alpha t} - e^{-\beta t})$';
title(titletxt , 'fontsize', 8);
xlabel('t (s)', 'fontsize', 7);
ylabel('U (V)', 'fontsize', 7, 'rotation', 0);

a = get("current_axes");
a.axes_visible = "off";
