d/dx of u[x](a*cos(u[x]) + b*cot(u[x]))