For instance: In (a+-b) * (c+-d), you're multiplying the intervals (a-b,a+b) and (c-d,c+d). If all those endpoints are positive, then you get the interval ((a-b) * (c-d),(a+b) * (c+d)).
The answer is then: ((a-b) * (c-d)+(a+b) * (c+d))/2 +- ((a+b) * (c+d) - (a-b) * (c-d))/2
If some of those endpoints have different signs, then you need to break it down into cases.
I can see that ((a-b) (c-d)+(a+b) * (c+d))/2 is different to a*b, and that makes sense when you think about it, but the second part is confusing me a bit.
Also, I'm a bit confused by your suggestion to break it down into cases, do you mean if an uncertainty is negative? Because at the moment I'm making all uncertainties positive, but if there are cases where negative ones are correct, I need to fix that up.
You also need to clear up your own confusion about what exactly you mean by "uncertainty". You seem to go back and forth between several different measures.
Uncertainties aren't negative. But solve the rest of the issues, and what I meant by cases should become more clear.