$fd[$n] = ( 0.25*$d[$n-2] + 0.5*$d[$n-1]
            + 1.0*$d[$n]
            + 0.5*$d[$n+1] + 0.25*$d[$n+2] )
             / (0.25+0.5+1.0+0.5+0.25);