function exped = cbexp(x,e)
% raise a sine to exponent > or < 1
s = sign(x);
exped = s .* ((s .* x) .^ e);
end
