import  scipy.optimize as opt

def reaktor(xlist):
    ksi1 = 490
    k1 = 0.016

    n_co = xlist[0]
    n_co2 = xlist[1]
    n_h2o = xlist[2]
    n_h2 = xlist[3]
    ksi2 = xlist[4]

    eq0 = k1*n_co2*n_h2 -(n_co*n_h2o)
    eq1 = n_co - (560 + ksi2 - ksi1)
    eq2 = n_h2o - (ksi2)
    eq3 = n_co2 - (140 - ksi2)
    eq4 = n_h2 -(5250 - 2*ksi1 - ksi2)

    balance = [eq0,eq1,eq2,eq3,eq4]

    return balance

gList=[200,200,200,200,200]
gkomp = ['co','co2','h2o','h2']

ans =opt.root(reaktor,gList)

print(ans.x)

n_ch4 = 1050

tot = n_ch4 + sum(ans.x)

for k,v in zip(gkomp,ans.x):
    print(k, ' : ', v/tot)

print(tot)