clc; clearvars;
u = 0.5:0.5:1.5;
lu = length(u);
t = 0:0.1:5;
lt = length(t);
Pt = zeros(lt,lu);
for m = 1:lu
    mu = u(m);
for i = 1:lt
    time = t(i);
    Pt(i,m) = exponential(mu,time);
    
end
end
plot(t',Pt);
xlabel('Time')
ylabel('P(t)')
legend({['u = ',num2str(u(1))],['u = ',num2str(u(2))],...
    ['u = ',num2str(u(3))]}, 'Location','northeast')
grid on