How can i use simulated annealing algorithm in matlab toolbox

Hello, i have a problem in calling my function in fsolve or simulannealbnd. My function looks like this:

function Fun=MyFun(q,k)

h = 1;

A = [-1 -3;2 -5];

Ad = [1.66 -0.697;0.93 -0.330];

Q = [q(1) q(3);q(2) q(4)];

H = Ad*h*Q;

[V,D]=eig(H);

Dw = diag(lambertw(k,diag(D)));

F = V*Dw/V*expm(V*Dw/V+A)-Ad;

Fun=[F(1);F(2);F(3);F(4)];

How can i use these function toolbox to solve my problem above. For instance, k=0. Thanks for you help.