function Z = mpower(X,Y)

%^   Matrix power.

%   Z = X^y is X to the y power if y is a scalar and X is square. If y

%   is an integer greater than one, the power is computed by repeated

%   squaring. For other values of y the calculation involves

%   eigenvalues and eigenvectors.

%

%   Z = x^Y is x to the Y power if Y is a square matrix and x is a scalar.

%   Computed using eigenvalues and eigenvectors.

%

%   Z = X^Y, where both X and Y are matrices, is an error.

%

%   C = MPOWER(A,B) is called for the syntax 'A ^ B' when A or B is an

%   object.

%

%   See also POWER.

%   Copyright 1984-2018 The MathWorks, Inc.

if isscalar(X) && isscalar(Y) 

    Z = X.^Y; 

else

    if isinteger(X) || isinteger(Y)

        error(message('MATLAB:mpower:integerNotSupported'));

    end

    if ~ismatrix(X) || ~ismatrix(Y)

        error(message('MATLAB:mpower:inputsMustBe2D'));

    end

    if isa(X,'single') || isa(Y, 'single')

        try

            X = single(X);

            Y = single(Y);

        catch err

            if (strcmp(err.identifier, 'MATLAB:unimplementedSparseType'))

                error(message('MATLAB:mpower:sparseSingleNotSupported'))

            else

                rethrow(err);

            end

        end

    else

        X = double(X);

        Y = double(Y);

    end

    if isscalar(Y) && isreal(Y) && isfinite(Y) && round(Y) == Y ... 

            && size(X,1) == size(X,2)

        Z = integerMpower(X,Y);

    else

        Z = generalMpower(X,Y);

    end

end