Info

This question is closed. Reopen it to edit or answer.

if Loop to print values into a variable

2 views (last 30 days)
chong kai sheng
chong kai sheng on 4 Jun 2020
Closed: MATLAB Answer Bot on 20 Aug 2021
How do i run though the loop and print the valid numbers in the 'value' variable?
The condition is to determine the 𝑟-value cases that achieve 𝐵 ≥ 10 at 𝑡 = 5.
colourmap = [228 26 28; 55 126 184; 77 175 74; ...
152 78 163; 255 127 0]/255;
%variables
r = [0.05 0.1 0.5 1 10];
t = (0:0.01:100);
k = 15
B0 = 1;
%
for i = 1: length(r)
dBdt =@(t,B) (r(i)).*B.*(1-(B./k));
[t,B] = ode45(dBdt,t,B0);
figure(6);
loglog(t,B,'color',colourmap(i,:))
grid on;
hold on;
if B>= 10
value = []
end
end
xlabel('Time');
ylabel('Bacteria growth');
title('B against t');
legend(num2str(r','r = %5.2f'),'location','northwest');

Answers (1)

Star Strider
Star Strider on 4 Jun 2020
Change the loop slightly, using the two new commented lines:
for i = 1: length(r)
dBdt =@(t,B) (r(i)).*B.*(1-(B./k));
[t,B] = ode45(dBdt,t,B0);
figure(6);
loglog(t,B,'color',colourmap(i,:))
grid on;
hold on;
idx = B>= 10; % Logic Index Vector
value{i,:} = [t(idx), B(idx)]; % Save To Cell Array
end
That ran without error whan I tested it just now.
  2 Comments
chong kai sheng
chong kai sheng on 4 Jun 2020
Tks for the help.
However the ans i needed to achieve was r = 1 , 10
I initally though i had to loop though to find the values that meet the condition but i ran into a wall and got stuck.
Star Strider
Star Strider on 4 Jun 2020
My pleasure.
I am not certain that I understand what you want to do.
Try this:
for i = 1: length(r)
dBdt =@(t,B) (r(i)).*B.*(1-(B./k));
[t,B] = ode45(dBdt,t,B0);
figure(6);
loglog(t,B,'color',colourmap(i,:))
grid on;
hold on;
Bmtx(:,i) = B; % Save To Matrix
end
t5idx = find(t >= 5, 1, 'first'); % Time Index
colidx = find(Bmtx(t5idx,:) >= 10, 1, 'first'); % Column Index
value = Bmtx(t5idx,colidx) % ‘value’
It retuns the column of ‘Bmtx’ that has the first instance that satisfies the conditions, and the associated ‘value’.

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!