Further changes in formula. Forces (in particular the magnetic ones) must be checked. WARNING: the constant "intensity" has now a random value, used for debugging.