#Vegard Jervell, 18.01-18
#Dette scriptet løser et likningssett som beskriver en mixer

from scipy import optimize as opt

def mix_modell(x_list):

    #definerer kjennte variabler
    w1_HNO3 = 0.23 #vektprosent
    w2_HNO3 = 0.9
    w4_HNO3 = 0.27

    w1_H2SO4 = 0.57 #vektprosent
    w3_H2SO4 = 0.93
    w4_H2SO4 = 0.6

    m4 = 100 #kg/s

    #definerer ukjente variabler
    m1 = x_list[0]
    m2 = x_list[1]
    m3 = x_list[2]

    #setter opp likninger

    eq1 = m1*w1_HNO3 + m2*w2_HNO3 - m4*w4_HNO3
    eq2 = m1*w1_H2SO4 + m3*w3_H2SO4 - m4*w4_H2SO4
    eq3 = m1 + m2 + m3 - m4

    balance = [eq1,eq2,eq3]

    return balance

#Mine gjettninger på m1,m2 og m3
Gm1 = 50
Gm2 = 10
Gm3 = 10
Glist = [Gm1,Gm2,Gm3]
#bruker scipy.optimize.root for å løse likningssettet
ans = opt.root(mix_modell, Glist)