[n,m,a,x,c]=cpevlr(n,m,a,x,c);

function [n,m,a,x,c]=cpevlr(n,m,a,x,c);
persistent ci cim1 i j mini np1 ; 

if isempty(ci), ci=0; end;
if isempty(cim1), cim1=0; end;
if isempty(i), i=0; end;
if isempty(j), j=0; end;
if isempty(mini), mini=0; end;
if isempty(np1), np1=0; end;
%***BEGIN PROLOGUE  CPEVLR
%***SUBSIDIARY
%***PURPOSE  Subsidiary to CPZERO
%***LIBRARY   SLATEC
%***TYPE      SINGLE PRECISION (CPEVLR-S)
%***AUTHOR  (UNKNOWN)
%***SEE ALSO  CPZERO
%***ROUTINES CALLED  (NONE)
%***REVISION HISTORY  (YYMMDD)
%   810223  DATE WRITTEN
%   890531  Changed all specific intrinsics to generic.  (WRB)
%   890831  Modified array declarations.  (WRB)
%   891214  Prologue converted to Version 4.0 format.  (BAB)
%   900402  Added TYPE section.  (WRB)
%***end PROLOGUE  CPEVLR
a_shape=size(a);a=reshape(a,1,[]);
c_shape=size(c);c=reshape(c,1,[]);
%***FIRST EXECUTABLE STATEMENT  CPEVLR
np1 = fix(n + 1);
for j = 1 : np1;
ci = 0.0;
cim1 = a(j);
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 + x.*ci;
end; i = fix(mini+1);
end; j = fix(np1+1);
a_shape=zeros(a_shape);a_shape(:)=a(1:numel(a_shape));a=a_shape;
c_shape=zeros(c_shape);c_shape(:)=c(1:numel(c_shape));c=c_shape;
end
%DECK CPOCO