nested_lsqnonlin.m

clc; clear;

a0 = 0;
b0 = 0;

solution = lsqnonlin(@(x) difffun_a_b(x), [a0, b0]);

disp(solution);

function difference = difffun_a(a,b)
    
    difference = (a*2 + b) - (1);
    disp('diffun_a')
    disp(difference);

end

function difference = difffun_a_b(x)
    
    a0 = x(1);
    b0 = x(2);

    difference = lsqnonlin(@(x) difffun_a(x,b0), a0);
    disp('diffun_a_b')
    disp(difference);
end