Zero phase digital filtering

```filtfilt``` function and associated tests are included.

This code is merely a wrapper for the work done by @mbauman
This could be `reverse(x, pad_length+1, length(x)-pad_length)`, and then the second `flipud` is not necessary.