why is my code not working?
显示 更早的评论
I'd like to get the index (i) for wich the condition I set in the foor loop appears on the results
Fm0= linspace(0.01,0.02,100);
%ciclo iterativo
for i=1:100
for Fm= (0.25).*(((log((2.7*10^-5)+(2.51.*((Rey.*sqrt(Fm0).^-1))))).^-1).*((log((2.7*10^-5)+(2.51.*((Rey.*sqrt(Fm0).^-1)))))).^-1)
if Fm(i)-Fm0<=10^-6
break
end
end
end
采纳的回答
Walter Roberson
2019-5-18
Fm0 = linspace(0.01,0.02,100);
Fm= (0.25).*(((log((2.7*10^-5)+(2.51.*((Rey.*sqrt(Fm0).^-1))))).^-1).*((log((2.7*10^-5)+(2.51.*((Rey.*sqrt(Fm0).^-1)))))).^-1)
i = find( Fm-Fm0 <= 10^-6, 1);
You will find that the answer is 1 -- Fm0(1) > Fm(1) so the subtraction is negative and so satisfies the condition.
You might perhaps have been thinking,
i = find( abs(Fm-Fm0) <= 10^-6, 1);
but you will find that this is empty.
2 个评论
Franco Belletti
2019-5-19
Walter Roberson
2019-5-19
If you are looking for the value of Fm such that
Fm = (0.25).*(((log((2.7*10^-5)+(2.51.*((Rey.*sqrt(Fm).^-1))))).^-1).*((log((2.7*10^-5)+(2.51.*((Rey.*sqrt(Fm).^-1)))))).^-1)
then the value turns out to be complex, at least for positive Rey (I did not experiment with negative Rey.)
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Mathematics 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
