home *** CD-ROM | disk | FTP | other *** search
/ QBasic & Borland Pascal & C / Delphi5.iso / C / Samples / C-SSP.ARJ / NACPLY.C < prev    next >
Encoding:
Text File  |  1984-08-09  |  757 b   |  43 lines
  1.    nacply(n,a,b,c)
  2.  
  3.       /*function to generate a chebyshev polynomial    */
  4.       /* of degree N.                                  */
  5.  
  6.       int n;
  7.       float a[],b[],c[];
  8.  
  9.    {
  10.       int i,ij,ijk;
  11.  
  12.       if(n <= 1)
  13.       {
  14.       a[0] = 1.;
  15.       return;
  16.       }
  17.       for(i = 0; i <= n; i++)
  18.       {
  19.        a[i] = 0.;
  20.        b[i] = 0.;
  21.        c[i] = 0.;
  22.       }
  23.       b[0] = 1.;
  24.       c[0] = 1.;
  25.  
  26.       for(ijk = 2; ijk <= n; ijk++)
  27.       {
  28.        ij = ijk + 1;
  29.        a[0] = 2. * b[0];
  30.        a[1] = 2. * b[1];
  31.  
  32.        for(i = 2; i <= ij-1; i++)
  33.        {
  34.         a[i] = (2. * b[i]) - c[i-2];
  35.        }
  36.        for(i = 0; i <= ij-1; i++)
  37.        {
  38.         c[i] = b[i];
  39.         b[i] = a[i];
  40.        }
  41.       }
  42.     }
  43.