function P = ceil(P)

% CEIP -- built-in CEIL, extended to matrix polynomials
%
%        P = ceil(P)
%
% Round P.coef towards +infinity.

% Copyright (c) 2004 by Fritz Keinert (keinert@iastate.edu),
% Dept. of Mathematics, Iowa State University, Ames, IA 50011.
% This software may be freely used and distributed for non-commercial
% purposes, provided this copyright statement is preserved, and
% appropriate credit for its use is given.
%
% Last update: Feb 20, 2004

if isnumeric(P)
    P.coef = ceil(P.coef);
elseif ~issymbolic(P)
% P contains only symbolic constants, not symbolic variables
% We convert to numbers, round, and convert back to symbolic form
% The "reshape" is necessary since Matlab cannot handle
% three-dimensional symbolic arrays correctly
    P.coef = reshape(sym(ceil(double(P.coef(:)))),size(P.coef));
else
    error('CEIL is not defined for symbolic variables');
end
P = trim(P);