function Cn = Cnew(nc,C)
Cn=C;
for i=1:nc
 if ((C(i)==2) && (rand() < 0.2))
  Cn(i)=0.2;
 end
end
end