Problem with multiple for loops

1 次查看(过去 30 天)
Harshit
Harshit 2011-12-31
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!

回答(2 个)

Jose Jeremias Caballero
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.

Dr. Seis
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?

类别

Help CenterFile Exchange 中查找有关 Particle & Nuclear Physics 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by