solve ln(s*p) = ln(c*10**12) - alpha*exp( -beta*ln(s / (24*6*365.25* 50/(2**(-0.471234+b/210000)))) ) for p