### Com gnuplot
## f(x) = x*sin((x + 6)/(x + 13))

u(x) = x + 6
du(x) = 1
v(x) = x + 13
dv(x) = 1
g(x) = u(x)/v(x)
dg(x) = ( v(x)*du(x) - u(x)*dv(x) )/( v(x)*v(x) )
f(x) = x*sin(g(x))
df(x) = sin(g(x)) + x*cos(g(x))*dg(x)
set xrange [-10:10]
plot f(x),  0


### com calc

define f(x) {return  x*sin((x + 6)/(x + 13));}

epsilon = 0.01
x = -10
delta = 1
while( x <= 10) {
	if ( f(x)*f(x+delta) <= 0 )
	{
		if (f(x) == 0)
			printf("%f  é raíz exata de f ", x);
			else if (f(x+delta) == 0)
				printf("%f  é raíz exata de f ", x+delta);
		else 
			printf("troca de sinal no intervalo [%f, %f] \n", x, x + delta);
		/* x = x + delta;  elimine este comentário  */
	}
	x = x + delta;
}