home *** CD-ROM | disk | FTP | other *** search
/ QBasic & Borland Pascal & C / Delphi5.iso / C / Samples / C-SSP.ARJ / CROSCO.C < prev    next >
Encoding:
Text File  |  1984-08-13  |  856 b   |  41 lines
  1.    crosco(a,b,n,l,bla,bda)
  2.  
  3.       /* finds the cross-covariances of series a with     */
  4.       /* series b.                                        */
  5.  
  6.       int l,n;
  7.       float a[],b[],bla[],bda[];
  8.  
  9.    {
  10.       int i,j,ij,nj;
  11.       float avea,aveb,aveab,fn,fnj,sr,ss;
  12.  
  13.       avea = 0.;
  14.       aveb = 0.;
  15.       fn = n;
  16.       for(i = 0; i <= n-1; i++)
  17.       {
  18.        avea = avea + a[i]/fn;
  19.        aveb = aveb + b[i]/fn;
  20.       }
  21.       aveab = avea * aveb;
  22.  
  23.       for(j = 1; j <= l; j++)
  24.       {
  25.        nj = n - j + 1;
  26.        fnj = nj;
  27.        fnj = 1./fnj;
  28.        sr = 0.;
  29.        ss = 0.;
  30.  
  31.        for(i = 1; i <= nj; i++)
  32.        {
  33.        ij = i + j - 1;
  34.        sr = sr + (a[i-1] * b[ij-1] * fnj);
  35.        ss = ss + (a[ij-1] * b[i-1] * fnj);
  36.        }
  37.       bla[j-1] = sr - aveab;
  38.       bda[j-1] = ss - aveab;
  39.       }
  40.    }
  41.