read qMultiZeilberger:

qMulZeil(qfac(n)/qfac(i)/qfac(j)/qfac(n-i-j),
q^i+q^j,binomial(i,2)+binomial(j,2),[i,j],n,N,{});

quit: