experiment :=
{
 pops : list of pop;
 mean : real := pluslist(pops.x)/length(pops.x);
 outmean : list of real := append(outmean_1, mean);
}

pop := 
{
 x       :  real
         := 2 * haselements (persons.att,1) 
            / length (persons) 
            - 1;
 persons :  list of person;
 outx    :  list of real := append (outx_1,x)
};
person :=
{ 
 att :  int
     := case att_1 of
        1 :0 if 
             uniform(1,0.0,1.0) <
             nu * exp(-(pi + kappa * elem(p.x_1,1)))
           else 1;
        default:1 if 
             uniform(1,0.0,1.0) <
             nu * exp(  pi + kappa * elem(p.x_1,1)) 
           else 0;
        end;
 p   :  list of pop
};