% fun1.m
% example function for fsolve
% Three equations in three unknowns
% exact solution at [-1,1,3/2]
%
function f = fun1(x)
f(1) = x(1).^2 + x(2).^2 + x(3).^2 - 17/4;
f(2) = 2*x(1)*x(2)+3*x(2)*x(3)+x(3)*x(1).^4 - 4;
f(3) = x(1).^2 + x(2).^4 + x(3).^4 - 113/16;