***behavior gen_evp
  **elasticity  isotropic
    young 2e5
    poisson 0.3
  **potential gen_evp ep
   *isotropic constant
    R0 300.
 ***return