derive 67.66*sqrt(0.0000918/a^4+0.3110082/a^3+0.6889)