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