function y=skewn_kurto(x,n)

%   Matlab code for computing skewness or kurtosis
%   Coded by Sungwon Shin, March 16, 2004
%   n = 3: Skewness
%   n = 4: Kurtosis
    
    if n ~=3 & n ~=4
        disp('Error !!! wrong input of n! n should be either 3 or 4 !');
        return
    end
    
    mu=mean(x);
    ic=size(x);        % size of matrix
    if ic(1)+ic(2) ==2
        y=0;
    elseif ic(1) ~= 1 & ic(2) ~=1
        for i=1:ic(2)
            Ex(i)=mean((x(:,i)-mu(i)).^n);
            sigm(i)=sqrt(mean((x(:,i)-mu(i)).^2));        
            y(i)=Ex(i)/sigm(i)^n;
        end
    else 
        Ex=mean((x-mu).^n);
        sigm=sqrt(mean((x-mu).^2));        
        y=Ex/sigm^n;
    end
    if n ==3
        msg=['Skewness === ' num2str(y)];
        disp(msg);
    elseif n ==4
        msg=['Kurtosis === ' num2str(y)];        
        disp(msg);
    end
    return