function [psi,X] = cpauwavf(varargin)
% cpauwavf Complex Paul wavelet.
%   [PSI,X] = CPAUWAVF(LB,UB,N,P) returns values of the Paul wavelet

% Check arguments.
if errargn(mfilename,nargin,[0,3 4],nargout,[0:2]), error('*'); end
if nargin==0
   if ~exist('maple') , nmax = 8; else , nmax = 45; end
   psi = nmax;
   % psi contains the number max for Gaussian wavelets.
   % This number depends of Symbolic Toolbox
   return
end
if nargin==4
   if ischar(varargin{4})
       [fam,num] = wavemngr('fam_num',varargin{4});
       num = wstr2num(num);
   else
       num = varargin{4};
   end
else
   num = 1;
end
if errargt(mfilename,num,'int') , error('*'); end

% Compute values of the Gauss wavelet.
X = linspace(varargin{1:3});  % wavelet support.


m=num;
eta = X;
a = (2.^m .* i.^m .* gamma(m+1));
b = (1 - i.*eta).^(-m-1);
c = sqrt(pi .* gamma(2.*m));

psi = (a .* b) ./ c;