(Subscript[K, e])={((B^T)(Subscript[D, s])Bif   (Subscript[ψ, i])(e)≤0(B^T)(Subscript[D, l])Bif   (Subscript[ψ, i])(e)>0)