(                (Subscript[s, 1,(Subscript[b, k])])=(∑_((x,y)∈(Subscript[b, k])))(∑_([i',j]∈Ω(x,y)))((((Subscript[τ, 1])∂(I^((1)))/∂x+(Subscript[τ, 0])∂(I^((0)))/∂x))^2);       (Subscript[s, 3,(Subscript[b, k])])=(∑_((x,y)∈(Subscript[b, k])))(∑_([i',j]∈Ω))((I^((1)))-(I^((0))))((Subscript[τ, 1])∂(I^((1)))/∂x+(Subscript[τ, 0])∂(I^((0)))/∂x);
(Subscript[S, 2,(Subscript[b, k])])=(∑_((x,y)∈(Subscript[b, k])))(∑_([i',j]∈Ω))((Subscript[τ, 1])∂(I^((1)))/∂x+(Subscript[τ, 0])∂(I^((0)))/∂x)((Subscript[τ, 1])∂(I^((1)))/∂y+(Subscript[τ, 0])∂(I^((0)))/∂y);
                (Subscript[s, 5,(Subscript[b, k])])=(∑_((x,y)∈(Subscript[b, k])))(∑_([i',j]∈Ω(x,y)))((((Subscript[τ, 1])∂(I^((1)))/∂y+(Subscript[τ, 0])∂(I^((0)))/∂y))^2);
                (Subscript[s, 6,(Subscript[b, k])])=(∑_((x,y)∈(Subscript[b, k])))(∑_([i',j]∈Ω))((I^((1)))-(I^((0))))((Subscript[τ, 1])∂(I^((1)))/∂y+(Subscript[τ, 0])∂(I^((0)))/∂y)Eq.  (14))