close; clear;

figure(1 ,'position',[0,0,845,400]);
t=0:0.1:25;

% Mortalità Infantile
beta = 0.5;
eta  = 10;
y1   = (beta/eta).*(t/eta).^(beta-1).*exp(1).^(-(t/eta).^beta);
p1 = plot(t,y1,"linestyle","-.","color","k",'linewidth',2);
hold on;


% Guasto aleatorio
beta = 1;
eta  = 10;
y2   = (beta/eta)*(t/eta).^(beta-1).*exp(1).^(-(t/eta).^beta);
plot(t,y2,"linestyle","-","color","r",'linewidth',2);
hold on;

% Usura iniziale
beta = 3.66;
eta  = 10;
y3   = (beta/eta)*(t/eta).^(beta-1).*exp(1).^(-(t/eta).^beta);
plot(t,y3,"linestyle","--","color","b",'linewidth',2);
hold on;


% Usura avanzata
beta = 5;
eta  = 10;
y4   = (beta/eta)*(t/eta).^(beta-1).*exp(1).^(-(t/eta).^beta);
plot(t,y4,"linestyle",":","color","m",'linewidth',2);


% Opzioni plot
grid on;
axis("auto");

xlabel("Invecchiamento","fontsize",22); ylabel("Tasso di guasto","fontsize",22);
axis([0,20,0,0.35]);
legend("Mortalità Infantile","Guasto aleatorio","Usura iniziale","Usura avanzata","fontsize", 18);
title("Probability Density Function", "fontsize", 32);

print('PDF.jpg','-opengl','-S845,500','-FEras ITC:20');