% T = taskset([1 1 1 2 2 2 3 3 4]);
% T.Name = {'t1' 't2' 't3' 't4' 't5' 't6' 't7' 't8' 't9' };
% p = problem('P|pmtn|Cmax');
% TS = mcnaughtonrule(T,p,3);
% plot(TS);
t1=task('t1',2);
t2=task('t2',4);
t3=task('t3',1);
t4=task('t4',3);
t5=task('t5',2);
t6=task('t6',1);
t7=task('t7',3);
t8=task('t8',1);
t9=task('t9',2);

prec = [0 1 0 0 0 0 0 0 0 ;...
        0 0 0 1 0 0 0 0 0 ;...
        0 0 0 1 0 0 0 0 0 ;...
        0 0 0 0 1 1 1 0 0 ;...
        0 0 0 0 0 0 0 0 1 ;...
        0 0 0 0 0 0 0 0 1 ;...
        0 0 0 0 0 0 0 1 0 ;...
        0 0 0 0 0 0 0 0 1 ;...
        0 0 0 0 0 0 0 0 0 ];

T = taskset([t1 t2 t3 t4 t5 t6 t7 t8 t9],prec);
p = problem('P|prec|Cmax');
TS = listsch(T,p,2);
plot(TS);