function [media, desv] = med_desv(x)

% ESTAD Media y desviación típica. Para un vector x, 
% estad(x) da la media y la desviación típica de x. 
% Para una matriz x, estad(x) da dos vectores fila conteniendo, 
% resp., la media y la desviación típica de cada columna. 

[m n] = size(x); 
if m == 1 
    m = n; % caso de un vector fila 
end 

media = sum(x)/m
desv = sqrt(sum(x.^2)/m - media.^2)