I finally got around to putting my errorbar calculation module on my website. <http://users.bigpond.net.au/gazzar/python.html> This should be of interest to Numeric Python users who deal with calculation of error-bounds on experimental data. Synopsis: A module providing a Python number type class and helper functions to ease the task of computing error bounds on experimental data values. This module defines an abstract data type, Err, which carries a central/prime value and upper and lower error bounds. Helper functions are provided to allow Python Numeric rank-1 arrays of Err objects to be constructed with ease and to apply Numeric Ufuncs. Written under Python 2.3.3 and Numeric 23.0 Example of usage: upperPressure = ArrayOfErr([909., 802., 677., 585., 560., 548.], 1.0) lowerPressure = ArrayOfErr([144., 246., 378., 469., 493., 505.], 1.0) pressureDiff = upperPressure - lowerPressure V_RStandard = ArrayOfErr([2.016, 2.016, 2.020, 2.017, 2.021, 2.019], 0.001) R = 100.2 # standard resistor value [Ohm] I = V_RStandard / R # current [A] V_RAB = ArrayOfErr([6.167, 6.168, 6.170, 6.160, (6.153, 0.02), (5.894, 0.01)], 0.002) R_AB = V_RAB / I logR_AB = ApplyUfuncToErr(R_AB, log10) # This means log10(R_AB) print PrimeVals(logR_AB) print MinVals(logR_AB) print MaxVals(logR_AB) Enjoy, Gary Ruben -- ___________________________________________________________ Sign-up for Ads Free at Mail.com http://promo.mail.com/adsfreejump.htm