usendarray::{RcArray, Dimension};/// calculate right hand side (rhs) of equation of motion from current state
pubtraitEOM<D: Dimension> {fnrhs(&self, RcArray<f64, D>)->RcArray<f64, D>;}/// calculate next step by integrating the equation of motion
pubtraitTimeEvolution<D: Dimension> {fniterate(&self, RcArray<f64, D>)->RcArray<f64, D>;}