# Variables
import math
s = 1.771*10**-4 #mole litre**-1
s1 = 0.3252*10**-2 #mole litre**-1
# Calculations
S = s*10**(0.509*(math.sqrt(s+s1)-math.sqrt(s)))
# Results
print 'Solubility = %.2e mole litre**-1'%(S)
# Variables
import math
a = 0.1
# Calculations
r = 10**(-0.509*math.sqrt(a)/(1+math.sqrt(a)))
# Results
print 'mean ionic acctivity coefficient = %.3f '%(r)