Producing outputs from loop

I have a matrix (xMatrix: 128 x 1068). I want to run my find peak code for all 1068 columns of the matrix to calculate the peak amplitude mean, peak amplitude std dev and number of peaks for each column. How can I do this so I create a new matrix (1068 x 3) displaying all these values for each coulmn?

Here is my find peak code:

for N=(1:1068)

[pks,locs]=findpeaks(xMatrix(:,N),'MinPeakDistance',20,'MinPeakHeight',750);

B=rmoutliers(pks,'mean');

[B,TF]=rmoutliers(pks);

peakampmean=mean(B);

peakstd=std(B);

numpeaks=numel(B);

end

Thanks

Sam

Answer 1

[Erivelton Gualter]

for N=1:1068
    [pks,locs] = findpeaks(xMatrix(:,N),'MinPeakDistance',20,'MinPeakHeight',750);
    B = rmoutliers(pks,'mean');
    [B,TF] = rmoutliers(pks);
    peakampmean = mean(B);
    peakstd = std(B);
    numpeaks = numel(B);
    
    % New Matrix 
    newMatrix(1,N) = peakampmean;    
    newMatrix(2,N) = peakstd;
    newMatrix(3,N) = numpeaks;
end

[sam van Bohemen]

Thanks again for your help!!!