Moron> (c665) 'at('diff(eta(q),q,1),q = %e^-(3*%pi)) [...] No! Just e^-pi: (c675) ('at('diff(eta(q),q,1),q = %e^-%pi) = -%e^(3*%pi)*thetaderiv[1](%pi/3,2,%e^-(%pi/6))/(24*sqrt(3))) = %e^(3*%pi)*(16*gamma(1/4)-gamma(1/4)^5/%pi^2)/(64*2^(7/8)*%pi^(7/4)) | d | (d675) -- (eta(q))| dq | - %pi |q = %e 3 %pi %pi - %pi/6 %e thetaderiv (---, 2, %e ) 1 3 = - -------------------------------------- = 24 sqrt(3) 5 1 gamma (-) 3 %pi 1 4 %e (16 gamma(-) - ---------) 4 2 %pi = --------------------------------- 7/8 7/4 64 2 %pi (c676) dfloat(%); | d | (d676) -- (eta(q))| dq | |q = 0.04321391826377d0 = - 77.8356998979772d0 = - 77.8356998979791d0 --rwg