// Ce script nécessite le module Atoms CASCI

clear;
clf;
// paramètres de la loi de Weibull

beta_forme = 0.845;
eta_echelle = 126;

// génération des données

Y = linspace(0, 1, 30)
Y = Y(2:$-1);

t_orig = floor(idfweibull(Y, beta_forme, eta_echelle))';

t = [2, 5, 9, 13, 17, 22, 27, 39, 39, 39, 52, 64, 64, 76, 86, 97, 108, 121,...
   135, 151, 169, 191, 215, 245, 282, 332]';
t_complet = [t ; 365 ; 365];
N0 = 28
tt = 14*(0:25);

// systèmes défaillants sur un intervalle

j = 1; // initialisation
// D(1) = 0
for i = 1:26 // par quinzaine
    D(i) = sum((t >= (i-1)*14)&(t < i*14))
end

// systèmes foncitonnant en entrée d'intervalle

N(1) = N0;
for i=2:26
    N(i) = N(i-1)-D(i-1);
end

// Probabilités conditionnelles

P = D./N;

// Fiabilité

for i = 1:26
    R(i)=prod(1-P(1:i));
end

// tracé

scf(0);
clf;

// plot(tt, [1; R(1:$-1)], "--");
plot2d2(tt, [1 ; R(1:$-1)]);
xtitle("Courbe de survie par la méthode actuarielle", "t (j)", "R")