extern int cpadeapprox(int cpn,int cpm,unsigned int * *cps); 

int cpn;
int cpm;
int cpnmax;            
int cpmmax;
int cpi,cpj;

complexvector cps_ram;unsigned int * * *cps=cps_ram.cv;
complexvector cpden_ram;unsigned int * * *cpden=cpden_ram.cv;
complexvector cpnum_ram;unsigned int * * *cpnum=cpnum_ram.cv;

int cpadeapprox(int cpn,int cpm,unsigned int * * *cps)
  {zchs(z); 
   for (cpi=1;cpi<cpm+1;cpi++) 
      {crecall(cps[cpn+cpi]);
       cchs(cz);
       cstore(cwr[cpi]);
       for (cpj=1;cpj<cpm+1;cpj++) 
         {crecall(cps[cpn+cpi-cpj]);
          cstore(cwm[cpi][cpj]);}
      }
   csolvlin(cpm,cwm,cwr);
   cconinttomp(1,0);
   cstore(cpden[0]);
   for (cpi=1;cpi<cpm+1;cpi++)
      {crecall(cwx[cpi]); 
       cstore(cpden[cpi]);
      }
   crecall(cps[0]);  
   cstore(cpnum[0]);
   for (cpi=1;cpi<cpn+1;cpi++)
      {czero(cpnum[cpi]);
       for (cpj=0;cpj<cpi+1;cpj++)         
         {cmult(cps[cpj],cpden[cpi-cpj]);
          cadd(cz,cpnum[cpi]);
          cstore(cpnum[cpi]);
         }
      }
   return 0;
  }