:- use_module(library(clpBNR)).
global_maximum/2
clpBNR_default_precision
?- X::real(0,3r4*pi), global_maximum(X*sin(4*X),Z,4). X:: 1.99..., Z:: 1.979... .