Untitled

def CalculateFlowX35(x23):  
    x35 = x23 * (1/200);
    return x35;











x= 14/1000
m_alucans = x*((7*0.1)/0.5) + x*((93*0.1)/0.5) + x*(50/0.5)
m_alucans










flow = 'x34'








A = np.array([[ 1, -1, -1, 0, 0, 0, 0, 0, 0, 0, 1], #1
              [ 0, 1, 0, -1, -1, 0, 0, 0, 0, 0, 0], #2
              [ 0, 0, 0, 1, 0, -1, -1, 0, 0, 0, 0], #3
              [ 0, 0, 0, 0, 0, 1, 0, -1, -1, 0, 0], #4
              [ 0, 0, 1, 0, 1, 0, 1, 0, 1, -1, -1], #5
              [ -0.72, 1, 0, 0, 0, 0, 0, 0, 0, 0, -0.72], #SHEET
              [ 0, -0.79, 0, 1, 0, 0, 0, 0, 0, 0, 0], #CAN
              [ 0, 0, 0, (-199/200), 0, (1), 0, 0, 0, 0, 0], #PD 
              [ 0, 0, 0, 0, 0, (-0.85), 0, 0, (1), 0, 0], #PD
              [ 0, 0, -0.975, 0, -0.975, 0, -0.975, 0, -0.975, 0, 1],  #PD
              [ 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0]]  #PD
             )

y=np.array([[0],
            [0],
            [0],
            [0],
            [0],
            [0],
            [0],
            [0],
            [0],
            [0],
            [m_alucans]])


x = np.linalg.inv(A)@y
print (A)
print (x)




















x51 = x[10]

















loss_share_consumers = (x[7]/(x[9]+x[7]))
loss_share_consumers




















x01 = x[0]
print (x01)



















x01_new = 0.28863762