n_x = 0;
n_y = 0;
n_z = 1;
L = 2.0;
B = 1.0;
H = 0.5;
N=2;
c = 343;
for n_x = 0:N
for n_y = 0:N
for n_z = 0: N
f = (c/2) * sqrt( ((n_x/L)^2) + ((n_y/B)^2) + ((n_z/H)^2) ) ;
%disp(f)
formatSpec = 'N_X is %4.2f, N_Y is %4.2f, N_Z is %4.2f, Frequency is %8.3f mm\n';
fprintf(formatSpec,n_x,n_y,n_z,f)
end
end
%MODE( 1, n_x+n_y+n_z) = f
end
