I guess you would need to check that usual range of each variable. And for each operation you perform measure the difference between the largest and smallest value. This will give you the precision of the op. Then you can try and compute the worst case scenario for your function.
I’m sure you can find some material on google from people studying these subjects for differential equation solvers.