Conversion to gsl shell of 1D root determination functions
Simplifiy wave particle demo using new matrix syntax