Matlab code for ternary model for unbleached molecules seems to work properly now. Problem before: volume fraction not defined properly and thus low mobility inside.