Error propagation

[szaghi]

If a physical quantity has an associated estimation of error, a more reliable computation should automatically propagate the errors. Possible pseudo code

type(qreal) :: q1
type(qreal) :: q2
type(qreal) :: q3

q1 = qreal(magnitude=1., unit=si%unit('second'), error=0.1)
q2 = qreal(magnitude=3., unit=si%unit('second'), error=0.2)

q3 = q1 + q2

print*, q3%stringify(with_error=.true.)
! output
! 4.0 +- 0.5 s

References

[1]rosettacode
[2]Measurements.jl
[3]uncertainties.py