function [tanho,eanho] = estado(L,R,Na)
% [tanho,eanho] = estado(L,R,Na)
%
% Funcion que nos devuelve dos vectores relacionados con el funcionamiento
% de un sistema reparable (p.ej. una linea electrica).
% Dadas una tasa de fallos (L) y tiempo medio de % reparacion (R), 
% la funcion nos devuelve el estado de la linea (eanho) cada hora 
% durante un numero de aņos (Na).
%
% DATOS ENTRADA:
%  - L : Tasa de fallos (fallos/anho)
%  - R : Tiempo medio de reparacion (horas/fallo)
%  - Na: Numero de aņos
%
% DATOS SALIDA:
%  - tanho: eje de tiempos con Na*8760 valores (nš de horas en Na aņos)
%  - enaho: estado de la linea cada hora (1: funcionamiento; 0:reparacion)
%

%Kmi 2005

if nargin~=3
    error('Nš incorrecto de parametros de entrada')
end
if Na <=0.9999 
    error('Nš de aņos ha de ser un entero mayor que uno')
end    

R = R/8760;

tanho = [0:1:8760*Na];
N = length(tanho);
eanho = zeros(N,1);

Estado = 0;
Alcance = 0;
taux = 0;
for k=1:length(tanho)
   if tanho(k) >= taux
       if Estado == 0
           T = exprnd(1/L)*8760;
           taux  = tanho(k) + T;
           Estado = 1;
       else   
           T = exprnd(R)*8760;           
           taux = tanho(k) + T;
           Estado = 0;
       end
   end    
   eanho(k) = Estado;        
end