from matplotlib.pyplot import *
from numpy import *

r = linspace(0.01e-6, 10e-6, 10e3)

#kelvins formel:
sigma = 0.072*10e-3     #J/m^2
n     = 3*10e25 #antall vannmolekyler per liter vann
k     = 1.38*10e-23     #J/K mol
T     = 273.15      #kelvin


i_naci = 2
i_amm = 3
Mw = 18.015     #g/mol
MnaCl = 58.44   #g/mol
Mamm = 132.13952  #g/mol
rho_na = 650   #Kg/m^3
rho_amm = 750   #3Kg/m^3

a = (2*sigma/(n*k*T))

m1 = m4 = 10e-19
m2 = m5 =10e-18
m3 = 10e-17
m6 = 10e-20


b1 = (i_naci*m1*Mw)/(MnaCl*((4/3)*pi*rho_na))
RH1 = 1 + a/r - b1/r**3

b2 = (i_naci*m2*Mw)/(MnaCl*((4/3)*pi*rho_na))
RH2 = 1 + a/r - b2/r**3

b3 = (i_naci*m3*Mw)/(MnaCl*((4/3)*pi*rho_na))
RH3 = 1 + a/r - b3/r**3

b4 = (i_amm*m4*Mw)/(Mamm*((4/3)*pi*rho_amm))
RH4 = 1 + a/r - b4/r**3

b5 = (i_amm*m5*Mw)/(Mamm*((4/3)*pi*rho_amm))
RH5 = 1 + a/r - b5/r**3

b6 = (i_amm*m6*Mw)/(Mamm*((4/3)*pi*rho_amm))
RH6 = 1 + a/r - b6/r**3


#Ex.1d)
def funk(r):
    return (1 + (a/r))*100

#Ex.2c)
def f1(r):
    return RH1

def f2(r):
    return RH2

def f3(r):
    return RH3

def f4(r):
    return RH4

def f5(r):
    return RH5

def f6(r):
    return RH6


semilogx(r/10e-6, funk(r))
hold(True)
semilogx(r/10e-6, f1(r))
hold(True)
semilogx(r/10e-6, f2(r))
hold(True)
semilogx(r/10e-6, f3(r))
hold(True)
semilogx(r/10e-6, f4(r))
hold(True)
semilogx(r/10e-6, f5(r))
hold(True)
semilogx(r/10e-6, f6(r))
ylim(1,1.0004)


show()