oc = load("occupancy_MM1_1.txt");
t = oc(:,1);
nw = oc(:,2);
ns = oc(:,3);
t_n = [t,nw+ns];
pr=probabilidad_estado(t_n)
n=(1:length(pr));
fun_geom=(1/2).^(n);
error=min(pr'-fun_geom);