_1_Changed definition of the multipoles. New definition is:
U1 = (x^2-y^2)/2
U2 = (2z^2-x^2-y^2)/2
U3 = xy
U4 = yz
U5 = zx

_2_Added synthetic trap simulation data to test on a high level how the code reads/analyzes/post-processes data from BEM solver.
There are three new scripts for this:
generate_synthetic_data.m: generate synthetic data and saves them to a text file
potential_from_multipoles: used by above script
project_parameters_synth: use this one instead of project_parameters if you want to play with the synthetic data. The latter is
                          meant to facilitate switching back and forth between real simulation and synthetic simulation
_3_Fixed a bug in get_trapping_field (plottingOption was not plotting in case of a single input text file)
_4_Removed forced axis scaling in plotN