S = 10;
TH = 7;
[Mu]e = 60;
[Mu]g = 30;
[Lambda]e = 30;

condition[S_, TH_][i_, j_] := sosa[i, j] Boole@Or[
[Mu]e*p[i, j] == [Lambda]e*p[i - 1, j] /; j == 0 && i == S (*No1*)
([Lambda]e + [Mu]e)*p[i, j] == [Mu]e*p[i + 1, j] + [Mu]g*p[i, j + 1] + [Lambda]e*p[i - 1, j] /; j == 0 && TH < i && i < S(*No2*)
([Lambda]e + [Mu]e)*p[i, j] == [Mu]e*p[i + 1, j] + [Mu]g*p[i, j + 1] /;j == 0 && i == TH(*No3*)
([Mu]e + [Mu]g)*p[i, j] == [Lambda]e*p[i - 1, j] /;j == S - i && S - TH <= i && i < S(*No4*)
([Lambda]e + [Mu]e)*p[i, j] == [Mu]e*p[i + 1, j - 1] + [Mu]e*p[i + 1, j] + [Mu]g*p[i, j + 1] /; j == TH - i && 0 < i && i < TH(*No5*)
[Lambda]e*p[i, j] == [Mu]e*p[i + 1, j - 1] + [Mu]e*p[i + 1, j] /; j == TH && i == 0(*No6*)
([Lambda]e + [Mu]e + [Mu]g)*p[i, j] == [Mu]e*p[i + 1, j] + [Lambda]e*p[i - 1, j] /; j == TH && 0 < i && i < S - TH(*No7*)
([Lambda]e + [Mu]e + [Mu]g)*
 p[i, j] == [Mu]e*p[i + 1, j] + [Lambda]e*p[i - 1, j] + [Mu]g*p[i, j + 1] /; 0 < j && j < TH && (TH - j) < i && i < (S - j)(*No8*)];
sosalist = SparseArray[{i_, j_} :> condition[S, TH][i - 1, j - 1], {S + 1,TH + 1}]