<*> <*> <*> Propagation of Error Bars <*> <*> <*>
H2init = setup(2.00, 0.01, 1)
	:  2.0000 ± 0.0100
I2init = setup(1.00, 0.01, 1)
	:  1.0000 ± 0.0100
Keq = setup(50.5, 0.1, 1)
	: 50.5000 ± 0.1000
a = (Keq - 4.)
	: 46.5000 ± 0.1000
b = (I2init + H2init)
	:  3.0000 ± 0.0141
b *= -Keq
	: -151.5000 ± 0.7746
c = (I2init * H2init)
	:  2.0000 ± 0.0224
c *= Keq
	: 101.0000 ± 1.1468
ac4 = (a*c*4.)
	: 18786.0000 ± 217.0950
b2 = (b*b)
	: 22952.2500 ± 165.9668
determinant = (b2 - ac4)
	: 4166.2500 ± 273.2676
surd = (sqrt(determinant))
	: 64.5465 ± 2.1168
nume1 = (-b + surd)
	: 216.0465 ± 2.2541
root1 = (nume1 / 2. / a)
	:  2.3231 ± 0.0247
nume2 = (-b - surd)
	: 86.9535 ± 2.2541
root2 = (nume2 / 2. / a)
	:  0.9350 ± 0.0243
I2end = (I2init - root2)
	:  0.0650 ± 0.0263

<*> <*> <*> Significant Figures <*> <*> <*>
H2init = setup(2.00, 0.01, 1)
	: 2.000
I2init = setup(1.00, 0.01, 1)
	: 1.000
Keq = setup(50.5, 0.1, 1)
	: 50.50
a = (Keq - 4.)
	: 46.50
b = (I2init + H2init)
	: 3.000
b *= -Keq
	: -151.5
c = (I2init * H2init)
	: 2.000
c *= Keq
	: 101.0
ac4 = (a*c*4.)
	: 18790
b2 = (b*b)
	: 22950
determinant = (b2 - ac4)
	: 4160
surd = (sqrt(determinant))
	: 64.5
nume1 = (-b + surd)
	: 216.0
root1 = (nume1 / 2. / a)
	: 2.323
nume2 = (-b - surd)
	: 87.0
root2 = (nume2 / 2. / a)
	: 0.935
I2end = (I2init - root2)
	: 0.065