Added files
Functional prototype
Apparently working prototype with tests
Symmetrical mode.

- Added a symmetrical misfit mode. This runs the lam2 misfit process with the two
  operators in both directions, and takes the average misfit of the two.
Initiated a splitting misfit example.
First stab at effective splitting using wavelet

Use [fast tlag] = MS_effective_splitting_N(0.1, ...
 55.0, [55.0 10.0], [1.0 1.0], 'mode', 'gaussianwavelet')
to do effective splitting calcualtion using this new
method.
Add tests, unwind output

Make the output consistent with the S&S method. Add tests
following some of those from S&S. Note that in three cases
(tests commented out) the result is a few degrees different.
Shows up more for longer periods. Is this close to a faliure
mode for S&S or is something wrong with the wavelets?
Remove some code duplication from new functions

And tidy up the comments.
Move splitting measurment into seperate function