complex(abs(real(q)),abs(imag(q))); end

function za=za(q), za= complex(abs(real(q)),abs(imag(q))); end
%***FIRST EXECUTABLE STATEMENT  CPEVL
if( d1==0.0 )
d1 = real(i1mach(10)).^(1-i1mach(11));
end;
np1 = fix(n + 1);
for j = 1 : np1;
ci = 0.0;
cim1 = a(j);
bi = 0.0;
bim1 = 0.0;
mini = fix(min(m+1,n+2-j));
for i = 1 : mini;
if( j~=1 )
ci = c(i);
end;
if( i~=1 )
cim1 = c(i-1);
end;
c(i) = cim1 + z.*ci;
if( kbd )
if( j~=1 )
bi = b(i);
end;
if( i~=1 )
bim1 = b(i-1);
end;
t = bi +(3..*d1+4..*d1.*d1).*za(ci);
r = real(za(z).*complex(real(t),-imag(t)));
s = imag(za(z).*t);
b(i) =(1.+8..*d1).*(bim1+d1.*za(cim1)+complex(r,s));
if( j==1 )
b(i) = 0.0;
end;
end;
end; i = fix(mini+1);
end; j = fix(np1+1);
a_shape=zeros(a_shape);a_shape(1:numel(a))=a;a=a_shape;
c_shape=zeros(c_shape);c_shape(1:numel(c))=c;c=c_shape;
b_shape=zeros(b_shape);b_shape(1:numel(b))=b;b=b_shape;
end
%*==cpevlr.f2   processed by SPAG 6.50Rc at 08:31 on 20 Jul 2006
%DECK CPEVLR