function model=MTRL_RR(Km,label,task_index,insIndex,ins_num)

    n=size(Km,1);

    task_num=length(insIndex);

    M12=zeros(n,task_num);

    for i=1:task_num

        M12(insIndex{i},i)=1;

    end

    tmp=[label,zeros(1,task_num)]/([[Km+diag(ins_num(task_index))/2,M12];[M12',zeros(task_num)]]);

    model.alpha=tmp(1:n);

    model.b=tmp((n+1):(n+task_num));

    clear tmp M12;

    clear Km label task_index insIndex ins_num;