home *** CD-ROM | disk | FTP | other *** search
- lincng(a,nda,b,c)
-
- /*function to make a linear change of variables*/
- /*in a given polynomial.*/
-
- int nda;
- float a[],b[],c[];
-
- {
- int i,nz;
- float b1pwr;
-
- b1pwr = 1.0;
- nz = nda;
-
- for(i = 0;i <= nz; i++)
- c[i] = a[i];
-
- while(nz > 0)
- {
- for(i = 1; i <= nz; i++)
- c[i] += c[i-1]*b[1];
-
- c[nz] = c[nz]*b1pwr;
- b1pwr *= b[0];
- nz--;
- }
- c[0] *= b1pwr;
- }
-
-