function q = xyz(p)
% function q = xyz(p)   Example function for testing fsolve
x = p(1); y = p(2); z = p(3);
q = zeros(3,1);
q(1) = sin(x) + y^2 + log(z) - 7;
q(2) = 3*x + 2^y - z^3 + 1;
q(3) = x + y + z - 5;
