Just spent half an hour chasing an error which resulted in my code giving answers that are a fraction of an arcsecond off. Turns out that, for some reason, exponentising reals can cause loss of precision in Fortran 🤷♂️
Once again Wolfram Alpha proves its incredible usefulness.