Using for loop for midpoint method
2 次查看(过去 30 天)
显示 更早的评论
Hi guys, this is my code
m = 1000;
f = 75000;
v0 = 0;
h = 0.1;
tspan = [0 10];
c = [50,100,500,1000,2000];
dvdt = cell(1,length(c));
for i = 1:length(c);
dvdt{i} = @(t,v) (f-c{i}*v)/m;
end
t =cell(1,length(c));
v =cell(1,length(c));
for k = 1:length(c)
[t{k},v{k}]=midpoint(dvdt{1,k},tspan,v0,h)
end
%% Question
I want to use midpoint to calculate the t and v, and my dvdt is relied on c(constant values run from 50 to 2000). The code works, but t and v are released just only based on c = 50, which means my for loop does not work. Can someone please help !!!
Thanks in advance
0 个评论
回答(0 个)
另请参阅
类别
在 Help Center 和 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!