hello
we should have freq the same for all 240 columns so we are not obliged to store 240 versions of the same vector
on the other side we can store the result amp in an 2D array - using the for loop :
for ci = 1:240
[amp(:,ci), freq]=periodogram(detrend(seis(:,ci)),[],[],Fs);
end
plot(freq,amp);