It depends on what you mean by "For each value of Jw and mu_o". If you mean for each of eight corresponding pairs of numbers from these two vectors, do this:
for k = 1:8
Hw(k) = fsolve(@(Hw) .....mu_o(k).....Jw(k).....
end
If you mean for each of the sixty-four possible combination of pairs from these two arrays, then do:
for k1 = 1:8
for k2 = 1:8
Hw(k1,k2) = fsolve(@(Hw) .....mu_o(k1).....Jw(k2).....
end
end
Let us hope that the initial estimate of 0.9 will lead to successful convergence for all cases.
There is also the possible problem of multiple solutions for this equation.
