d((1-f)^2*s^2 + f^2*h^2 + 2*f*(1-f)*c*s*h)/df