how to convert this coding to broyden method coding

x0=[1 1 1]';

x =x0;

for it=0:2

it=it+1;

F=[10*x(1)-2*x(2)^2+x(2)-2*x(3)-5 8*x(2)^2+4*x(3)^2-9 8*x(2)*x(3)+4];

J=[10 -4*x(2)+1 -2; 0 16*x(2) 8*x(3); 0 8*x(3) 8*x(2)];

J_inv=inv(J);

y=-J_inv*F';

x=x+y;

string=['x(',num2str(it),'): ',num2str(x')];

disp(string)

end

Hi guys this is the coding for newton method but i have to use the same coding and change it to broyden method so can you help me...please