clc //initialisation of variables n= 2 C= 0.01 //M //CALCULATIONS r= 10^(-0.509*n*sqrt(C)) //RESULTS printf ('mean ionic activity coefficient = %.2f ',r)