generalize filter_cubic

[albop]

Currently the filter_cubic function has the following restrictions:

1. only d<=4
2. ignores terminal conditions (always natural splines)
   3.- does only cubic splines

It would be good to remove some or all of these restrictions. 1. and 2. are probably easy.