fixed simulation of concentration matrix