expr = 1 / cos(10 * x) + 5 * sin(x)
plot(expr, ylim=(-10, 10), exclude=solveset(cos(10 * x)))
# Expected:
## Plot object containing:
## [0]: cartesian line: 5*sin(x) + 1/cos(10*x) for x over (-10.0, 10.0)