From energycons follows the equation of motion (for details
see [1]}, [2]}) for the phase deviation of unit \(i\)
\(\ddot{\phi }_i = P_i - \kappa \dot{\phi }_i +\sum _j K_{ji}^{\rm MAX} \sin (\phi _j - \phi _i),\) 
