Or, don't even write the fortran manually, just transpile the R function to fortran: https://github.com/t-kalinowski/quickr