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