Here is an example of a simple 1D function that you can run BO directly to find the optimal.
The tricky part is to define the function handle that can deal with the BO variable, which is implicitly explained in the document.
var = optimizableVariable('theta',[1,10],'Type','real');
fun = @(x)objfunxx(x.theta); %Create Function Handles
results = bayesopt(fun,var)
function y = objfunxx(x)
y=x*sin(x);
end
