MATLAB: How to apply ‘for’ cycle for two arrays

programming by matlab

Hello,
I have these parameters:
D=0.030;
A=7.0685*10^-4;
Jo=0.763;
Jw=[0.778090833 0.998156927 1.194644511 1.395061847 1.583689928 1.772318008 1.964875841 2.369640264];
ro_o=910;
ro_w=1000;
mu_w=0.001;
mu_o=[1.057 1.079 1.057 1.117 1.086 1.050 1.057 0.993];
and the function as:
Hw = fsolve(@(Hw) (((8*ro_o*((D-(2*sqrt(Hw*A/pi)))*(Jo/(1-Hw))*ro_o/mu_o)^-1.0)*(((Jo/(1-Hw))-(Jw/Hw))*((Jo/(1-Hw))-(Jw/Hw))))*(pi*(D-(2*sqrt(Hw*A./pi)))*(1/(1-Hw))))-((0.023*(((Jw*D*ro_w)/mu_w)^-.2)*ro_w*(Jw/Hw)^2*(pi*D))), 0.9)
For each value of Jw and mu_o, there must be one value of Hw, how can I write 'for' loop to take into consideration of two vectors Jw and mu_o?

Best Answer

  • 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.