extern int padeapprox(int pn,int pm,unsigned int * *ps); 

int pn;
int pm;
int pnmax;            
int pmmax;
int pi,pj;

realvector ps_ram;unsigned int * *ps=ps_ram.v;
realvector pden_ram;unsigned int * *pden=pden_ram.v;
realvector pnum_ram;unsigned int * *pnum=pnum_ram.v;

int padeapprox(int pn,int pm,unsigned int * *ps)
  {zchs(z);
   for (pi=1;pi<pm+1;pi++)
      {zrecall(ps[pn+pi]);
       zchs(z);
       zstore(wr[pi]);
       for (pj=1;pj<pm+1;pj++) 
         {zrecall(ps[pn+pi-pj]);
          zstore(wm[pi][pj]);}
      }
   wgenlin(pm,wm,wr);
   zconinttomp(1);
   zstore(pden[0]);
   for (pi=1;pi<pm+1;pi++)
      {zrecall(wx[pi]); 
       zstore(pden[pi]);
      }
   zrecall(ps[0]);  
   zstore(pnum[0]);
   for (pi=1;pi<pn+1;pi++)
      {zzero(pnum[pi]);
       for (pj=0;pj<pi+1;pj++)
         {zmult(ps[pj],pden[pi-pj]);
          zadd(z,pnum[pi]);
          zstore(pnum[pi]);
         }
      }
   return 0;
  }