function q=CreateRandomSolution(model,n,s)

    I=model.I;    
    q=randperm(I(s)+n-1);

end