int main (int argc, char *argv[])
{
  int nmax, l, lc, n, m;
  Cheby12WorkSpace *wc=0;
  double *ptr,*phi;

  if (argc <3) {
    printf("chebyshev12 requires nmax as input argument followed by phi values\n");
    return 0;
  }

  nmax=atoi(argv[1]);
  l=argc-2;
  if ((phi=(double*)calloc(l,sizeof(double)))==0) {
    printf("chebyshev12 could not allocate memory for %d phi-values\n",l);
    return 0;
  }
  if ((wc=chebyshev12_alloc(nmax,l))==0) {
    printf("chebyshev12 could not allocate memory for n=%d\n and l=%d\n",nmax,
           l);
    return 0;
  }
  for (n=0; n<l; n++) {
    phi[n]=atof(argv[n+2]);
  }
  chebyshev12(phi,wc);

  ptr=wc->t;
  for (lc=0; lc<l; lc++) {
    printf("pt %d:\n",lc);
    for (m=-nmax; m<=nmax; m++) {
      printf("T[%2d](%7.4f)=%7.4f\n",m,phi[lc],*ptr++);
    }
  }
  chebyshev12_free(wc);
  free(phi);
  return 1;
}
static void allocMTXChdata (MTXCh *x) 
{
   x->phi=(double*)calloc(x->l,sizeof(double));
   x->wc=chebyshev12_alloc(x->nmax,x->l);
   x->list_ch=(t_atom*)calloc(x->l*(2*x->nmax+1)+2,sizeof(t_atom));
}