isDeficient -- determines deficiency of a network

i1 : N = reactionNetwork "A <--> 2B, A + C <--> D, B + E --> A + C, D --> B+E"

o1 = A-->2B
     2B-->A
     A+C-->D
     D-->A+C
     D-->B+E
     B+E-->A+C

o1 : ReactionNetwork

i2 : isDeficient N

o2 = 0

i3 : W = wnt()

o3 = X_1-->X_2
     X_2-->X_1
     X_2-->X_3
     X_2+X_4-->Y_4
     Y_4-->X_2+X_4
     Y_4-->X_2+X_5
     X_5+X_8-->Y_6
     Y_6-->X_5+X_8
     Y_6-->X_4+X_8
     X_4+Y_0-->Y_8
     Y_8-->X_4+Y_0
     Y_8-->X_4
     0-->Y_0
     Y_0-->0
     Y_0-->Y_1
     X_3+X_6-->Y_5
     Y_5-->X_3+X_6
     Y_5-->X_3+X_7
     X_7+X_9-->Y_7
     Y_7-->X_7+X_9
     Y_7-->X_6+X_9
     X_6+Y_1-->Y_9
     Y_9-->X_6+Y_1
     Y_9-->X_6
     Y_1-->0
     Y_1-->Y_0
     Y_1+Y_2-->Y_3
     Y_3-->Y_1+Y_2
     X_3-->X_2
     X_5-->X_7
     X_7-->X_5

o3 : ReactionNetwork

i4 : isDeficient W

o4 = 4