for loops not working

1 view (last 30 days)
sriram amarnath
sriram amarnath on 24 Jun 2019
Commented: sriram amarnath on 24 Jun 2019
The below code is to calculate the modes on a plate. The for loops are not using all the values of m & n
% Modes on a plate
Lx=0.286; % length of plate(m)
Lz=0.198; % breadth of plate(m)
for m = 1:1000 % mode numbers
for n = 1:10
func=((n.*pi)./Lx).^2 +((m.*pi)./Lz).^2;% summation function n=1:10 & m= 1:1000
end
end
  2 Comments
dpb
dpb on 24 Jun 2019
I'm sure the loop runs 10,000 times--only you don't do anything inside the loop except write over the same variable with each computed value in succession so all you get in the end is the very last value.
sriram amarnath
sriram amarnath on 24 Jun 2019
So how do i get a matrix of every value?

Sign in to comment.

Accepted Answer

madhan ravi
madhan ravi on 24 Jun 2019
Edited: madhan ravi on 24 Jun 2019
The same goal can be achieved trivially using meshgrid() or ndgrid() without a loop:
% Modes on a plate
Lx=0.286; % length of plate(m)
Lz=0.198; % breadth of plate(m)
n=10;
m=1000;
[Y,X]= ndgrid(1:m,1:n);
func=(X*pi/Lx).^2 +(Y*pi/Lz).^2;
But the modication in your code is as follows:
func = zeros(1000,10); % before loop
func(m,n) = ... %just add (m,n) next to func inside the loop
  1 Comment
sriram amarnath
sriram amarnath on 24 Jun 2019
Thanks madhan ravi your a saviour!

Sign in to comment.

More Answers (1)

dpb
dpb on 24 Jun 2019
"The MATLAB way" is to avoid the loops altogether...
N=10; M=1000;
[X,Y]=meshgrid([1:N],1:M]);
func=(X*pi/Lx).^2 +(Y*pi/Lz).^2;

Categories

Find more on Loops and Conditional Statements in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!