Problem with multiple for loops
1 次查看(过去 30 天)
显示 更早的评论
for i=1:n
for j=1:m
for k=1:m
if(k~=j)
rij(i,k)=x(i,j)-x(i,k);
fij(k) = double((48 /((rij(i,k))^13))-(24 /((rij(i,k))^7)));
p=p+1;
elseif (k==j)
fij(k)=0;
else
disp('Something went wrong')
end
end
f(j)=sum(fij);
xddot(i,j)=f(j)/m;
x(i+1,j)=x(i,j)+ ((xdot(i,j))*delt);
xdot(i+1,j)= xdot(i,j) + ((xddot(i,j))*delt);
end
end
Above shown is my code and the problem is, in the innermost loop the program is not running the if command properly and just runs the elseif command and gets out of the inner for loop and never returns as it should. Anyone with a solution!
0 个评论
回答(2 个)
Jose Jeremias Caballero
2011-12-31
Observation.
Since k and j are natural numbers,
they may be identical or different.
Now if k is different from j then we do.
rij (i, k) = x (i, j)-x (i, k);
fij (k) = double ((48 / ((rij (i, k)) ^ 13)) - (24 / ((rij (i, k)) ^ 7)));
p = p +1;
contrast model is done
fij (k) = 0
So the other option is meaningless.
0 个评论
Dr. Seis
2011-12-31
So after you run the code above (and assuming you initialized "p = 0;" before the code) you get "p" still equal to 0 after the code has completed?
p should be equal to n*(m^2-m), so the only time p should equal 0 in the code above is if m = 1. Is m > 1?
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Particle & Nuclear Physics 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!