I have the following integration:

Int(c-sqrt(a+b*(v*x+u)^2), x = -(b*u+sqrt(b*c^2-a*b))/(b*v) .. (-b*u+sqrt(b*c^2-a*b))/(b*v));

This integration is equivalent to the following integration (between dashed lines):

-----------------------------

NumericEventHandler(invalid_operation = `Heaviside/EventHandler`(value_at_zero = 1)):

Heaviside(0)

Heaviside(x) = convert(Heaviside(x), piecewise)

Int(Heaviside(x+(b*u+sqrt(b*c^2-a*b))/(b*v))*(c-sqrt(a+b*(v*x+u)^2)), x = -infinity .. (-b*u+sqrt(b*c^2-a*b))/(b*v));

-----------------------------

In this form the function, which is to be integrated, is taken as equal to zero for x < -(b*u+sqrt(b*c^2-a*b))/(b*v).

Now I want to write a form such that the function is also taken as equal to zero for x > (-b*u+sqrt(b*c^2-a*b))/(b*v).

How should I do this? Can this also be done with Heaviside?

With the following equation

eqn:=y=1/2+(1/2)*erf((1/2)*sqrt(2)*(x-mu)/sigma)-exp(-lambda*(x-mu)+(1/2)*lambda^2*sigma^2+ln(1/2-(1/2)*erf((1/2)*sqrt(2)*(lambda^2*sigma^2-lambda*(x-mu))/(lambda*sigma))));

and with

x:=solve(eqn,x) assuming sigma > 0, lambda > 0;

I got the following solution

x := -(1/2)*(-lambda^2*sigma^2-2*lambda*mu+2*RootOf(-exp(_Z)*erf((1/4)*sqrt(2)*(lambda^2*sigma^2+2*_Z)/(lambda*sigma))+exp(_Z)+erf((1/4)*sqrt(2)*(-lambda^2*sigma^2+2*_Z)/(lambda*sigma))+2*y-1))/lambda;

In order to get rid of RootOf I gave the command:

allvalues(%);

However, RootOf did not disappear. How should I proceed?