Added mathematica scripts for semianalytical solving of 4 dipoles approximation