二项式平滑的Fortran程序

请教二项式平滑（binomialsmooth)的算法或程序！subroutinerunning2(aa,kkk,bb,key)*****binomialrunning.suittosmall'kkk'******aa:originaldata,1-diemsion,aa(key)*bb:smootheddata*kkk:lengthofrunningwindow*key:lengthofaarealaa(key),bb(key)reala(0:kkk*2+1),b(0:kkk*2+1)kw=kkk*2+1!movingwindowdok=0,kwa(k)=0.b(k)=0.enddoa(1)=1.b(1)=1.dok=1,kwk2=(k+1)/2yy=0.doj=1,kb(j)=a(j)+a(j-1)enddodoj=1,kwa(j)=b(j)b(j)=b(j)/(2**(k-1))yy=yy+b(j)enddoif(mod(k,2).ne.0)thenxx=0.dok0=1,kxx=xx+aa(k0)*b(k0)/yyenddobb(k2)=xxendifenddoyb=yydok=kkk+2,key-(kkk+1)xx=0.j=1dok0=k-kkk,k+kkkxx=xx+aa(k0)*b(j)/ybj=j+1enddobb(k)=xxenddodok=0,kwa(k)=0.b(k)=0.enddoa(1)=1.b(1)=1.dok=1,kwk2=key-k/2yy=0.doj=1,kb(j)=a(j)+a(j-1)enddodoj=1,kwa(j)=b(j)b(j)=b(j)/(2**(k-1))yy=yy+b(j)enddoif(mod(k,2).ne.0)thenxx=0.dok0=1,kxx=xx+aa(key+1-k0)*b(k0)/yyenddobb(k2)=xxprint*,k2endifenddoreturnend