Heterogeneous elastic properties

Hi, there are several ways of doing this. You could for instance define spatially varying subdomains, or use a DG-0 function.