d/dn ((1/2)*sqrt(pi)/sqrt(b)*exp(-(a-1)^2/(4*b))*erfi((2*b*log(n)-a+1)/(2*sqrt(b))))