$qn2= -$cys/(1+10^(9.0-$ph));