For loops for multiple variables
显示 更早的评论
Hello there!
I want to generate 3 graphs of y vs x with three different boundaries (C). However, I need to assume that changing C also changes 2 other parameters (m and rho), I mean for each of C there are certain values for those two parameters. In reality, I have to take 3 graphs, but I obtain around 27 graphs. What I made so far is the following:
C = [0.342 0.599 0.856];
m =[0.000755 0.000870 0.000969]; % viscosity [Pa s]
rho = [1022 1028 1035]; % density of seawater [kg m-3]
D = 1.5e-9; % diffusion coefficient [m2 s-1]
k = 1.42e-5; % membrane permeability [m2 s-1]
for i = 1:length(C)
for j = 1:length(m)
for l = 1:length(rho)
b = - m(j) / (rho(l)*k*D);
fun1 = @(x,Y)[Y(2);-(Y(2).*(x.^2.*b+1.0))./x];
R1 = 1.63e-4; % inner radius [m]
R2 = 2.5e-4; % outer radius [m]
bc = @(ya, yb) [ya(1) - (C(i)-0.9998*C(i)); yb(1)- C(i)];
xmesh = linspace(R1, R2, 5000);
solinit = bvpinit(xmesh, [1 0]);
sol = bvp4c(fun1, bc, solinit);
figure(1);
plot(sol.x, sol.y(1,:),'LineWidth',2); hold on;
set(gca, 'XDir','reverse')
end
end
end
Thanks in advance. Any comment is highly appreciated.
采纳的回答
You only need a single loop:
for n = 1:length(C)
m_=m(n);
rho_=rho(n);
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
