https://andymaypetrophysicist.com/r-code-3/