some symbolic computations are not going well
neo-hookean hyperelasticity working well