auto differentiation for up to 4th derivatives

Build

1. Download files
2. use make run to build and run the tests

Tested with gfortran 10.1.0