Let $S=k[x_1,...,x_9]$ be a polynomial ring over field $k$. Set $q_1=(x_1,x_2,x_5,x_6)$, $q_2=(x_1,x_2,x_6,x_7)$, $q_3=(x_2,x_3,x_7,x_8)$, $q_4=(x_1,x_5,x_6,x_7)$, $q_5=(x_1,x_6,x_7,x_8)$, $q_6=(x_2,x_6,x_7,x_8)$, $I=(x_1,..,x_4)$. How to compute $\operatorname{Depth}_I(S/\cap_{i=1}^{6}q_i)$? in general, Is there a way to compute $\operatorname{Depth}_I(S/\cap_{i=1}^{t}q_i)$ when $S=k[x_1,...,x_n]$, $q_i=(x_{i_1},...,x_{i_r})$? Is there any algorithm to compute it (using Macaulay2)?
background:
the common length of the maximal $M$-sequences in $I$ is called the grade of $I$ on $S$-module $M$ denoted by $\operatorname{Depth}_I(M)$.