A = rand(10,1);
B = rand(10,1);
f = @(a,b)@(x)[x(1).^a - b; x(2).^b - a];
n = numel(A);
op = optimset('display','off');
x = rand(2,n);
for i = 1:n
x(:,n) = fsolve(f(A(i),B(i)),rand(2,1),op);
end
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!